restart: Digits:=24: f := n->n/(n!)^(1/n); for k from 0 to 15 do approx[k] := evalf(f(2^k)); od; for k from 0 to 15 do err[k] := evalf(exp(1)-approx[k]); od;