|
PROG
|
(PARI) {a(n)=local(A=x+x^2); for(i=1, n, A=x+subst(A, x, subst(A, x, A^2+x*O(x^n)))); polcoeff(A, n)}
(PARI) {a(n)=local(A=x); if(n<1, 0, for(i=1, n, A=serreverse(x - subst(A, x, subst(A, x, x^2+x*O(x^n))) )); polcoeff(A, n))}
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+x^2+x*O(x^n)); for(i=1, n, A=x+sum(m=1, n, Dx(m-1, subst(A, x, subst(A, x, x^2+x*O(x^n)))^m)/m!)+x*O(x^n)); polcoeff(A, n)}
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+x^2+x*O(x^n)); for(i=1, n, A=x*exp(sum(m=1, n, Dx(m-1, subst(A, x, subst(A, x, x^2+x*O(x^n)))^m/x)/m!)+x*O(x^n))); polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))
|