program2_1b.m

此为matlab的代码。包括数值o45算法及o23算法

%MATLAB Program2_1b
x0=[0;0;0];
t0=0;tf=8;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')