fhuabai.m

来自「利用matlab进行动力学仿真之一」· M 代码 · 共 10 行

function ydot=fhuabai(t,y)
global a w
x1=[a+1 cos(y(2))
   cos(y(2)) 1];
x2=[y(4)*y(4)*sin(y(2))
   -w*w*sin(y(2))];
ydot=[y(3)
   y(4)
   inv(x1)*x2];