clc;clear;
E = 1e-4; kmax = 40; X0 = [2 2 2]';
X = zeros(3,1); J = zeros(3,3);
k = 0; Fn = ones(3,1)*2*E;
J(1,2) = -2;
J(2,1) = 1;
J(2,3) = 0;
J(3,1) = 1;
J(3,2) = 1;
J(3,3) = 1;
while abs(max(Fn)) >= E & k <= kmax
    k = k+1;
    if k == 1
        Xp = X0;
    else
        Xp = X(:,k-1);
    end
    Fn(1) = exp(Xp(1)/2)-2*Xp(2)+3*Xp(3)^2-24.65;
    Fn(2) = Xp(1)-log(Xp(2))-0.699;
    Fn(3) = Xp(1)+Xp(2)+Xp(3)-6;
    J(1,1) = 0.5*exp(Xp(1)/2);
    J(1,3) = 6*Xp(3);
    J(2,2) = -Xp(2)^(-1);
    X(:,k) = Xp-inv(J)*Fn;
end
disp('x1');disp(X(1,k));
disp('x2');disp(X(2,k));
disp('x3');disp(X(3,k));
disp('liczba iteracji:');disp(k);
disp('dokładność ');disp(abs(max(Fn)));