OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..500
Eric Weisstein's World of Mathematics, Hermite Polynomial
Wikipedia, Hermite polynomials
FORMULA
E.g.f.: exp(exp(2*x)-x^2-1).
a(n) ~ 2^n * n^n * exp(n/LambertW(n) - LambertW(n)^2 / 4 - n - 1) / (sqrt(1 + LambertW(n)) * LambertW(n)^n). - Vaclav Kotesovec, Jun 29 2022
MAPLE
g := exp(exp(2*x)-x^2-1): ser := series(g, x, 32):
seq(n!*coeff(ser, x, n), n=0..21); # Peter Luschny, Oct 19 2016
MATHEMATICA
Table[Sum[HermiteH[n, k]/k!, {k, 0, Infinity}]/E, {n, 0, 20}] (* or *)
Table[SeriesCoefficient[Exp[Exp[2 x] - x^2 - 1], {x, 0, n}] n!, {n, 0, 20}]
With[{nmax = 50}, CoefficientList[Series[Exp[Exp[2*x] - x^2 - 1], {x, 0, nmax}], x]*Range[0, nmax]!] (* G. C. Greubel, Jul 12 2018 *)
Table[Sum[Binomial[n, k] * 2^k * BellB[k] * HermiteH[n-k, 0], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Jul 13 2018 *)
PROG
(PARI) x='x+O('x^50); Vec(serlaplace(exp(exp(2*x)-x^2-1))) \\ G. C. Greubel, Jul 12 2018
(Magma) m:=50; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(Exp(Exp(2*x) -x^2 -1))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Jul 12 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Reshetnikov, Oct 11 2016
STATUS
approved