clear all; close all; tspan=[0 30]; tspan2=0:0.2:30; y0=[1 0 0]; A=2; [t,y]=ode45('rhs',tspan2,y0,[],A); plot(t,y(:,1),t,y(:,2),'r',t,y(:,3),'g')