OFFSET
0,1
COMMENTS
MATHEMATICA
LogBaseBStirling[b_, n_] := Block[{}, N[ Log[b, 2*Pi*n]/2 + n*Log[b, n/E] + Log[b, 1 + 1/(12n) + 1/(288n^2) - 139/(51840n^3) - 571/(2488320n^4) + 163879/(209018880n^5)], 64]]; f[0] = 2; f[n_] := f[n] = Block[{k = 10*g[n - 1]}, While[ LogBaseBStirling[10^n, k] <= k, k++ ]; k]; Table[ f[n], {n, 1, 18}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Jul 13 2003
STATUS
approved