OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..300
Eric Weisstein's World of Mathematics, Laguerre Polynomial
Wikipedia, Laguerre polynomials
FORMULA
E.g.f.: 2*K(-16*x)/Pi, where K() is the complete elliptic integral of the first kind.
a(n) ~ (-1)^n * 16^n * (n-1)! / Pi. - Vaclav Kotesovec, Nov 21 2017
From Peter Luschny, Nov 21 2017: (Start)
a(n) = (-16)^n*Gamma(n+1/2)^2/(Pi*Gamma(n+1)).
a(n) = (-16)^n*binomial(n-1/2,-1/2)*Gamma(n+1/2)/sqrt(Pi).
a(n) ~ (-exp(-1)*n*16)^n/sqrt(n*Pi/2). (End)
a(n) = (-1)^n*binomial(2*n,n)^2*n!. - Alois P. Heinz, Oct 02 2021
MAPLE
a := n -> (-16)^n*GAMMA(n+1/2)^2/(Pi*GAMMA(n+1)):
seq(a(n), n=0..15); # Peter Luschny, Nov 21 2017
MATHEMATICA
Table[(2 n)! SeriesCoefficient[(-x/(1 - x))^n /((1 - x) n!), {x, 0, 2 n}], {n, 0, 15}]
nmax = 15; CoefficientList[Series[2 EllipticK[-16 x]/Pi, {x, 0, nmax}], x] Range[0, nmax]!
Table[(-16)^n*Gamma[n + 1/2]^2/(Pi*Gamma[n + 1]), {n, 0, 50}] (* G. C. Greubel, Feb 06 2018 *)
PROG
(PARI) for(n=0, 30, print1(round((-16)^n*gamma(n+1/2)^2/(Pi*gamma(n+1))), ", ")) \\ G. C. Greubel, Feb 06 2018
(Magma) R:= RealField(); [Round((-16)^n*Gamma(n+1/2)^2/(Pi(R)*Gamma(n+1) )): n in [0..30]]; // G. C. Greubel, Feb 06 2018
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Nov 21 2017
STATUS
approved