求助 用ode45编的小程序有问题
function dy=rigid(t,y)dy=zeros(3,1); % a column vector
dy(1)=y(2)*y(3);
dy(2)=-y(1)*y(3);
dy(3)=-0.51*y(1)*y(2);
options=odeset('RelTol',1e-4,'AbsTol',);
=ode(@rigid,,,options);
plot(T,y(:,1),'-',T,y(:,2),'-.',T,y(:,3),'.')
以上程序运行时总是出现错误， Input argument 'y' is undefined.
请高手指点，谢谢:handshake 分成两部分即可.
%%%%% 存为函数 %%%%%%%
function dy=rigid(t,y)
dy=zeros(3,1); % a column vector
dy(1)=y(2)*y(3);
dy(2)=-y(1)*y(3);
dy(3)=-0.51*y(1)*y(2);
%%%%%% 在命令窗口运行下列代码 %%%%%%%
options=odeset('RelTol',1e-4,'AbsTol',);
=ode(@rigid,,,options);
plot(T,y(:,1),'-',T,y(:,2),'-.',T,y(:,3),'.')
%%%%%%%%%%%%%%%%%%
建议认真看看书.
回复 #2 xjzuo 的帖子
运行出来了谢谢指点，看来我真的要认真看看书:loveliness: :loveliness:
页:
[1]