|
%I #19 Sep 08 2022 08:45:45
%S 1,46,194,-167900,-8842004,884083016,125639477176,-4415829390416,
%T -1893481677885040,-19202364475675424,31870137298174352416,
%U 1835095760938501860416,-589384037754831073199936,-69436314367007836275831680,11532279106459848726285343616
%N Numerator of Hermite(n, 23/31).
%H Vincenzo Librandi, <a href="/A160334/b160334.txt">Table of n, a(n) for n = 0..200</a>
%F From _G. C. Greubel_, Jul 12 2018: (Start)
%F a(n) = 31^n * Hermite(n, 23/31).
%F E.g.f.: exp(46*x - 961*x^2).
%F a(n) = numerator(Sum_{k=0..floor(n/2)} (-1)^k*n!*(46/31)^(n-2*k)/(k!*(n-2*k)!)). (End)
%e Numerators of 1, 46/31, 194/961, -167900/29791, -8842004/923521, ...
%t Numerator[HermiteH[Range[0,20],23/31]] (* _Harvey P. Dale_, Aug 10 2014 *)
%t Table[31^n*HermiteH[n, 23/31], {n,0,30}] (* _G. C. Greubel_, Jul 12 2018 *)
%o (PARI) a(n)=numerator(polhermite(n,23/31)) \\ _Charles R Greathouse IV_, Jan 29 2016
%o (Magma) [Numerator((&+[(-1)^k*Factorial(n)*(46/31)^(n-2*k)/( Factorial(k) *Factorial(n-2*k)): k in [0..Floor(n/2)]])): n in [0..30]]; // _G. C. Greubel_, Jul 12 2018
%Y Cf. A009975 (denominators).
%K sign,frac
%O 0,2
%A _N. J. A. Sloane_, Nov 12 2009
| 