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