OFFSET
0,1
COMMENTS
All terms, except a(0) and a(1), are integer multiples of 6.
LINKS
Harry J. Smith, Table of n, a(n) for n = 0..100
FORMULA
a(n) = Integral_{x=0..oo} (x^n*(exp(-sqrt(x)) * (-1+sqrt(x)+2/sqrt(x)) + x*Ei(-sqrt(x))) ), n=0, 1..., where Ei(y) is the exponential integral. Representation as the n-th moment of a positive function on a positive half-axis, in Maple notation. This representation is unique.
MATHEMATICA
Table[6*(2*n)!/(n+2), {n, 0, 20}] (* G. C. Greubel, May 03 2019 *)
PROG
(PARI) { s=6; for (n=0, 100, if (n, s*=2*n*(2*n - 1)); a=s/(n + 2); write("b064335.txt", n, " ", a) ) } \\ Harry J. Smith, Sep 12 2009
(PARI) a(n) = 6*(2*n)!/(n+2); \\ Michel Marcus, Jun 24 2018
(Magma) [6*Factorial(2*n)/(n+2): n in [0..20]]; // G. C. Greubel, May 03 2019
(Sage) [6*factorial(2*n)/(n+2) for n in (0..20)] # G. C. Greubel, May 03 2019
(GAP) List([0..20], n-> 6*Factorial(2*n)/(n+2)) # G. C. Greubel, May 03 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Karol A. Penson, Sep 13 2001
STATUS
approved