Source code in the dc programming language

1000 sn                 [ n = number of iterations ]sx
2574k                   [ precision to use during computation ]sx
1 d sf                  [ set e (kept on the stack) and f to 1 ]sx
[                       [ p = begin ]sx
  lf li 1 + d si / d sf   [ f = f / ++i ]sx
  +                       [ e = e + f ]sx
  ln li 

]sp                     [ end ]sx

lpx                     [ call p ]sx

[After ]n lnn           [ output header. uses n + 10P for newline ]sx
[ iterations, e =]n     [ instead of p, so nothing is left hanging ]sx
10P                     [ around on the stack.]sx

2570k                   [ now reset precision to match correct digits ]sx
1 /                     [ get result truncated to that precision ]sx
n 10P                   [ and print it out, again with n + 10P ]sx