OFFSET
0,4
MAPLE
A := 1; f := proc(n) global A; local A2, A3; A2 := subs(x=x^2, A); A3 := subs(x=x^3, A);
coeff(series( 1+x*( (A^4+3*A2^2+8*A*A3)/12), x, n+1), x, n); end;
for n from 1 to 50 do A := series(A+f(n)*x^n, x, n +1); od: A;
MATHEMATICA
a = 1; f[n_] := Module[{a2, a3}, a2 = a /. x -> x^2; a3 = a /. x -> x^3; Coefficient[ Series[1 + x*(a^4 + 3*a2^2 + 8*a*a3)/12, {x, 0, n + 1}] // Normal, x, n]]; For[n = 1, n <= 30, n++, a = Series[a + f[n]*x^n, {x, 0, n + 1}] // Normal]; CoefficientList[a, x] (* Jean-François Alcover, Jan 16 2013, after Maple *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved