exm051442_1.m

这是本人收集的一些程序源代码

%exm051442_1.m
G=6.672e-11;ME=5.97e24;vy0=4000;x0=-4.2e7;t0=0;tf=60*60*24*9;
tspan=[t0,tf];Y0=[x0;0;0;vy0];
% 以下四条指令，为相平面图预置一个坐标轴范围
clf,shg,set(gca,'xlim',[-5 25]*1e7,'ylim',[-3 3]*1e3);%设置适当坐标范围		<5>
box on			%框形坐标
hold on;		%使中间结果绘在同一幅图上
ode45('DYDt4',[],[],[],G,ME,tspan,Y0);hold off