function J = nljac(x)
    J(1,:) = [-exp(x(2)-x(1)), exp(x(2)-x(1)), 0];
    J(2,:) = [x(2), x(1), 1];
    J(3,:) = [2*x(1), x(3)-1, x(2)];
end