Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A123966
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A123966 Jacobi tridiagonal matrices with 4 on the main diagonal used to get characteristic polynomials and a triangular sequence. +0
2
4, 4, -1, 15, -8, 1, 56, -46, 12, -1, 209, -232, 93, -16, 1, 780, -1091, 592, -156, 20, -1, 2911, -4912, 3366, -1200, 235, -24, 1, 10864, -21468, 17784, -8010, 2120, -330, 28, -1, 40545, -91824, 89238, -48624, 16255, -3416, 441, -32, 1, 151316, -386373, 430992, -275724, 111524, -29589, 5152, -568, 36 (list; table; graph; listen)
OFFSET

1,1

COMMENT

Matrices: {{4}}, {{4, -1}, {-1, 4}}, {{4, -1, 0}, {-1, 4, -1}, {0, -1, 4}}, {{4, -1, 0, 0}. {-1, 4, -1, 0}, {0, -1, 4, -1}, {0, 0, -1, 4}}, {{4, -1, 0, 0, 0}, {-1, 4, -1, 0, 0}, {0, -1, 4, -1, 0}, {0, 0, -1, 4, -1}, {0, 0, 0, -1, 4}}, {{4, -1, 0, 0, 0, 0}, {-1, 4, -1, 0, 0, 0}, {0, -1, 4, -1, 0, 0}, {0, 0, -1, 4, -1, 0}, {0, 0, 0, -1, 4, -1}, {0, 0, 0, 0, -1, 4}}

REFERENCES

Eric Weisstein's World of Mathematics, "Tridiagonal Matrix." http://mathworld.wolfram.com/TridiagonalMatrix.html

FORMULA

a(n,m)=If[ n == m, 4, If[n == m - 1 || n == m + 1, -1, 0]]] p(n,x)=CharacteristicPolynomial(a(n,m)) p(n,x)->t(n,m)

EXAMPLE

Polynomials:

4

4 - x,

15 - 8 x + x^2,

56 - 46 x + 12 x^2 - x^3,

209 - 232 x + 93 x^2 - 16 x^3 + x^4,

780 - 1091 x + 592 x^2 - 156 x^3 + 20 x^4 - x^5,

2911 - 4912 x + 3366 x^2 - 1200 x^3 + 235 x^4 - 24 x^5 + x^6

Triangular sequence:

{4},

{4, -1},

{15, -8, 1},

{56, -46,12, -1},

{209, -232, 93, -16, 1},

{780, -1091, 592, -156, 20, -1},

{2911, -4912, 3366, -1200, 235, -24, 1},

{10864, -21468, 17784, -8010, 2120, -330, 28, -1}

MATHEMATICA

Clear[M, T, d, a, x]; T[n_, m_] = If[ n == m, 4, If[n == m - 1 || n == m + 1, -1, 0]]; M[d_] := Table[T[n, m], {n, 1, d}, {m, 1, d}]; Table[M[d], {d, 1, 10}]; Table[Det[M[d] - x*IdentityMatrix[d]], {d, 1, 10}]; a = Join[{{3}}, Table[CoefficientList[Det[M[d] - x*IdentityMatrix[d]], x], {d, 1, 10}]]; Flatten[a]

CROSSREFS

Cf. A123343.

Sequence in context: A102323 A145902 A124028 this_sequence A079507 A098364 A116866

Adjacent sequences: A123963 A123964 A123965 this_sequence A123967 A123968 A123969

KEYWORD

uned,probation,tabl,sign

AUTHOR

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

EXTENSIONS

Looking at the triangle suggests that the very first term should be 1, not 4. - N. J. A. Sloane (njas(AT)research.att.com), 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 December 16 17:18 EST 2009. Contains 170825 sequences.


AT&T Labs Research