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