|
PROG
|
(PARI) {a(n)=local(A=1+x, B); for(i=0, n, A = 1 + x*A * subst(A, x, x*A+x*O(x^n))); B=A; B=subst(B, x, x*A+x*O(x^n)); polcoeff(B, n)}
for(n=0, 30, print1(a(n), ", "))
(PARI) {a(n)=local(A=1+x+x*O(x^n)); for(i=1, n, A = 1/subst(1-x*A, x, x/(1-x*A +x*O(x^n))) ); polcoeff(A, n)}
for(n=0, 30, print1(a(n), ", "))
|