|
|
A277002
|
|
Numerators of an asymptotic series for the Gamma function (odd power series).
|
|
4
|
|
|
-1, 7, -31, 127, -511, 1414477, -8191, 118518239, -5749691557, 91546277357, -23273283019, 1982765468311237, -22076500342261, 455371239541065869, -925118910976041358111, 16555640865486520478399, -1302480594081611886641, 904185845619475242495834469891
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Let y = x+1/2 then Gamma(x+1) ~ sqrt(2*Pi)*(y/E)^y*exp(Sum_{k>=1} r(k)/y^(2*k-1)) as x -> oo and r(k) = A277002(k)/A277003(k) (see example 7.1 in the Wang reference).
See also theorem 2 and formula (58) in Borwein and Corless. - Peter Luschny, Mar 31 2017
|
|
LINKS
|
|
|
FORMULA
|
a(n) = numerator(b(2*n-1)) with b(n) = Bernoulli(n+1, 1/2)/(n*(n+1)) for n>=1, b(0)=0.
|
|
EXAMPLE
|
The underlying rational sequence b(n) starts:
0, -1/24, 0, 7/2880, 0, -31/40320, 0, 127/215040, 0, -511/608256, ...
|
|
MAPLE
|
b := n -> `if`(n=0, 0, bernoulli(n+1, 1/2)/(n*(n+1))):
a := n -> numer(b(2*n-1)):
seq(a(n), n=1..18);
|
|
MATHEMATICA
|
b[n_] := BernoulliB[n+1, 1/2]/(n(n+1));
a[n_] := Numerator[b[2n-1]];
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign,frac
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|