[seqfan] Re: A132870 : a Mathematica line

[Paul Hanna]

Wouter,
      Thanks for verifying the fractional elements.
"Bummer" is right ...
I'll ask Neil what to do with the non-integral triangle.
Thank you!
      Paul

On Sat, Nov 18, 2017 at 11:50 AM, Wouter Meeussen <
wouter.meeussen at telenet.be> wrote:

> pt = {{1}}; Table[rhs = CoefficientList[(k^2 + x)^k, x];
> qt = Join[pt, {vars = Array[Subscript[a, #] &, k + 1]}];
> b = MatrixPower[PadRight[qt], k] ;
> {out} = vars /. Solve[Thread[Reverse[Last[b]] == Reverse[rhs]], vars];
> pt = Join[pt, {out}]; out, {k, 20}]
>
>
> not the fastest of codes : the first 20 rows take 43 seconds on my old
> machine.
>
>
> W.
>
> --
> Seqfan Mailing list - http://list.seqfan.eu/
>