clear all 
close all

x0  = 0.5;
f  = @(x) log(x);
df = @(x) 1./x;

[x k] = NewtonIteration(x0,f,df,1.0e-03,1000);

fprintf('Ergebnis = %.2e, res = %.2e nach %d Iterationen\n',x,abs(f(x)),k);

