% Example 1.1: nummathI: Approximation of e
for i=1:15
  n = 10^i; r = (1+1/n)^n;
  fprintf('10^%2d %20.15f %20.15f\n',i,r,r-exp(1));
end