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