[t, x] = ode45('test2', [0 10], [0 1]);
plot(t, x)
title('The solution to x1'' = x2 - x1^2 and x2'' = -x1 -2*x1*x2')
xlabel('t')
ylabel('x')
grid
fprintf('Press <Enter> to continue\n')
pause
plot(x(:,1), x(:,2))
pause
plot3(x(:,1), x(:,2), t)
grid



