function xut=NR(f,x,tol)
h=0.0001;
fx=feval(f,x);
fprimx=(feval(f,x+h)-feval(f,x))/h;
format long
disp('        x_n                f(x_n)            f/fprim')
disp([x fx fx/fprimx])
while abs(fx/fprimx)>tol
x=x-fx/fprimx;
fx=feval(f,x);
fprimx=(feval(f,x+h)-feval(f,x))/h;
disp([x fx fx/fprimx])
end
format short
