DF = @(X) [15 2*X(2) -1; -1 30 -2*X(3); -2*X(1) 1 100];
F = @(X) [15*X(1)+X(2).^2-X(3)-30;-X(1)+30*X(2)-X(3).^2-30;-X(1).^2+X(2)+100*X(3)-20];

X0 = [2 2 2]';
tol = 1e-7;
err = 1;
while err > tol
    X = X0 - DF(X0)\F(X0);
    err = norm(X0-X);
    X0 = X;
end