Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A123971
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A123971 tridiagonal matrices of the type 3 center with terminal 2 as characteristic polynomial triangular sequence. +0
1
2, 2, -1, 5, -5, 1, 13, -19, 8, -1, 34, -65, 42, -11, 1, 89, -210, 183, -74, 14, -1, 233, -654, 717, -394, 115, -17, 1, 610, -1985, 2622, -1825, 725, -165, 20, -1, 1597, -5911, 9134, -7703, 3885, -1203, 224, -23, 1, 4181, -17345, 30691, -30418, 18633, -7329, 1855, -292, 26, -1, 10946, -50305, 100284, -114242 (list; table; graph; listen)
OFFSET

1,1

COMMENT

Major roots are odd Tutte-Beraha constants plus one: 2., 3.61803, 4.24698, 4.53209, 4.68251, 4.77091, 4.82709, 4.86494, 4.89163, 4.91115} Matrices: {{2}}, {{3, -1}, {-1, 2}}, {{3, -1, 0}, {-1, 3, -1}, {0, -1, 2}}, {{3, -1, 0, 0}, {-1, 3, -1, 0}, {0, -1, 3, -1}, {0, 0, -1, 2}}, {{3, -1, 0, 0, 0}, {-1, 3, -1, 0, 0}, {0, -1, 3, -1, 0}, {0, 0, -1, 3, -1}, {0, 0, 0, -1, 2}}, {{3, -1, 0, 0, 0, 0}, {-1, 3, -1, 0, 0, 0}, {0, -1, 3, -1, 0, 0}, {0, 0, -1, 3, -1, 0}, {0, 0, 0, -1, 3, -1}, {0, 0, 0, 0, -1, 2}} The determinants are increasing: {2, 5, 13, 34, 89, 233, 610, 1597, 4181, 10946}

FORMULA

M(n,m,d)=If[ n == m && n < d && m < d, k, If[n == m - 1 || n == m + 1, -1, If[n == m == d, k - 1, 0]]];

EXAMPLE

Triangular sequence:

{2},

{2, -1},

{5, -5, 1},

{13, -19, 8, -1},

{34, -65, 42, -11,1},

{89, -210, 183, -74, 14, -1}

{233, -654, 717, -394, 115, -17, 1}

MATHEMATICA

Mathematica ( general k th center) Clear[M, T, d, a, x, k] k = 3 T[n_, m_, d_] := If[ n == m && n < d && m < d, k, If[n == m - 1 || n == m + 1, -1, If[n == m == d, k - 1, 0]]] M[d_] := Table[T[n, m, d], {n, 1, d}, {m, 1, d}] Table[M[d], {d, 1, 10}] Table[Det[M[d]], {d, 1, 10}] Table[Det[M[d] - x*IdentityMatrix[d]], {d, 1, 10}] a = Join[{M[1]}, Table[CoefficientList[ Det[M[d] - x*IdentityMatrix[d]], x], {d, 1, 10}]] Flatten[a] MatrixForm[a] Table[NSolve[Det[M[d] - x*IdentityMatrix[d]] == 0, x], {d, 1, 10}] Table[x /. NSolve[Det[M[d] - x*IdentityMatrix[d]] == 0, x][[d]], {d, 1, 10}]

CROSSREFS

Adjacent sequences: A123968 A123969 A123970 this_sequence A123972 A123973 A123974

Sequence in context: A099605 A079218 A079220 this_sequence A114292 A141751 A079222

KEYWORD

uned,probation,tabl,sign

AUTHOR

Gary Adamson and Roger Bagula (rlbagulatftn(AT)yahoo.com), Oct 30 2006

EXTENSIONS

Looking at the triangle suggests that the very first term should be 1, not 2. - njas, Nov 01, 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 October 10 20:39 EDT 2008. Contains 144831 sequences.


AT&T Labs Research