examp7_12.m

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

x0=[1.2; 0; 0; -1.04935751];
tic, [t1,y1]=rk_4('apolloeq',[0,20,0.01],x0); toc
plot(y1(:,1),y1(:,3))     % 绘制出轨迹曲线

tic, [t2,y2]=rk_4('apolloeq',[0,20,0.001],x0); toc
plot(y2(:,1),y2(:,3))     % 绘制出轨迹曲线