如图 g和k是常数 最后要求y对z作图 t是参数
k=2;g=3
myode=@(t,x)[x(2)
2*x(4)
x(4)
-3-2*x(2)];
[t,S]=ode45(myode,[0,20],[0 0 1 2]);
plot(t,S),figure
plot(S(:,1),S(:,3))%S(:,1)=y,S(:,3)=z