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

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A158617 a(n) = Hermite(n, 16). 1

%I

%S 1,32,1022,32576,1036300,32900992,1042468744,32964187904,

%T 1040259450512,32760875409920,1029623343008224,32292729468064768,

%U 1010715629431891648,31567874634586978304,983893381941554122880,30600687732361296539648

%N a(n) = Hermite(n, 16).

%C First negative term is a(140). - _Georg Fischer_, Feb 15 2019

%H G. C. Greubel, <a href="/A158617/b158617.txt">Table of n, a(n) for n = 0..695</a>

%F From _G. C. Greubel_, Jul 13 2018: (Start)

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

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

%t With[{nmax = 50}, CoefficientList[Series[Exp[32*x - x^2], {x, 0, nmax}], x]*Range[0, nmax]!] (* _G. C. Greubel_, Jul 13 2018 *)

%o (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

%o (PARI) x='x+O('x^30); Vec(serlaplace(exp(32*x - x^2))) \\ _G. C. Greubel_, Jul 13 2018

%o (PARI) for(n=0,30, print1(polhermite(n, 16), ", ")) \\ _G. C. Greubel_, Jul 13 2018

%K sign

%O 0,2

%A _N. J. A. Sloane_, Nov 11 2009

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.

Last modified June 15 15:09 EDT 2021. Contains 345049 sequences. (Running on oeis4.)