OFFSET
0,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Hermite Polynomial.
Index entries for linear recurrences with constant coefficients, signature (7,-21,35,-35,21,-7,1).
FORMULA
G.f.: (-120 + 1024*x - 4632*x^2 + 27968*x^3 + 23768*x^4 -2112* x^5 + 184*x^6)/(1-x)^7.
a(n) = 7*a(n-1)-21*a(n-2)+35*a(n-3)-35*a(n-4)+21*a(n-5)-7*a(n-6)+a(n-7).
MATHEMATICA
Table[64 n^6 - 480 n^4 + 720 n^2 - 120, {n, 0, 30}] (* or *) CoefficientList[Series[(-120 + 1024 x -4632 x^2 + 27968 x^3 + 23768 x^4 - 2112 x^5 + 184 x^6)/(1-x)^7, {x, 0, 30}], x]
LinearRecurrence[{7, -21, 35, -35, 21, -7, 1}, {-120, 184, -824, 14136, 150664, 717880, 2389704}, 30] (* Harvey P. Dale, Apr 08 2019 *)
PROG
(Magma) [64*n^6-480*n^4+720*n^2-120: n in [0..30]]; /* or */ I:=[-120, 184, -824, 14136, 150664, 717880, 2389704]; [n le 7 select I[n] else 7*Self(n-1)-21*Self(n-2)+35*Self(n-3)-35*Self(n-4)+21*Self(n-5)-7*Self(n-6)+Self(n-7): n in [1..30]]
(PARI) a(n)=polhermite(6, n) \\ Charles R Greathouse IV, Jan 29 2016
(Python)
from sympy import hermite
def A247851(n): return hermite(6, n) # Chai Wah Wu, Jan 06 2022
CROSSREFS
KEYWORD
sign,easy
AUTHOR
Vincenzo Librandi, Sep 25 2014
STATUS
approved