|
| |
|
|
A107402
|
|
a(n)= -a(n-1) +5*a(n-2) +5*a(n-3) -a(n-4) -a(n-5).
|
|
0
| |
|
|
0, 1, 1, 2, 3, 11, 12, 55, 55, 266, 261, 1277, 1248, 6121, 5977, 29330, 28635, 140531, 137196, 673327, 657343, 3226106, 3149517, 15457205, 15090240, 74059921, 72301681, 354842402, 346418163, 1700152091, 1659789132, 8145918055
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 0,4
|
|
|
FORMULA
| G.f.: -x*(5*x^3+2*x^2-2*x-1)/((x+1)*(x^4-5*x^2+1)). [From Maksym Voznyy (voznyy(AT)mail.ru), Aug 12 2009]
a(0)=0, a(2n) = (1/3)*(A055271(n)+2).
|
|
|
MATHEMATICA
| m = 5 M = {{0, 1, 0, 0, 0}, {0, 0, 1, 0, 0}, {0, 0, 0, 1, 0}, {0, 0, 0, 0, 1}, {-1, -1, m, m, -1}} Expand[Det[M - x*IdentityMatrix[5]]] NSolve[Det[M - x*IdentityMatrix[5]] == 0, x] v[1] = {0, 1, 1, 2, 3} digits = 50 v[n_] := v[n] = M.v[n - 1] a = Table[v[n][[1]], {n, 1, digits}]
|
|
|
CROSSREFS
| Sequence in context: A062936 A136972 A135115 * A020610 A030482 A019366
Adjacent sequences: A107399 A107400 A107401 * A107403 A107404 A107405
|
|
|
KEYWORD
| nonn,easy
|
|
|
AUTHOR
| Roger L. Bagula (rlbagulatftn(AT)yahoo.com), May 25 2005
|
|
|
EXTENSIONS
| Definition replaced by recurrence by the Associate Editors of the OEIS, Sep 28 2009
|
| |
|
|
