|
|
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
|
|
|
LINKS
|
|
|
FORMULA
|
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
|
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|