求高手指导 谢谢 我的小程序运行不出来
程序如下 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=;
主程序：u=;
for j=1:2
= ode23(@fun1, , ,[],u(j));
figure
plot(y(:,1),y(:,2));
end
页:
[1]