function J=Joct12(x)
J=[2*x(1)+x(2)^3 3*x(1)*x(2)^2 ; 6*x(1)*x(2) 3*(x(1)+x(2))*(x(1)-x(2))];