Question

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 中间加上分号吗？我试过了，还是不行啊

Answer 1

[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);