j.m

用牛顿法解方程组

function A=J(X)
A=[5 3 1 1;X(2) X(1)+X(2) X(2)+X(4) 0;2*X(1) 2*X(2) 2*X(3) -2*X(4);3*X(1)^2+X(3) X(4) X(1) X(2)];
end