Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A122074
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A122074 a(0)=1, a(1)=6, a(n)=7a(n-1)-2a(n-2). +0
1
1, 6, 40, 268, 1796, 12036, 80660, 540548, 3622516, 24276516, 162690580, 1090281028, 7306586036, 48965540196, 328145609300, 2199088184708, 14737326074356, 98763106151076, 661867090908820, 4435543424059588 (list; graph; listen)
OFFSET

0,2

COMMENT

First row sum of the matrix M^n, where M is the 3 X 3 matrix {{2,2,2},{2,3,2},{2,2,3}}.

REFERENCES

Peter Steinbach, "Golden Fields: A Case for the Heptagon", Mathematics Magazine, Vol. 70, No. 1, Feb. 1997.

FORMULA

a(n)=8a(n-1)-9a(n-2)+2a(n-3); a(0)=1,a(1)=6,a(2)=40 (follows from the minimal polynomial x^3-8x^2+9x-2 of M).

EXAMPLE

a(2)=40 because M^2={{12,14,14},{14,17,16},{14,16,17}} and 12+14+14=40.

MAPLE

with(linalg): M[1]:=matrix(3, 3, [2, 2, 2, 2, 3, 2, 2, 2, 3]): for n from 2 to 20 do M[n]:=multiply(M[n-1], M[1]) od: 1, seq(M[n][1, 1]+M[n][1, 2]+M[n][1, 3], n=1..20);

MATHEMATICA

M = {{2, 2, 2}, {2, 3, 2}, {2, 2, 3}}; v[1] = {1, 1, 1}; v[n_] := v[n] = M.v[n - 1]; a1 = Table[v[n][[1]], {n, 1, 50}]

CROSSREFS

Sequence in context: A065113 A052518 A135032 this_sequence A123357 A081016 A083426

Adjacent sequences: A122071 A122072 A122073 this_sequence A122075 A122076 A122077

KEYWORD

nonn

AUTHOR

Gary Adamson (qntmpkt(AT)yahoo.com), Oct 16 2006

EXTENSIONS

Edited by N. J. A. Sloane (njas(AT)research.att.com), Oct 29 2006 and Dec 04 2006

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified November 25 20:09 EST 2009. Contains 167514 sequences.


AT&T Labs Research