Решение систем нелинейных уравнений. Метод Ньютона

clc
clear all
xn(1)=0; xn(2)=0; xn(3)=0; xn(4)=0;
x(1)=-10; x(2)=-10; x(3)=-10; x(4)=-10;
 
eps=1e-6;
fl_stop=true;
while  fl_stop
    xn(1)=sqrt(-x(2)-7*x(3)-3*x(4))/2; %x1 из 1 ур
    xn(3)=-sqrt(-x(1)-3*x(2)-6)/10; %x3 из 2 ур
    xn(2)=(-12*x(1)-9*x(3)^2-x(4))/2; %x2 из 3 ур
    xn(4)=-sqrt(-3*x(1)-5*x(2)-3*x(3)); %x4 из 4 ур
    if sum(abs(x-xn))<eps
        fl_stop=false;
    end
    x=xn;
end
x

clc
clear all
x(1)=-10; x(2)=-10; x(3)=-10; x(4)=-10;
 
eps=1e-6;
fl_stop=true;
in=1000;
while  fl_stop && in>0
    xn(1)=-sqrt((-x(2)-7*x(3)-3*x(4))/2); %x1 из 1 ур
    %xn(2)=-2*x(1)^2-7*x(3)-3*x(4); %x2 из 1 ур
    %xn(3)=(-2*x(1)^2-x(2)-3*x(4))/7; %x3 из 1 ур
    %xn(4)=(-2*x(1)^2-x(2)-7*x(3))/3; %x4 из 1 ур
 
 
    %xn(1)=-3*x(2)-10*x(3)^2-6; %x1 из 2 ур
    xn(2)=(-x(1)-10*x(3)^2-6)/3; %x2 из 2 ур
    %xn(3)=-sqrt((-x(1)-3*x(2)-6)/10); %x3 из 2 ур
 
    %xn(1)=(-2*x(2)-9*x(3)^2-x(4))/12; %x1 из 3 ур
    %xn(2)=-12*x(1)-2*x(2)-9*x(3)^2; %x2 из 3 ур
    xn(3)=-sqrt((-12*x(1)-2*x(2)-x(4))/9); %x3 из 3 ур
    %xn(4)=-12*x(1)-2*x(2)-9*x(3)^2; %x4 из 3 ур
 
    %xn(1)=(-5*x(2)-3*x(3)-x(4)^2)/3; %x1 из 4 ур
    %xn(2)=(-3*x(1)-3*x(3)-x(4)^2)/5; %x2 из 4 ур
    %xn(3)=(-3*x(1)-5*x(2)-x(4)^2)/3; %x3 из 4 ур
    xn(4)=-sqrt(-3*x(1)-5*x(2)-3*x(3)); %x4 из 4 ур
    if sum(abs(x-xn))<eps
        fl_stop=false;
    end
    x=xn;
    in=in-1;
end
x
a1=2*x(1)*x(1)+x(2)+7*x(3)+3*x(4)
a2=x(1)+3*x(2)+10*x(3)*x(3)+6
a3=12*x(1)+2*x(2)+9*x(3)*x(3)+x(4)
a4=3*x(1)+5*x(2)+3*x(3)+x(4)*x(4)