|
|
A324313
|
|
Expansion of the unique formal power series R(t) with constant term 0 satisfying t = Sum_{n>=0} (1/(n+1))*binomial(2n,n)*binomial(3n,n)/*R(t)^(n+1).
|
|
3
|
|
|
1, -3, -12, -105, -1206, -16002, -232416, -3592377, -58113000, -973372686, -16757906256, -295003183410, -5289613240464, -96324638802300, -1777361586269760, -33170533400455353, -625219596117283590, -11887661990007138006, -227778723332612670600, -4394621278195058676150
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
MATHEMATICA
|
m = 20;
aa = Array[a, m];
R[t_] = aa.t^Range[m];
eq = Thread[CoefficientList[t-Sum[1/(n+1) Binomial[2n, n] Binomial[3n, n]* R[t]^(n+1), {n, 0, m}] + O[t]^(m+1), t] == 0];
|
|
PROG
|
(PARI) lista(nn) = {my(v = vector(nn), R, P, c, r, s); kill(y); for (n=1, nn, v[n] = y; R = sum(k=1, n, v[k]*t^k); P = sum(k=0, n, binomial(2*k, k)*binomial(3*k, k)/(k+1)*R^(k+1)); c = polcoef(P, n, t); r = polcoef(c, 0, y); s = polcoef(c, 1, y); if (n==1, v[n] = (1-r)/s, v[n] = -r/s); ); R = sum(k=1, #v, v[k]*t^k); vector(nn, k, polcoef(R, k, t)); }
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|