求高手指导 谢谢 我的小程序运行不出来

yuexiaoming

程序如下 m文件function dy = fun1(t,y,flag,u)

if y(1)>=1
    a=0.1*cos(u*t)-y(1)+1;
elseif y(1)<1&y(1)>-1
    a=0.1*cos(u*t);
elseif y(1)<=-1
    a=0.1*cos(u*t)-y(1)-1;
end
dy=zeros(2,1);
dy=[y(2);a];

主程序：u=[0.5, 2.0];
for j=1:2
[t,y] = ode23(@fun1,[0 5] , [2 0],[],u(j));
figure
plot(y(:,1),y(:,2));
end