OFFSET
0,4
COMMENTS
Characteristic polynomial of the matrix M is x^6-(x+1)^5.
LINKS
Index entries for linear recurrences with constant coefficients, signature (1, 5, 10, 10, 5, 1).
FORMULA
Recurrence relation: a(n)=a(n-1)+5a(n-2)+10a(n-3)+10a(n-4)+5a(n-5)+a(n-6) for n>=6; a(0)=0,a(1)=a(2)=1,a(3)=2,a(4)=3,a(5)=5.
O.g.f.: x*(2*x+1)*(14*x^3+2*x-1)/(-1+x+5*x^2+10*x^3+10*x^4+5*x^5+x^6) . - R. J. Mathar, Dec 05 2007
MAPLE
a[0]:=0:a[1]:=1:a[2]:=1:a[3]:=2:a[4]:=3:a[5]:=5: for n from 6 to 25 do a[n]:=a[n-1]+5*a[n-2]+10*a[n-3]+10*a[n-4]+5*a[n-5]+a[n-6] od: seq(a[n], n=0..25);
MATHEMATICA
M = {{0, 1, 0, 0, 0}, {0, 0, 1, 0, 0}, {0, 0, 0, 1, 0}, {0, 0, 0, 0, 1}, {1, 4, 6, 4, 1}}; w[0] = {0, 1, 1, 2, 3}; w[n_] := w[n] = M.w[n - 1] a = Flatten[Table[w[n][[1]], {n, 0, 25}]]
LinearRecurrence[{1, 5, 10, 10, 5, 1}, {0, 1, 1, 2, 3, 5}, 30] (* Harvey P. Dale, Jun 15 2014 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Feb 18 2006
EXTENSIONS
Edited by N. J. A. Sloane, May 13 2006
STATUS
approved