ex5_24.m

来自「源代码集合」· M 代码 · 共 7 行

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)),