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