login
Denominators of numbers occurring in continued fraction connected with expansion of gamma function.
(Formerly M4831)
4

%I M4831 #32 May 08 2018 15:11:54

%S 12,30,210,371,22737,19733142,48264275462,9769214287853155785,

%T 113084128923675014537885725485,

%U 5271244267917980801966553649147604697542,24274291553105128438297398108902195365373879212227726,13346384670164266280033479022693768890138348905413621178450736182873

%N Denominators of numbers occurring in continued fraction connected with expansion of gamma function.

%D M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math.Series 55, Tenth Printing, 1972, p. 258.

%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

%H Vincenzo Librandi, <a href="/A005147/b005147.txt">Table of n, a(n) for n = 0..30</a>

%H 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 [alternative scanned copy].

%H 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. 258.

%H B. W. Char, <a href="http://www.jstor.org/stable/2006103">On Stieltjes' continued fraction for the gamma function</a>, Math. Comp., 34 (1980), 547-551.

%H Peter Luschny, <a href="/A005146/a005146.txt">Maple program for A005146/A005147</a>

%H Peter Luschny, <a href="http://www.luschny.de/math/factorial/approx/continuedfraction.html">Continued fraction</a>

%H 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).

%t len = 12; s[p_] := (-1)^p * BernoulliB[2p+2]/(2p+1)/(2p+2); Do[m[n, 1] = 0, {n, 0, len}]; Do[m[n, 2] = s[n+1]/s[n], {n, 0, len-1}]; Do[m[n, k] =

%t If[OddQ[k], m[n+1, k-2]+m[n+1, k-1]-m[n, k-1],

%t m[n+1, k-2]*m[n+1, k-1]/m[n, k-1]], {k, 3, len}, {n, 0,

%t len-k+1}]; Do[m[n, 1] = s[n], {n, 0, len}];

%t Table[m[0, k], {k, 1, len}] // Denominator

%t (* _Jean-François Alcover_, May 24 2011, after _Peter Luschny_ *)

%Y Cf. A005146.

%K nonn,frac,nice

%O 0,1

%A _Simon Plouffe_ and _N. J. A. Sloane_

%E More terms from _Rainer Rosenthal_, Jan 11 2007