login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A098055 Let M={{0,1,0},{0,0,1},{1,1,0}}, M0=MatrixPower[(M-IdentityMatrix[3]),2], Det[M0], a[n_]:=M0.a[n-1]; a[0]:={{0,1,1},{1,1,1},{1,1,2}}. 0
0, 1, 1, 1, 1, 1, 1, 1, 2, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 3, 4, 3, 4, 3, 4, 3, 1, 10, 14, 11, 14, 11, 4, 11, 4, 3, 49, 40, 16, 40, 16, 9, 16, 9, 24, 145, 63, 26, 63, 26, 82, 26, 82, 89, 245, 71, 279, 71, 279, 316, 279, 316, 208, 176, 945, 1119, 945, 1119, 769, 1119, 769, 174, 3185 (list; graph; refs; listen; history; internal format)
OFFSET

0,9

COMMENTS

3 X 3 matrix from the minimal Pisot generator matrix by: (M-I)^2.

MATHEMATICA

(*square matrix 3 X 3 Markov sequence*) Clear[x, M, A] digits=21 M={{0, 1, 0}, {0, 0, 1}, {1, 1, 0}} M0=MatrixPower[(M-IdentityMatrix[3]), 2] Det[M0] A[n_]:=M0.A[n-1]; A[0]:={{0, 1, 1}, {1, 1, 1}, {1, 1, 2}}; (* flattened sequence of 3 X 3 matrices made with an alternating recurrence*) b=Flatten[Table[Abs[A[n]], {n, 0, digits}]] ListPlot[b, PlotJoined->True]

CROSSREFS

Sequence in context: A079748 A073368 A037889 * A092111 A050317 A141095

Adjacent sequences:  A098052 A098053 A098054 * A098056 A098057 A098058

KEYWORD

nonn

AUTHOR

Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Sep 11 2004

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified February 15 20:26 EST 2012. Contains 205852 sequences.