function [f J]=GaNew(x)
 f(1)=exp(x(1)-2*x(2))-0.5;
 f(2)=exp(x(1)-1.0*x(2))-1.0;
 f(3)=exp(x(1))-2.0;
 f(4)=exp(x(1)+x(2))-4.0;
 if nargout<2, return, end
 J(1,1)=exp(x(1)-2.0*x(2)); J(1,2)=-2*exp(x(1)-2*x(2));
 J(2,1)=exp(x(1)-x(2));     J(2,2)=-exp(x(1)-x(2));
 J(3,1)=exp(x(1));          J(3,2)=0;
 J(4,1)=exp(x(1)+x(2));     J(4,2)=exp(x(1)+x(2));
end