(PARI){a(n)=局部(Egf=1,X=X+X*O(X^n),R=serreverse(X/(exp(4*X)*cosh(4*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,(4*k+1)^(k-1)*cosh((4*k+1)*X)*X^k/k!);n!*polcoeff(Egf,n)}
对于(n=0,25,打印1(a(n),“,”)