求求了！不知道后面哪里出了问题matlab

clear all;
close all;
syms T;
syms k;
syms k1;
syms k2;
syms i;
A=3;
H=40;
syms t;
y1=t*(2*pi)^(-0.5)*exp(-(t-10)^2/2);
y2=(2*pi)^(-0.5)*exp(-(t-10)^2/2);
fy1=int(y1,k*T,k*T+T);
fy2=int(y2,k*T,k*T+T);
a=ceil((k*T-T)/0.5)*A/T;
b=ceil(k*T/0.5)*A;
c=ceil(k*T/0.5)*A/T;
for T=0.05:0.05:5
    k2=ceil(H*0.5/(A*T));
    if(T>=0.5)
        k1=1;
        i=round(T/0.05);  % 此行增加round函数，保证i为整数
        E(i)=symsum((a*(fy1-k*T)+b*fy2+c*(fy1-k*T)),k,k1,k2)+0+40*int((2*pi)^(-0.5)*exp(-(t-1)^2/2),k*T,k*T+T);
        disp(E(i));
    else
        k1=ceil(0.5/T);
        i=round(T/0.05);  % 此行增加round函数，保证i为整数
        E(i)=symsum((a*(fy1-k*T)+b*fy2+c*(fy1-k*T)),k,k1,k2)+0+40*int((2*pi)^(-0.5)*exp(-(t-1)^2/2),k*T,k*T+T);
        disp(E(i));
    end;
end
T = round(0.05:0.05:5); % 绘图x轴修改
plot( T,double(E) )     % 绘图S转换为double量，符号量不能绘图
xlim([0.05 5])