OFFSET
1,2
LINKS
Index entries for linear recurrences with constant coefficients, signature (4,0,-8,9).
FORMULA
From Colin Barker, Oct 19 2012: (Start)
a(n) = 4*a(n-1) -8*a(n-3) +9*a(n-4).
G.f.: -x*(x-1)*(3*x^2-1)/(9*x^4-8*x^3+4*x-1). (End)
MAPLE
seq(coeff(series(-x*(x-1)*(3*x^2-1)/(9*x^4-8*x^3+4*x-1), x, n+1), x, n), n = 1 .. 25); # Muniru A Asiru, Oct 13 2018
MATHEMATICA
LinearRecurrence[{4, 0, -8, 9}, {1, 3, 9, 31}, 30] (* Harvey P. Dale, Jul 26 2018 *)
PROG
(PARI) x='x+O('x^30); Vec(-x*(x-1)*(3*x^2-1)/(9*x^4-8*x^3+4*x-1)) \\ G. C. Greubel, Oct 12 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(-x*(x-1)*(3*x^2-1)/(9*x^4-8*x^3+4*x-1))); // G. C. Greubel, Oct 12 2018
CROSSREFS
KEYWORD
nonn,less
AUTHOR
Roger L. Bagula, Oct 05 2006
EXTENSIONS
Sequence edited by Joerg Arndt and Colin Barker, Oct 19 2012
STATUS
approved