OFFSET
1,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..10000
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).
FORMULA
a(n) = n *Sum_{k=0..n-1} (-1)^k * (2*n-2*k-1)^3.
a(n+1) = (n+1) * ((2n+1)^3 - a(n)/n).
From G. C. Greubel, Jan 19 2019: (Start)
G.f.: x*(1 +47*x +47*x^2 +x^3)/(1-x)^5.
E.g.f.: x*(1 +25*x +24*x^2 +4*x^3)*exp(x). (End)
MATHEMATICA
Table[4*n^4-3*n^2, {n, 1, 40}] (* G. C. Greubel, Jan 19 2019 *)
PROG
(PARI) vector(40, n, 4*n^4-3*n^2) \\ G. C. Greubel, Jan 19 2019
(Magma) [4*n^4-3*n^2: n in [1..40]]; // G. C. Greubel, Jan 19 2019
(Sage) [4*n^4-3*n^2 for n in (1..40)] # G. C. Greubel, Jan 19 2019
(GAP) List([1..40], n -> 4*n^4-3*n^2); # G. C. Greubel, Jan 19 2019
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Kit Vongmahadlek (kit119(AT)yahoo.com), Jan 07 2003
EXTENSIONS
Edited by Don Reble, Nov 02 2005
STATUS
approved