function J = Jfun(x)
J = [2*x(1),              2*x(2);...
     cos(pi*x(1)/2)*pi/2, 3*x(2)^2];
