%I #19 Aug 17 2017 05:54:40
%S 0,1,21,501,14455,496770,19911486,913839031,47303189361,2727741976785,
%T 173455231572865,12060173714421756,910301022642409476,
%U 74134150415555474881,6479678618270868170265,605042444997867941987385,60110944381660549838273911
%N Bessel polynomial {y_n}'(3).
%D J. Riordan, Combinatorial Identities, Wiley, 1968, p. 77.
%H G. C. Greubel, <a href="/A065921/b065921.txt">Table of n, a(n) for n = 0..345</a>
%H <a href="/index/Be#Bessel">Index entries for sequences related to Bessel functions or polynomials</a>
%F Recurrence: (n-1)^2*a(n) = (2*n - 1)*(3*n^2 - 3*n + 1)*a(n-1) + n^2*a(n-2). - _Vaclav Kotesovec_, Jul 22 2015
%F a(n) ~ 2^(n+1/2) * 3^(n-1) * n^(n+1) / exp(n-1/3). - _Vaclav Kotesovec_, Jul 22 2015
%F From _G. C. Greubel_, Aug 14 2017: (Start)
%F a(n) = 2*n*(1/2)_{n}*6^(n - 1)* hypergeometric1f1(1 - n, -2*n, 2/3).
%F E.g.f.: ((1 - 6*x)^(3/2) + 3*x*(1 - 6*x)^(1/2) + 15*x - 1) * exp((1 - sqrt(1 - 6*x))/3)/(9*(1 - 6*x)^(3/2)). (End)
%F G.f.: (t/(1-t)^3)*hypergeometric2f0(2,3/2; - ; 6*t/(1-t)^2). - _G. C. Greubel_, Aug 16 2017
%t Table[Sum[(n+k+1)!*3^k/((n-k-1)!*k!*2^(k+1)), {k,0,n-1}], {n,0,20}] (* _Vaclav Kotesovec_, Jul 22 2015 *)
%t Join[{0}, Table[2*n*Pochhammer[1/2, n]*6^(n - 1)* Hypergeometric1F1[1 - n, -2*n, 2/3], {n,1,50}]] (* _G. C. Greubel_, Aug 14 2017 *)
%o (PARI) for(n=0,50, print1(sum(k=0,n-1, (n+k+1)!*3^k/((n-k-1)!*k! *2^(k+1))), ", ")) \\ _G. C. Greubel_, Aug 14 2017
%Y Cf. A001514, A065920, A065922.
%K nonn
%O 0,3
%A _N. J. A. Sloane_, Dec 08 2001