(PARI)/*A(x)=x+(p+q)*G((p*A(x)+q*x)/(p+q));G(x)=积分A(x)dx:*/
{a(n,p=2,q=-1)=my(a=x,G)
对于(n=1,30,打印1(a(n,2,-1),“,”)
(PARI)/*A(x-积分p*A(x)dx)=x+积分q*A(x)dx:*/
{a(n,p=2,q=-1)=my(a=[1],F=x);对于(i=1,n,a=concat(a,0);F=x*Ser(a);G=intformal(F);a[#a]=-polcoeff(subst(F,x,x-p*G)-q*G,#a));n!*a[n]}
对于(n=1,30,打印1(a(n,2,-1),“,”)
(PARI)/*生成前N个术语的非正式代码:*/
{N=20;p=2;q=-1;A=x;对于(i=1,N,G=整数(A+x*O(x^N));A=x+(p+q)*子集(G,x,(p*A+q*x)/(p+q));Vec(serlaplace(A))}