clc
x=[4;1];
fx=F(x);
k=0;
while norm(fx)>=10^(-2) && k<=250
    A=J(x);
    b=-fx;
    s=A\b;
    fprintf('%1.0f \t %2.5f \t %2.5f \t %2.5f \t %2.5f \n',k,x(1),x(2),fx(1),fx(2))
    x=x+s;
    k=k+1;
    fx=F(x);
    
end
    fprintf('%1.0f \t %2.5f \t %2.5f \t %2.5f \t %2.5f \n',k,x(1),x(2),fx(1),fx(2))