|
PROG
|
(PARI) {a(n)=local(A=x^2, G=x); for(i=0, n, A=intformal(2*G); G = serreverse(x - 3*A^2*G +O(x^(4*n)))); polcoeff(A, 4*n-2)}
for(n=1, 25, print1(a(n), ", "))
(PARI) {a(n)=local(A=x^2, G=x); for(i=1, n, A=intformal(2*G); G = serreverse(x - 3*A^2*G +O(x^(4*n)))); polcoeff(A, 4*n-2)}
for(n=1, 25, print1(a(n), ", "))
(PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}
{a(n)=local(A=x^2); for(i=1, n, A = 2*intformal(x + sum(m=1, n+1, Dx(m-1, (3/2)^m*A^(2*m)*(A')^m/m!)) +O(x^(4*n)))); polcoeff(A, 4*n-2)}
for(n=1, 25, print1(a(n), ", "))
(PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}
{a(n)=local(A=x^2); for(i=1, n, A = 2*intformal(x*exp(sum(m=1, n, Dx(m-1, (3/2)^m*A^(2*m)*(A')^m/(m!*x))) +O(x^(4*n))))); polcoeff(A, 4*n-2)}
for(n=1, 25, print1(a(n), ", "))
|