# 
restart;
f:=x->exp(x)-2+x;
D(f);
g:=x->x-f(x)/D(f)(x);
x0:=1.0;
x1:=g(x0);
plot(f(x),x=0..2);
line:=y-f(x0)=D(f)(x0)*(x-x0);
a:=solve(line,y);
plot([f(x),a],x=0..2);
x1:=g(x0);
x2:=g(x1);
x3:=g(x2);
x4:=g(x3);
x5:=g(x4);
Digits:=60;
x0:=1.0;
x1:=g(x0);
x2:=g(x1);
x3:=g(x2);
x4:=g(x3);
x5:=g(x4);
x6:=g(x5);
x7:=g(x6);
x8:=g(x7);

#