program2_1b.asv

来自「此为matlab的代码。包括数值o45算法及o23算法」· ASV 代码 · 共 10 行

%MATLAB Program2_1b
x0=[0;0;0];
t0=0;tf=4;tspan=[t0,tf]
[t,x]=ode45('ode_example',tspan,x0);
y=(3333333/10000)*x(:,2)|+(666667/10000*x(:,3);
subplot(1,2,1)
plot(t,x(:,1),'k',t,x(:,2),'-.r',t,x(:,3),':b')

subplot(1,2,2)
plot(t,y,'k')