我按附件中的公式中的内容编了一个程序，对于两次的求和用matlab不知道应该怎样写，程序可以运行，但是得到的图却不对，最终得到的图应该是峰变窄，面积不变，大家帮我看看应该怎么修改，在下感激不尽！程序如下：
clear
clc
x(1)=40;x(2)=60;
m=2;
x=1:0.1:100;
f=50*exp(-(x-40).^2/100)+70*exp(-(x-60).^2/100);

N1=zeros(size(x));
N2=zeros(size(x));
N3=zeros(size(x));
N4=zeros(size(x));
t=zeros(size(x));

N=length(x);
for k=1:1:m
    t=x-x(k)
    for j=1:N
        if t(j)>=0 & t(j)<1;
            N1(k,j)=1/2*(t(j)).^2
        elseif t(j)>=1 & t(j)<2
            N2(k,j)=3/4-(t(j)-3/2).^2
        elseif t(j)>=2 & t(j)<=3
            N3(k,j)=1/2*(3-t(j)).^2
        else
            N4(k,j)=0
        end
    end
    N=N1+N2+N3+N4;
    F=N.*f;
end
plot(F)

在线等高人啊，急啊

2次微分应该有两个定义域的限制？
问题应该在
k=1:1:m

谢谢，还有就是公式里两个求和的运算我不知道怎么写进去，你能指点我一下吗？