function Jv=Jfun(x)
  Jv(1,1)=2*x(1);
  Jv(1,2)=2*x(2);
  Jv(2,1)=pi/2.0*cos(pi/2.0*x(1));
  Jv(2,2)=3*x(2)^2;
endfunction