rigid.m

来自「用来求解线性或者非线性数学方程 并在坐标系下观察其取值.」· M 代码 · 共 6 行

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