program2_2a.m

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

%MATLAB Program2_2a
x0=[1;0;0];
t0=0;tf=6;tspan=[t0,tf]
[t,x]=ode45('ode_example2_2a',tspan,x0);
y1=x(:,2)-x(:,3)
y2=x(:,1)+2*x(:,2)+x(:,3)
subplot(1,3,1)
plot(t,x(:,1),'k',t,x(:,2),'-.r',t,x(:,3),':b')
subplot(1,3,2)
plot(t,y1,'k')
subplot(1,3,3)
plot(t,y2,'k')