求助 用ode45编的小程序有问题

lily1218

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',[1e-4 1e-4 1e-5]);
[t,y]=ode(@rigid,[0 12],[0 1 1],options);
plot(T,y(:,1),'-',T,y(:,2),'-.',T,y(:,3),'.')
以上程序运行时总是出现错误， Input argument 'y' is undefined.
请高手指点，谢谢:handshake

xjzuo

分成两部分即可.
%%%%% 存为函数 %%%%%%%
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',[1e-4 1e-4 1e-5]);
[t,y]=ode(@rigid,[0 12],[0 1 1],options);
plot(T,y(:,1),'-',T,y(:,2),'-.',T,y(:,3),'.')
%%%%%%%%%%%%%%%%%%
建议认真看看书.

lily1218

运行出来了
谢谢指点，看来我真的要认真看看书:loveliness: :loveliness: