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)=W11(X(1,1), X(2,1), X(3,1));
  Z(1,2)=W12(X(1,1), X(2,1), X(3,1));
  Z(1,3)=W13(X(1,1), X(2,1), X(3,1));

  Z(2,1)=W21(X(1,1), X(2,1), X(3,1));
  Z(2,2)=W22(X(1,1), X(2,1), X(3,1));
  Z(2,3)=W23(X(1,1), X(2,1), X(3,1));

  Z(3,1)=W31(X(1,1), X(2,1), X(3,1));
  Z(3,2)=W32(X(1,1), X(2,1), X(3,1));
  Z(3,3)=W33(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)));