（PARI）/*对于此序列，运行limnpq（nmax，1，2））：*/
limnpq（nmax，p，q）={
f=矢量（nmax+1）；f[1]=q/（p+q）；r=-p/（p+q）；
对于（i=2，nmax+1，p1=i-1；bc=1；m=p1；s=0；
对于（j=1，i-1，p2=j-1；如果（p2，bc=bc*m/p2；m=m-1；）；
s=s+bc*f[j]；）；f[i]=r*s；）；
fac=（p+q）/q；
对于（i=1，nmax+1，f[i]=f[i]*fac；fac=（p+q）*fac；
写入（“输出文件”，i-1，“”，f[i]）；）；}
（PARI）x='x+O（'x^66）；Vec（塞拉普拉斯（3/（2+exp（3*x）））\\乔格·阿恩特2013年4月21日
（PARI）/*O.g.f.：*/
{a（n）=polcoeff（总和（m=0，n，m！*（-x）^m/prod（k=1，m，1-3*k*x+x*O（x^n）），n）}\\保罗·D·汉纳2013年5月30日
