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