OFFSET
0,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..100
FORMULA
Sum_{n>=0} a(n)*x^n/n!^2 = BesselI(0,2*sqrt(x/(1-x))).
Special values of hypergeometric function of type 1F2. In Maple notation: a(n)=((n!)^2)*hypergeom([1-n],[2,2],-1), n=0,1... . This sequence arises in exponentiating the operator D=d(x^2)(d^2), where d=d/dx. - Karol A. Penson, Nov 22 2008
Recurrence: a(n) = (3*n^2-5*n+3)*a(n-1) + (n-3)*(n-1)^2*(n-2)^3*a(n-3) - (n-1)^2*(3*n-4)*(n-2)*a(n-2). - Vaclav Kotesovec, Jun 03 2013
a(n) ~ n^(2*n)*exp(3*n^(1/3)-2*n)/sqrt(3). - Vaclav Kotesovec, Jun 03 2013
MATHEMATICA
CoefficientList[Series[BesselI[0, 2*Sqrt[x/(1-x)]], {x, 0, 20}], x]* Range[0, 20]!^2 (* Vaclav Kotesovec, Jun 03 2013 *)
PROG
(PARI) a(n)={if(n<1, n==0, sum(k=0, n, (n!/k!)^2*binomial(n-1, k-1)))} \\ Andrew Howroyd, Jan 08 2020
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Vladeta Jovovic, Jul 25 2006
EXTENSIONS
Terms a(15) and beyond from Andrew Howroyd, Jan 08 2020
STATUS
approved