login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A136455
Characteristic polynomials of the inverse beta function based matrices as a triangle of integer coefficients: n*IM(i,j) = n*Inverse(1/Gamma(i,j)); i,j>=n.
1
1, 0, 1, -1, 1, 1, -48, 28, 25, -1, 233280, -91368, -60993, 2305, 1, 222953472000, -65503641600, -33198846720, 985867696, 446161, -1, -69132994560000000000, 16249035196800000000, 6593300559405000000, -157196644177875000, -59060479175425, 144069601, 1
OFFSET
1,7
LINKS
FORMULA
M(i,j)=1/Gamma[i+j]; i,j<=n IM(i,j)=Inverse(M(i,j))
EXAMPLE
{1},
{0, 1},
{-1, 1, 1},
{-48, 28, 25, -1},
{233280, -91368, -60993, 2305, 1},
{222953472000, -65503641600, -33198846720, 985867696, 446161, -1}
MATHEMATICA
M[w_] := Table[Table[1/Gamma[n + m], {n, 0, w}], {m, 0, w}]
IM[w_] := Inverse[M[w]] Join[{1, x}, Table[CharacteristicPolynomial[n*IM[n], x], {n, 1, 10}]]
a = Join[{{1}, {0, 1}}, Table[CoefficientList[CharacteristicPolynomial[n*IM[n], x], x], {n, 1, 10}]];
Flatten[a]
Join[{1, 1}, Table[Apply[Plus, CoefficientList[CharacteristicPolynomial[n*IM[n], x], x]], {n, 1, 10}]]
CROSSREFS
Sequence in context: A128380 A298451 A094658 * A085517 A299525 A033694
KEYWORD
uned,tabl,sign
AUTHOR
Roger L. Bagula, Mar 20 2008
STATUS
approved