dc/.dc

   1 5k
   2 10o
   3
   4 # [ 1-d2+/
   5 #   d2*Sk              # Initialize multiplier
   6 #   d*Sy               # Initialize multiplier factor
   7 #   0Ss                # Initialize accumulator
   8 #   10K^Sp             # Initialize 10^k power
   9 #   [ d1r/lk*ls+lsrdss   # Update accumulator
  10 #     -lp*d*1 [s_q]s_>_  # Check precision
  11 #     lkly*sk 2+         # Update multiplier and counter
  12 #     lfx
  13 #   ]Sf
  14 #   1lfxLs
  15 # LkLyLpLfs_s_s_s_] sL
  16 # natural logarithm: 1000lLx 10lLx /p
  17 [1-d2+/d2*Skd*Sy0Ss10K^Sp[d1r/lk*ls+lsrdss-lp*d*1[s_q]s_>_lkly*sk2+lfx]Sf1lfxLsLkLyLpLfs_s_s_s_]sL
  18
  19 # sum of everything in stack: lSxp
  20 [0d[+2z>a]salax]sS