function dy=dong22(t,y)
m=1;
k=2;
g= %注意输入g l0
l0=
dy=[y(2);
y(1)*y(4)^2+g*cos(y(3))-(k/m)*(y(1)-l0);
y(4);
-(2/y(1))*y(2)*y(4)-(g/y(1))*sin(y(3))];
然后和一次微分方程的一样。
