OFFSET
0,1
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (1,6,-6,-1,1).
FORMULA
a(n) = a(n-1) + 6*a(n-2) - 6*a(n-3) - a(n-4) + a(n-5).
G.f.: 216*(3 - 28*z + 78*z^2 + 4*z^3 - 13*z^4)/((1 - z)*(1 + 2*z - z^2) *(1 - 2*z - z^2)).
E.g.f.: 216*(cosh(x)*(14*cosh(sqrt(2)*x) - 4*sqrt(2)*sinh(sqrt(2)*x) - 11) - sinh(x)*(6*cosh(sqrt(2)*x) - 10*sqrt(2)*sinh(sqrt(2)*x) + 11)). - Stefano Spezia, Oct 03 2022
MATHEMATICA
uu = {648, -5400, 15336, -20088, 100872}; a1 = aa[[1]]; a2 = aa[[2]]; a3 = aa[[3]]; a4 = aa[[4]]; a5 = aa[[5]]; Do[an = a5 + 6 a4 - 6 a3 - a2 + a1; a1 = a2; a2 = a3; a3 = a4; a4 = a5; a5 = an; AppendTo[uu, an], {nn, 1, 20}]; uu
PROG
(PARI) my(x='x+O('x^30)); Vec(216*(3-28*x+78*x^2+4*x^3-13*x^4)/((1-x)*(1+2*x-x^2)*(1-2*x-x^2))) \\ G. C. Greubel, Aug 18 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(216*(3-28*x+78*x^2+4*x^3-13*x^4)/((1-x)*(1+2*x-x^2)*(1-2*x-x^2)))); // G. C. Greubel, Aug 18 2018
CROSSREFS
KEYWORD
sign,easy,changed
AUTHOR
Artur Jasinski, Nov 25 2011
STATUS
approved