求助
浏览:85991 回答:1
function dydt=funode2(t,y)
dydt=[y(2)*y(3);-y(1)*y(3);-0.51*y(1)*y(2)]
y(1)=0;y(2)=1;y(3)=1;
y=[y(1);y(2);y(3))]
[t,y]=ode45('funode2',[0,12],y)
当y(1)从正值变成负值时,停止计算,输出y(1)为负时,那组y值
dydt=[y(2)*y(3);-y(1)*y(3);-0.51*y(1)*y(2)]
y(1)=0;y(2)=1;y(3)=1;
y=[y(1);y(2);y(3))]
[t,y]=ode45('funode2',[0,12],y)
当y(1)从正值变成负值时,停止计算,输出y(1)为负时,那组y值