%I #11 Jun 09 2019 18:18:52
%S 1,1,5,49,789,18741,612383,26218956,1419303189,94531262917,
%T 7582017897795,719690829785016,79691175192777855,10170046938232956048,
%U 1480481369981439216732,243659154929530351237884,44987315567879408248084629,9254611189980167520327621253
%N Doubly exponentiate the Bessel function I(0,2*sqrt(z)).
%H <a href="/index/Be#Bessel">Index entries for sequences related to Bessel functions or polynomials</a>
%F Hypergeometric generating function for a(n), in Maple notation: exp(exp(BesselI(0, 2*sqrt(z))-1)-1)=sum(a(n)*z^n/((n!)^2), n=0..infinity)
%t nmax = 20; CoefficientList[Series[E^(-1 + E^(-1 + BesselI[0, 2*Sqrt[x]])), {x, 0, nmax}], x] * Range[0, nmax]!^2 (* _Vaclav Kotesovec_, Jun 09 2019 *)
%Y Cf. A023998.
%K easy,nonn
%O 0,3
%A _Karol A. Penson_, Jun 28 2001
%E More terms from _Vaclav Kotesovec_, Jun 09 2019
|