I1=inline("2*x*y + 2*x*z + 2*y*z - 22");
I2=inline("x^2 - 3*y*z - 2*x*y^2 + 25");
I3=inline("-3*y^3 + 2*x^2*z^2 + 6");
blad=0.5;
niedokladnosc=blad+1;
Z11=inline("2*y + 2*z");
Z21=inline("2*x - 2*y^2");
Z31=inline("4*x * z^2");
Z12=inline("2*x + 2*z");
Z22=inline("-4*x*y - 3*z");
Z32=inline("-9*y^2");
Z13=inline("2*x + 2*y");
Z23=inline("-3*y");
Z33=inline("4*x^2 * z");
X=[2; 2.5; 3];
I=zeros(3,1);
I=zeros(3,3);
while niedokladnosc>blad
Z(1,1)=Z11(X(1,1), X(2,1), X(3,1));
Z(1,2)=Z12(X(1,1), X(2,1), X(3,1));
Z(1,3)=Z13(X(1,1), X(2,1), X(3,1));
Z(2,1)=Z21(X(1,1), X(2,1), X(3,1));
Z(2,2)=Z22(X(1,1), X(2,1), X(3,1));
Z(2,3)=Z23(X(1,1), X(2,1), X(3,1));
Z(3,1)=Z31(X(1,1), X(2,1), X(3,1));
Z(3,2)=Z32(X(1,1), X(2,1), X(3,1));
Z(3,3)=Z33(X(1,1), X(2,1), X(3,1));
I(1,1)=I1(X(1,1), X(2,1), X(3,1));
I(2,1)=I2(X(1,1), X(2,1), X(3,1));
I(3,1)=I3(X(1,1), X(2,1), X(3,1));
I=I*(-1);
dX=inv(Z)*I;
X=X+dX;
disp(X);
niedokladnosc=abs(I1(X(1,1), X(2,1), X(3,1)))+abs(I2(X(1,1), X(2,1), X(3,1)))+abs(I3(X(1,1),
X(2,1), X(3,1)))
disp("***");
end
disp("wynik")
disp(X);
disp("Sprawdzenie")
disp(I1(X(1,1), X(2,1), X(3,1)));
disp(I2(X(1,1), X(2,1), X(3,1)));
disp(I3(X(1,1), X(2,1), X(3,1)));