%I #15 Sep 08 2022 08:45:43
%S 1,16,254,4000,62476,967616,14857096,226102144,3409634960,50936525056,
%T 753610971616,11039045044736,160045279340224,2295787388369920,
%U 32571420951072896,456860688342808576,6332628384952750336,86702512132274130944,1171930829027992583680
%N a(n) = Hermite(n,8).
%C The first negative term is a(40) = -15216726248438640506371403737952889391284224. - _Robert Israel_, Apr 30 2018
%H Robert Israel, <a href="/A158531/b158531.txt">Table of n, a(n) for n = 0..722</a>
%F E.g.f.: exp(-x^2+16*x). a(n+2)=16*a(n+1)-(2*n+2)*a(n). - _Robert Israel_, Apr 30 2018
%p seq(orthopoly[H](n,8),n=0..100); # _Robert Israel_, Apr 30 2018
%t Table[HermiteH[n, 8], {n, 0, 40}] (* _Vincenzo Librandi_, May 01 2018 *)
%o (PARI) a(n) = polhermite(n, 8); \\ _Michel Marcus_, May 01 2018
%o (Magma) I:=[1,16]; [n le 2 select I[n] else 16*Self(n-1)-(2*n-4)*Self(n-2): n in [1..25]]; // _Vincenzo Librandi_, May 01 2018
%K sign
%O 0,2
%A _N. J. A. Sloane_, Nov 11 2009
|