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 ==> fevalUndefi

来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/24 08:57:46
matlab求解常微分方程组数值解,总出现这个,急死了我输入functiondy=rigid(t,y);dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(3);d

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 ==> fevalUndefi
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求解常微分方程组数值解,总出现这个,急死了我输入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 ==> fevalUndefi
[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);