examp7_11.m

来自「东北大学 薛定谔的高等应用数学问题MATLAB求解的源代码」· M 代码 · 共 12 行

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

options=odeset; options.RelTol=1e-6;
tic, [t1,y1]=ode45('apolloeq',[0,20],x0,options); toc
length(t1), plot(y1(:,1),y1(:,3)),

min(diff(t1))

plot(t1(1:end-1),diff(t1))