|
Search: id:A106851
|
|
|
| A106851 |
|
Let M = {{0, 0, 0, 1}, {1, 4, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 4}}, v[1] = {0, 1, 1, 2}', v[n]=M.v[n-1]; then a(n) = v[n][[1]] |
|
+0
1
|
|
| 0, 2, 9, 37, 152, 626, 2585, 10701, 44400, 184610, 769065, 3209461, 13415048, 56153618, 235357241, 987609501, 4148575200, 17443003202, 73402179657, 309116995525, 1302649664888, 5492768393906, 23173154692697, 97810060234605
(list; graph; listen)
|
|
|
OFFSET
|
0,2
|
|
|
COMMENT
|
A 4 X 4 vector Markov chain with characteristic polynomial x^4-8*x^3+16*x^2-1..
Real-valued roots: {{x -> -0.236068}, {x -> 0.267949}, {x -> 3.73205}, {x -> 4.23607}}
|
|
MATHEMATICA
|
M = {{0, 0, 0, 1}, {1, 4, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 4}} v[1] = {0, 1, 1, 2}; v[n_] := v[n] = M.v[n - 1]; digits = 50; a = Table[v[n][[1]], {n, 1, digits}]
|
|
CROSSREFS
|
Sequence in context: A111601 A037553 A012493 this_sequence A129169 A162548 A150983
Adjacent sequences: A106848 A106849 A106850 this_sequence A106852 A106853 A106854
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Roger L. Bagula (rlbagulatftn(AT)yahoo.com), May 30 2005
|
|
EXTENSIONS
|
Edited by N. J. A. Sloane (njas(AT)research.att.com), Apr 09 2007
|
|
|
Search completed in 0.002 seconds
|
