|
MATHEMATICA
|
nmax = 23; co = Coefficient; ex = Exponent;
b[n_, i_] := b[n, i] = If[n == 0, {0}, If[i<1, {}, Flatten[Table[Function[ {p}, p + j x^i] /@ b[n - i j, i - 1], {j, 0, n/i}]]]];
g[n_, k_] := g[n, k] = Sum[Sum[2^Sum[Sum[GCD[i, j] co[s, x, i] co[t, x, j], {j, 1, ex[t, x]}], {i, 1, ex[s, x]}]/Product[i^co[s, x, i]*co[s, x, i]!, {i, 1, ex[s, x]}]/Product[i^co[t, x, i] co[t, x, i]!, {i, 1, ex[t, x]}], {t, b[n + k, n + k]}], {s, b[n, n]}];
A[n_, k_] := g[Min[n, k], Abs[n - k]];
A[d_] := Sum[A[n, d - n], {n, 0, d}];
B[x_] = Sum[A[n] x^n, {n, 0, nmax}];
S[_, _] = 0; Do[S[c_, t_] = Series[1 + (c/(1 + c)) S[c, t]^2 + t S[c, t]^3, {c, 0, nmax}, {t, 0, nmax}] // Normal, {nmax}];
T[x_] = 1 - S[x/(1 - x), 1 - 2x - 1/B[x]];
|