OFFSET
0,3
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..200
FORMULA
a(n) = 1+sum((-1)^(n-k)*C(n,k)*C(2k-1,k-1)*(2k)!/k!, k=0..n).
Recurrence: n>=3, a(n) = (2*(-1)^n + (32 - 48*n + 16*n^2)*a(n-3) + (33 - 65*n + 32*n^2)*a(n-2) + (5 - 18*n + 16*n^2)*a(n-1))/n
E.g.f.: exp(-x) (1/2 + 1/pi K(16x) ), where K(z) is the elliptic integral of the first kind (defined as in Mathematica).
a(n) ~ 16^n*n^(n-1/2)/(sqrt(2*Pi)*exp(n+1/16)). - Vaclav Kotesovec, Aug 10 2013
MAPLE
seq((-1)^n+add((-1)^(n-k)*binomial(n, k)*binomial(2*k-1, k-1)*(2*k)!/k!, k=1..n), n=0..20);
MATHEMATICA
Table[(-1)^n + Sum[(-1)^(n-k)Binomial[n, k]Binomial[2k-1, k-1](2k)!/k!, {k, 1, n}], {n, 0, 20}]
PROG
(Maxima) makelist((-1)^n+sum((-1)^(n-k)*binomial(n, k)*binomial(2*k-1, k-1) *(2*k)!/k!, k, 1, n), n, 0, 12);
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Emanuele Munarini, Mar 11 2011
STATUS
approved
