function J = J3(x)
   N = length(x);
   J = zeros(N,N);
   J(1,1) = 2*x(1);
   J(1,2) = -1;
   J(2,1) = -3*sin(x(1));
   J(2,2) = -1;
end