function dx = derivs(t,x)
dx(1) = -2*x(1).^2 + 2*x(1) + x(2) -1;
dx(2) = -x(1) - 3*x(2).^2 + 2*x(2) + 2;
dx = dx';
