clear all, clc;

F=@(X) [(X(1)-X(2))^2+X(2)^3-2;... 
X(1)^3+X(2)^4-2*X(1)*X(2)*X(3);...
X(1)^2+X(2)^2+X(3)^2-6];

J=@(X) [2*(X(1)-X(2)) -2*(X(1)-X(2))+3*X(2)^2 0;...
    3*(X(1))^2-2*X(2)*X(3) 4*(X(2))^3-2*X(1)*X(3) -2*X(1)*X(2);...
    2*X(1) 2*X(2) 2*X(3)];

X0=[1 1 -1]';
delta=inf;
while norm(delta)>1e-5
    F0=F(X0);
    J0=J(X0);
    delta=J0\F0;
    X0=X0-delta
end
F(X0)
