ex5_24.m

源代码集合

x0=[1.2; 0; 0; -1.04935751];
tic, [t,y]=ode45('apolloeq',[0,20],x0); toc
length(t), plot(y(:,1),y(:,3))

tic, [t1,y1]=rk_4('apolloeq',0,0.001,20,x0); toc
length(t1), plot(y(:,1),y(:,3)),