Correction: new procedure for converting a new recursive polynomial set into matrices

• To: mathgroup at smc.vnet.net
• Subject: [mg70731] Correction: new procedure for converting a new recursive polynomial set into matrices
• From: Roger Bagula <rlbagula at sbcglobal.net>
• Date: Wed, 25 Oct 2006 01:39:56 -0400 (EDT)
• References: <ehkcmn\$jn9\$1@smc.vnet.net>

```I was getting a "shift" function type of output.
It took some work to get the numbers lined up right.
The signs in the two triangular sequences still come up different!

Clear[p, a, b]
p[0, x] = 1; p[1, x] = x - 1;
p[k_, x_] := p[k, x] = p[k - 1, x] - x^2*p[k - 2, x]
Table[Expand[p[n, x]], {n, 0, 10}]
Table[Sum[CoefficientList[p[n, x], x][[m]], {m, 1, n + 1}], {n, 0, 20}]
w = Table[CoefficientList[p[n, x], x], {n, 0, 20}]
Flatten[w]
MatrixForm[w]
An[d_] := Table[If[n == d, -w[[n]][[m]], If[m == n, 1, 0]], {n, 2, d}, {m,
1, d - 1}]
Table[An[d], {d, 2, 19}]
Table[CharacteristicPolynomial[An[d], x], {d, 2, 19}]
b = Join[{{1}},
Table[CoefficientList[CharacteristicPolynomial[An[d], x], x], {d,
2, 19}]]
Flatten[%]
MatrixForm[b]

>
>
>

```

• Prev by Date: Re: ListPlot and ListPlot3D
• Next by Date: Re: Solve with dot products
• Previous by thread: Re: what's the command to load a specific "page" in the help browser?
• Next by thread: Defining two functions at once?