% Script fastroot.m

% quickly covergent iteration to find root of
% x - exp(-x) = 0 using Newton's method

x0 = 0.5;
x = x0;
for n = 1:6
   e = exp(-x);
   x = x + (e - x)/(e + 1);
   fprintf(1, '%5d %25.15e\n', n, x);
end