function res = sys1(t,x)
res = [x(2)^2-x(1)^2;
       0.5-x(1)^2-x(2)^2];