The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A094943 A sequence generated from a semi-magic square. 2
1, 13, 72, 429, 2601, 15534, 93339, 559845, 3359232, 20155473, 120932109, 725594598, 4353563943, 26121388761, 156728328192, 940369966965, 5642219821473, 33853318876350, 203119913356515, 1218719480001309 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
The 3 rows: 1 3 2, 2 1 3 and 3 2 1 form a semi-magic square; with rows, columns and a diagonal having a sum of 6. a(n)/a(n-1) tends to 6, an eigenvalue of the matrix. E.g.: a(7)/a(6) = 93339/15534 = 6.0086... A094944 uses the same format and operations but has different terms.
LINKS
FORMULA
G.f.: (1+10*x+18*x^2)/(1-3*x-18*x^3-15*x^2); a(n+3)=3*a(n+2)+15*a(n+1)+18*a(n), a(0) = 1, a(1) = 13, a(2) = 72. - Alec Mihailovs, Aug 28 2005
Let [1 3 2 / 2 1 3 / 3 2 1] = the 3 X 3 matrix M. Take M^n * [1 0 0] = [p q r]; then a(n) = p.
EXAMPLE
a(4) = 429 since M^4 * [1 0 0] = [429 q r]
MATHEMATICA
a[n_] := (MatrixPower[{{1, 3, 2}, {2, 1, 3}, {3, 2, 1}}, n].{{1}, {0}, {0}})[[1, 1]]; Table[ a[n], {n, 10}] (* Robert G. Wilson v, May 29 2004 *)
CROSSREFS
Sequence in context: A096913 A026916 A106173 * A189605 A000470 A297905
KEYWORD
nonn
AUTHOR
Gary W. Adamson, May 25 2004
EXTENSIONS
Edited and extended by Robert G. Wilson v, May 29 2004
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 20 11:46 EDT 2024. Contains 372712 sequences. (Running on oeis4.)