x0 = [2 0]';
tout = 0:0.1:5;
[t x] = ode45(@derivs,tout,x0);
plot(t,x(:,1),'m-',t,x(:,2),'g-') ; grid
xlabel('t')
ylabel('x1 and x2')
legend('x1','x2')
