login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A158617 a(n) = Hermite(n, 16). 1
1, 32, 1022, 32576, 1036300, 32900992, 1042468744, 32964187904, 1040259450512, 32760875409920, 1029623343008224, 32292729468064768, 1010715629431891648, 31567874634586978304, 983893381941554122880, 30600687732361296539648 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

First negative term is a(140). - Georg Fischer, Feb 15 2019

LINKS

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

FORMULA

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

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

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

MATHEMATICA

With[{nmax = 50}, CoefficientList[Series[Exp[32*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(32*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(32*x - x^2))) \\ G. C. Greubel, Jul 13 2018

(PARI) for(n=0, 30, print1(polhermite(n, 16), ", ")) \\ G. C. Greubel, Jul 13 2018

CROSSREFS

Sequence in context: A299673 A299751 A300265 * A171337 A029548 A016745

Adjacent sequences:  A158614 A158615 A158616 * A158618 A158619 A158620

KEYWORD

sign

AUTHOR

N. J. A. Sloane, Nov 11 2009

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 6 15:28 EDT 2021. Contains 343586 sequences. (Running on oeis4.)