|
Search: id:A130636
|
|
|
| A130636 |
|
A vector matrix Markov doubly recursive sequence with characteristic polynomial: -1 - n x + x^4 and matrix: M(n)={{0, 0, n, 1}, {1, 0, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 0}}. |
|
+0 1
|
|
| 0, 1, 0, 0, 4, 1, 0, 28, 12, 1, 280, 160, 24, 3641, 2520, 520, 58280, 46481, 11880, 1107840, 987900, 295961, 24384360, 23829540, 8090964, 609904961, 643952400, 242285568, 17085429872, 19284524561, 7912519440
(list; graph; listen)
|
|
|
OFFSET
|
1,5
|
|
|
COMMENT
|
Determinants are all -1. The toral inverses appear to be Pisot: f[x_] = CharacteristicPolynomial[M[n], x] Table[NSolve[x^4*f[1/x] == 0, x], {n, 0, 10}] Table[Expand[x^4*f[1/x]], {n, 0, 10}]
|
|
FORMULA
|
M(0) = {{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {1, 0, 0, 1}}; M(n)={{0, 0, n, 1}, {1, 0, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 0}}; v(n)=M(n).v(n-1) a(n) = v(n)[[1]]
|
|
MATHEMATICA
|
M[0] = {{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {1, 0, 0, 1}}; M[n_] := {{0, 0, n, 1}, {1, 0, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 0}}; v[0] = {0, 0, 0, 1}; v[n_] := v[n] = M[n].v[n - 1]; a = Table[v[n][[1]], {n, 0, 30}]
|
|
CROSSREFS
|
Sequence in context: A007789 A081114 A069018 this_sequence A117414 A085639 A135302
Adjacent sequences: A130633 A130634 A130635 this_sequence A130637 A130638 A130639
|
|
KEYWORD
|
nonn,uned
|
|
AUTHOR
|
Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Jun 19 2007
|
|
|
Search completed in 0.002 seconds
|