A158752 a(n) = Hermite(n, 23).

1, 46, 2114, 97060, 4452076, 204019016, 9340353976, 427208054704, 19520805560720, 891121726917856, 40640224938128416, 1851627912615550016, 84280799031676475584, 3832477685554344676480, 174102672760676266752896

Offset

0,2

Comments

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

Links

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

Formula

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

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

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

Mathematica

Table[HermiteH[n, 23], {n, 0, 50}] (* or *) With[{nmax = 50}, CoefficientList[ Series[ Exp[ 46*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(46*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(46*x - x^2))) \\ G. C. Greubel, Jul 13 2018
(PARI) for(n=0, 30, print1(polhermite(n, 23), ", ")) \\ G. C. Greubel, Jul 13 2018

Crossrefs

Sequence in context: A170727 A170765 A218748 * A223738 A223721 A223827

Adjacent sequences: A158749 A158750 A158751 * A158753 A158754 A158755

Keyword

sign

Author

N. J. A. Sloane, Nov 11 2009

Status

approved