function [x k]=NewtonIteration(x0,f,df,TOL,ItMax)

x   = x0;
Res = abs(f(x));
k   = 0;
while Res>TOL & k<ItMax
   
    x   = x - f(x)/df(x);
    Res = abs(f(x));
    k   = k+1;
    %fprintf('x_%d = %.8f mit Res = %.2e\n',k,x,Res);
        
end


end