|
PROG
|
(PARI) {a(n) = my(A=[1]); for(i=1, n, A=concat(A, 0);
A[#A] = polcoeff(-1 + sum(m=-#A, #A, x^m * (2*Ser(A) - x^m)^(4*m-1) ), #A-1)/2); A[n+1]}
for(n=0, 30, print1(a(n), ", "))
(PARI) {a(n) = my(A=[1]); for(i=1, n, A=concat(A, 0);
A[#A] = polcoeff(1 - sum(m=-#A, #A, x^(4*m^2)/(1 - 2*Ser(A)*x^m)^(4*m+1) ), #A-1)/2); A[n+1]}
for(n=0, 30, print1(a(n), ", "))
|