function xut=fixit(F,x,tol)
n=0;
format long
disp('                        x_n                        x_n-F(x_n)')
while abs(x-feval(F,x))>tol
n=n+1;
x=feval(F,x);
disp(n) , disp([x x-feval(F,x)])
end
format short
