(PARI){a(n)=局部(Egf=1,X=X+X*O(X^n),R=serreverse(X/(exp(3*X)*cosh(3*X)));Egf=exp(R)*coph(R);n!*polcoff(Egf,n)}
对于(n=0,25,print1(a(n),“,”)
(PARI)/*从LambertW恒等式导出的公式:*/
{a(n)=局部(Egf=1,X=X+X*O(X^n));Egf=和(k=0,n,(3*k+1)^(k-1)*cosh((3*k+1)*X)*X^k/k!);n!*polcoeff(Egf,n)}
对于(n=0,25,print1(a(n),“,”)