1,7

Based on:

Beta[n,m]=Gamma[n]*Gamma[m]/Gamma[n+m]=Integate[x^n&(1-x)^m,{x,0,1}];

f[x,n]=x^n/Gamma[n]

g[x,n]=(1-x)^n/Gamma[n]

Integral:

Matrix[n,m]=Integrate[f[x,n]*g[x,m],{x,0,1}]=1/Gamma[n,m]

IM[n]=n*Inverse[Matrix[n,m]]

These matrices are made to be like the transorthogonal or simplex coding :

-1/(2^n-1)

1/Gamma[n+m] is mostly less than that.

The row sums are:

1, 1, 1, 4, 83225, 125237297536, -46447914508307980823,

-7674877258831704425297015668284,

826995023940325227060455484795500305924239025,

78707026870454795917065346205715405165749494298851023916447744,

-8564339911597252330069530473576149288702875987041232549818479682054967832279248559, -1336289257117044917306250532620084989502533702227664929930431988712778911134146754043951495782103516647100

Weisstein, Eric W. "Beta Function." http : // mathworld.wolfram.com/BetaFunction.html

Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Mar 20 2008, Table of n, a(n) for n = 1..78

M(i,j)=1/Gamma[i+j];i,j<=n IM(i,j)=Inverse(M(i,j))

{1},

{0, 1},

{-1, 1, 1},

{-48, 28, 25, -1},

{233280, -91368, -60993, 2305, 1},

{222953472000, -65503641600, -33198846720, 985867696, 446161, -1}

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}]]

Sequence in context: A033368 A128380 A094658 this_sequence A085517 A033694 A156461

Adjacent sequences: A136452 A136453 A136454 this_sequence A136456 A136457 A136458

uned,tabl,sign

Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Mar 20 2008