OFFSET
0,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..899
FORMULA
a(n) = Sum_{k=0..n} 3^k*binomial(n,k)*binomial(2*k,k).
a(n) = [x^n] (1+7*x+9*x^2)^n.
n * a(n) = 7 * (2*n-1) * a(n-1) - 13 * (n-1) * a(n-2) for n > 1.
E.g.f.: exp(7*x) * BesselI(0,6*x). - Ilya Gutkovskiy, Feb 01 2021
a(n) ~ 13^(n + 1/2) / (2 * sqrt(3*Pi*n)). - Vaclav Kotesovec, Nov 13 2021
MATHEMATICA
a[n_] := Sum[3^k * Binomial[n, k] * Binomial[2*k, k], {k, 0, n}]; Array[a, 18, 0] (* Amiram Eldar, Feb 01 2021 *)
nxt[{n_, a_, b_}]:={n+1, b, (7*b(2n+1)-13*n*a)/(n+1)}; Join[{1}, NestList[nxt, {2, 7, 67}, 20] [[All, 2]]] (* Harvey P. Dale, Apr 27 2022 *)
PROG
(PARI) my(N=20, x='x+O('x^N)); Vec(1/sqrt((1-x)*(1-13*x)))
(PARI) {a(n) = sum(k=0, n, 3^k*binomial(n, k)*binomial(2*k, k))}
(PARI) {a(n) = polcoef((1+7*x+9*x^2)^n, n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Feb 01 2021
STATUS
approved