OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..449 (terms 0..200 from Alois P. Heinz)
Peter Luschny, Approximations to the Factorial Function.
Eric Weisstein's World of Mathematics, Stirling's Approximation.
MAPLE
Digits:= 2000;
a:= n-> round(sqrt((2*n+1/3)*Pi)*n^n/exp(n)):
seq(a(n), n=0..30); # Alois P. Heinz, Feb 04 2013
MATHEMATICA
Join[{1}, Table[Round[Sqrt[(2*n + 1/3)*Pi]*n^n/Exp[n]], {n, 1, 50}]] (* G. C. Greubel, Nov 28 2017 *)
PROG
(PARI) a(n) = round(sqrt((2*n+1/3)*Pi)*n^n/exp(n)); \\ Bill McEachen, Aug 16 2014
(Magma) C<i> := ComplexField(); [Round(Sqrt((2*n + 1/3)*Pi(C))*n^n/Exp(n)): n in [0..30]]; // G. C. Greubel, Nov 28 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, Jan 10 2004
STATUS
approved