Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A124020
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A124020 New tetradiagonal form matrix as triangular sequence from solution of : X(n,m)=Steinbach(n,m)^(-1).tri-Antidiagonal_1(n,n). +0
1
1, 1, -1, 0, -3, 1, -1, -4, 5, -1, -1, -2, 12, -7, 1, 0, 2, 17, -24, 9, -1, 1, 4, 13, -52, 40, -11, 1, 1, 1, 0, -76, 115, -60, 13, -1, 0, -5, -10, -72, 235, -214, 84, -15, 1, -1, -8, -2, -34, 352, -554, 357, -112, 17, -1, -1, -4, 24, 2, 383, -1092, 1113, -552, 144, -19, 1, 0, 4, 46, -24, 295, -1673, 2688, -2008, 807, -180, 21, -1 (list; table; graph; listen)
OFFSET

1,5

COMMENT

Matrices {{1}}, {{1, -1}, {-2, 2}}, {{1, -1, 0}, {-2, 2, -1}, {1, -2, 2}}, {{1, -1, 0, 0}, {-2, 2, -1, 0}, {1, -2, 2, -1}, {0, 1, -2, 2}}, {{1, -1, 0, 0, 0}, {-2, 2, -1, 0, 0}, {1, -2, 2, -1, 0}, {0, 1, -2, 2, -1}, {0, 0, 1, -2, 2}} Large root sequence: a0 = Table[x /. NSolve[CharacteristicPolynomial[Xn[d], x] == 0, x][[d]], {d, 1, 20}] {1., 3., 3.91223, 4.37167, 4.62826, 4.78478, 4.88683, 4.95691, 5.00703, 5.0441, 5.07226, 5.09415, 5.1115, 5.12547, 5.1369, 5.14636, 5.15428, 5.16097, 5.16668, 5.17159} Determinant sequence is: Table[Det[Xn[d]], {d, 1, 20}] {1,0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0}

FORMULA

Steinbach(n,m)=If[n + m - 1 > d, 0, 1]; tri-Antidiagonal_1(n,m)=If[n + m - 1 == d, 1, If[n + m == d, -1, If[n + m - 2 == d, -1, 0]]]; X(n,m)=Steinbach(n,m)^(-1).tri-Antidiagonal2(n,n)

EXAMPLE

triangular sequence:

{1},

{1, -1},

{0, -3, 1},

{-1, -4,5, -1},

{-1, -2, 12, -7, 1},

{0, 2, 17, -24, 9, -1},

{1, 4,13, -52, 40, -11,1},

{1, 1, 0, -76, 115, -60, 13, -1},

{0, -5, -10, -72, 235, -214, 84, -15,1}

MATHEMATICA

An[d_] := Table[If[n + m - 1 > d, 0, 1], {n, 1, d}, {m, 1, d}]; Bn[d_] := Table[If[n + m - 1 == d, 1, If[n + m == d, -1, If[n + m - 2 == d, -1, 0]]], {n, 1, d}, {m, 1, d}]; Xn[d_] := MatrixPower[An[d], -1].Bn[d]; a = Join[Xn[1], Table[CoefficientList[CharacteristicPolynomial[Xn[d], x], x], {d, 1, 20}]]; Flatten[%]

CROSSREFS

Sequence in context: A026374 A102716 A134510 this_sequence A124234 A135226 A104730

Adjacent sequences: A124017 A124018 A124019 this_sequence A124021 A124022 A124023

KEYWORD

sign,tabl,uned,probation

AUTHOR

Roger Bagula and Gary Adamson (rlbagulatftn(AT)yahoo.com), Oct 31 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 August 19 23:53 EDT 2008. Contains 142930 sequences.


AT&T Labs Research