(PARI){RR(n,w,z,p,po,i,m,h,h1,j,w1,h2)=w=1+O(x^(n+1));p=1;po=1;对于(i=1,n,w=p-po*x*q^i;po=p;p=w);m=极度数(w);w1=0;对于(i=0,m,h=polceoff(w,i);h1=0,对于(j=1,n-1+i,h1=h1+polcooff(h,j)*q^j);w1=w1+h1*x^i);q=0;对于(i=1,n-1,q=q+s[i]/x^i*polceoff(h2,0)*(-1)}s=向量(30);s[1]=1;打印(s[1]);对于(j=2,30,s[j]=RR(j);打印(s[j]);\\Kok Seng Chua(chuaks(AT)ihpc.nus.edu.sg),2000年6月13日