A158783 a(n) = Hermite(n, 24).

1, 48, 2302, 110304, 5280780, 252595008, 12071752584, 576412983936, 27498818692752, 1310720689509120, 62419614359968224, 2969927075488292352, 141183268107518731968, 6705518619349180118016, 318194128757965158633600

Offset

0,2

Comments

The first negative term is a(304). - Georg Fischer, Feb 16 2019

Links

G. C. Greubel, Table of n, a(n) for n = 0..650

Formula

From G. C. Greubel, Jul 13 2018: (Start)

E.g.f.: exp(48*x - x^2).

a(n) = 48*a(n-1) - 2*(n-1)*a(n-2). (End)

Mathematica

HermiteH[Range[0, 20], 24] (* Harvey P. Dale, Aug 10 2011 *)
With[{nmax = 50}, CoefficientList[Series[Exp[48*x - x^2], {x, 0, nmax}], x]*Range[0, nmax]!] (* G. C. Greubel, Jul 13 2018 *)

Prog

(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(Exp(48*x - x^2))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Jul 13 2018
(PARI) x='x+O('x^30); Vec(serlaplace(exp(48*x - x^2))) \\ G. C. Greubel, Jul 13 2018
(PARI) for(n=0, 30, print1(polhermite(n, 24), ", ")) \\ G. C. Greubel, Jul 13 2018

Crossrefs

Sequence in context: A170767 A218750 A263504 * A227139 A009992 A042105

Adjacent sequences: A158780 A158781 A158782 * A158784 A158785 A158786

Keyword

sign

Author

N. J. A. Sloane, Nov 11 2009

Status

approved