作业帮 > 综合 > 作业

matlab求解常微分方程组数值解,总出现这个,急死了

来源:学生作业帮 编辑:神马作文网作业帮 分类:综合作业 时间:2024/11/11 09:31:28
matlab求解常微分方程组数值解,总出现这个,急死了
我输入function dy=rigid(t,y) ;
dy=zeros(3,1);
dy(1)=y(2)*y(3);
dy(2)=-y(1)*y(3);
dy(3)=-0.51*y(1)*y(2);
运行[T,Y]=ode45('rigid',[0 12],[0 1 1]);后总出现
Error using ==> feval
Undefined function or method 'rigid' for input arguments of type 'double'.
Error in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.
Error in ==> ode45 at 173
[neq,tspan,ntspan,next,t0,tfinal,tdir,y0,f0,odeArgs,odeFcn,...
在0 1 1 中间加上分号吗?我试过了,还是不行啊
matlab求解常微分方程组数值解,总出现这个,急死了
[T,Y]=ode45('rigid',[0 12],[0;1;1]); %标点符号状态
plot(T,Y)
----------------------
function dy=rigid(t,y)
dy=zeros(3,1);
dy(1)=y(2)*y(3);
dy(2)=-y(1)*y(3);
dy(3)=-0.51*y(1)*y(2);