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