用MATLAB进行星球轨迹仿真程序

dt=0.01; x=[100 -1 50 2]'; y=x; t=0; for i=1:100000 t=[t i*dt]; x=x+dt*fn(x)'; y=[y x]; end figure plot(y(1,:),y(3,:)); figure plot(t,y(2,:)); figure plot(t,y(4,:));