function z = Hessian(x,y,r)
z=[jacobian(-2*x-2*y-2-r*((-4*x^3/(4-x^4-y^2))-(1/x^2)),x);jacobian(-2*x-2*y-2-r*((-4*x^3/(4-x^4-y^2))-(1/x^2)),y);jacobian(-2*x-2*y+1-r*((-2*y/(4-x^4-y^2))-(1/y^2)),x);jacobian(-2*x-2*y+1-r*((-2*y/(4-x^4-y^2))-(1/y^2)),x)];
