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