cp's OEIS Frontend

A001164 Stirling's formula: denominators of asymptotic series for Gamma function.

%I A001164 M4878 N2091 #77 Mar 31 2025 04:47:57
%S A001164 1,12,288,51840,2488320,209018880,75246796800,902961561600,
%T A001164 86684309913600,514904800886784000,86504006548979712000,
%U A001164 13494625021640835072000,9716130015581401251840000,116593560186976815022080000,2798245444487443560529920000,299692087104605205332754432000000,57540880724084199423888850944000000
%N A001164 Stirling's formula: denominators of asymptotic series for Gamma function.
%D A001164 L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 267, #23.
%D A001164 N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
%D A001164 N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%D A001164 Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapters 2 and 43, equations 2:6:1 and 43:6:6 at pages 21, 415.
%H A001164 Seiichi Manyama, <a href="/A001164/b001164.txt">Table of n, a(n) for n = 0..295</a> (terms 0..100 from T. D. Noe)
%H A001164 M. Abramowitz and I. A. Stegun, eds., <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of Mathematical Functions</a>, National Bureau of Standards Applied Math.Series 55, Tenth Printing, 1972, p. 257, Eq. 6.1.37.
%H A001164 S. Brassesco and M. A. Méndez, <a href="http://arxiv.org/abs/1002.3894">The asymptotic expansion for the factorial and Lagrange inversion formula</a>, arXiv:1002.3894 [math.CA], 2010.
%H A001164 V. De Angelis, <a href="http://www.jstor.org/stable/40391303">Stirling's series revisited</a>, Amer. Math. Monthly, 116 (2009), 839-843.
%H A001164 Peter Luschny, <a href="http://oeis.org/wiki/User:Peter_Luschny/FactorialFunction">Approximations to the factorial function</a>.
%H A001164 G. Marsaglia and J. C. W. Marsaglia, <a href="http://www.jstor.org/stable/2324749">A new derivation of Stirling's approximation to n!</a>, Amer. Math. Monthly, 97 (1990), 827-829. MR1080390 (92b:41049)
%H A001164 T. Müller, <a href="http://dx.doi.org/10.1007/BF01195003">Finite group actions and asymptotic expansion of e^P(z)</a>, Combinatorica, 17 (4) (1997), 523-554.
%H A001164 Richard M. Slevinsky, <a href="https://arxiv.org/abs/1602.02618">On the use of Hahn's asymptotic formula and stabilized recurrence for a fast, simple, and stable Chebyshev-Jacobi transform</a>, arXiv preprint arXiv:1602.02618 [math.NA], 2016.
%H A001164 N. M. Temme, <a href="http://dx.doi.org/10.1137/0510071">The asymptotic expansion of the incomplete gamma function</a>, SIAM J. Math. Anal., 10 (1979), 757-766. [From _N. J. A. Sloane_, Feb 20 2012]
%H A001164 Nico Temme, <a href="http://arxiv.org/abs/math/9603218">Uniform Asymptotics for the incomplete gamma functions starting from negative values of the parameters</a>, arXiv:math/9603218 [math.CA], 1996.
%H A001164 W. Wang, <a href="http://dx.doi.org/10.1016/j.jnt.2015.12.016">Unified approaches to the approximations of the gamma function</a>, J. Number Theory (2016).
%H A001164 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/StirlingsSeries.html">Stirlings Series</a>
%H A001164 J. W. Wrench, Jr., <a href="http://dx.doi.org/10.1090/S0025-5718-1968-0237078-4">Concerning two series for the gamma function</a>, Math. Comp., 22 (1968), 617-626.
%F A001164 The coefficients c_k have g.f. 1 + Sum_{k >= 1} c_k/z^k = exp( Sum_{k >= 1} B_{2k}/(2k(2k-1)z^(2k-1)) ).
%F A001164 Numerators/denominators: A001163(n)/A001164(n) = (6*n+1)!!/(4^n*(2*n)!) * Sum_{i=0..2*n} Sum_{j=0..i} Sum_{k=0..j} (-1)^k*2^i*k^(2*n+i+j)*C(2*n,i)* C(i,j)*C(j,k)/((2*n+2*i+1)*(2*n+i+j)!), assuming 0^0 = 1 (when n = 0), n!! = A006882(n), C(n,k) = A007318(n,k) are binomial coefficients. - _Vladimir Reshetnikov_, Nov 04 2015
%F A001164 a(n) = denominator(h(2*n)*doublefactorial(2*n-1)) where h(k) = (h(k-1)/k - Sum_{j=1..k-1} h(k-j)*h(j)/(j+1))/(1+1/(k+1)) and h(0)=1. - _Peter Luschny_, Nov 05 2015
%F A001164 Let B_n be the Bernoulli number, and define the sequence {c_n} by the recurrence
%F A001164 c_0 = 1, c_n = (1/n) * Sum_{k=0..n-1} B_{n-k+1}*c_k/(n-k+1) for n > 0. Then a(n) is the denominator of c_n. - _Seiichi Manyama_, Sep 01 2018
%e A001164 Gamma(z) ~ sqrt(2 Pi) z^(z-1/2) e^(-z) (1 + 1/(12 z) + 1/(288 z^2) - 139/(51840 z^3) - 571/(2488320 z^4) + ... ), z -> infinity in |arg z| < Pi.
%p A001164 h := proc(k) option remember; local j; `if`(k=0, 1,
%p A001164 (h(k-1)/k-add((h(k-j)*h(j))/(j+1),j=1..k-1))/(1+1/(k+1))) end:
%p A001164 coeffStirling := n -> h(2*n)*doublefactorial(2*n-1):
%p A001164 seq(denom(coeffStirling(n)), n=0..16); # _Peter Luschny_, Nov 05 2015
%t A001164 Denominator[ Reverse[ Drop[ CoefficientList[ Simplify[ PowerExpand[ Normal[ Series[n!, {n, Infinity, 17}]]Exp[n]/(Sqrt[2Pi n]*n^(n - 17))]], n], 1]]]
%t A001164 h[k_] := h[k] = If[k==0, 1, (h[k-1]/k-Sum[h[k-j]*h[j]/(j+1), {j, 1, k-1}]) / (1+1/(k+1))]; StirlingAsympt[n_] := h[2n]*2^n*Pochhammer[1/2, n]; a[n_] := StirlingAsympt[n] // Denominator; Table[a[n], {n, 0, 30}] (* _Jean-François Alcover_, Oct 12 2015, after _Peter Luschny_ *)
%o A001164 (PARI) a(n)=local(A,m); if(n<1,n==0,A=vector(m=2*n+1,k,1); for(k=2,m,A[k]=(A[k-1]-sum(i=2,k-1,i*A[i]*A[k+1-i]))/(k+1)); denominator(A[m]*m!/2^n/n!)) /* _Michael Somos_, Jun 09 2004 */
%Y A001164 Product_{z=1..n} z^(z^m): A001163/A001164 (m=0), A143475/A143476 (m=1), A317747/A317796(m=2).
%K A001164 nonn,frac,nice
%O A001164 0,2
%A A001164 _N. J. A. Sloane_
%E A001164 More terms from _Vladeta Jovovic_, Nov 14 2001