(PARI)/*根据定义:*/
{a(n)=n!*polcoeff(和(k=0,n,(exp((k+1)*x+x*O(x^n))-1)^k),n)}
对于(n=0,25,print1(a(n),“,”)
(PARI)/*例如从无穷级数:*/
\p100\\设置精度
{A=Vec(serlaplace(总和(n=0,500,1.*exp(n*(n+1)*x+O(x^26))/(1+exp(n*x+O(x^26)
对于(n=0,#A-1,打印1(圆形(A[n+1]),“,”)
(PARI)/*来自o.g.f.:*/
{a(n)=极坐标(和(m=0,n,(m+1)^m*m!*x^m/prod(k=1,m,1-(m+1,*k*x+x*O(x^n)),n)}
对于(n=0,25,print1(a(n),“,”)
(PARI){斯特林2(n,k)=n!*polceoff(((exp(x+x*O(x^n))-1)^k)/k!,n)}
{a(n)=和(k=0,n,(k+1)^n*k!*Stirling2(n,k))}
对于(n=0,25,print1(a(n),“,”)