OFFSET
0,4
LINKS
Index entries for linear recurrences with constant coefficients, signature (1, 19, 10, 1).
FORMULA
G.f.: (-18*x^3 + x)/(-x^4 - 10*x^3 - 19*x^2 - x + 1).
a(0)=0, a(1)=a(2)=1, a(3)=2, a(n)=a(n-1)+19a(n-2)+10a(n-3)+a(n-4) [From Harvey P. Dale, Apr 26 2011]
MATHEMATICA
M={{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {1, 10, 19, 1}} v[0] = {0, 1, 1, 2}; v[n_] := v[n] = M.v[n - 1] a = Table[v[n][[1]], {n, 0, 50}]
LinearRecurrence[{1, 19, 10, 1}, {0, 1, 1, 2}, 40] (* or *) CoefficientList[ Series[(-18x^3+x)/(-x^4-10x^3-19x^2-x+1), {x, 0, 40}], x] (* Harvey P. Dale, Apr 26 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Mar 31 2005
STATUS
approved