重新启动; nmax（最大值） :=21 20; 具有(组合):jn公司cfn1型 :=nmax（最大值）: 感应电动机 程序(n个, k个): 总和((-1)^j个*斯特林1(n个+1, n个+1-k个+j个) * 斯特林1(n个+1, n个+1-k个-j个), j个=nmax（最大值）-k个..k个) 结束 程序：欧米茄（0） := 1：对于从1到nmax的n do对于 j个 从 1 到 jn公司 做 cfn1型(1, j个):=1 结束 做: 对于 我 从 2 到 感应电动机 做 cfn1型(我, 1):=0 结束 做: 对于 j个 从 2 到 jn公司 做 对于 我 从 2 到 感应电动机 做 cfn1型(我,j个):=cfn1型(我-1,j个-1)*(j个-1)^2+cfn1型(我,j个-1) 结束 做 结束 做: 欧米茄（n） ：=（总和（（-1））^(k个千1+n+1）*（伯努利（2*k个千1)/(2*k个千1))*cfn1（n-k个+1个，n个-千1),k个千1=1..n））/（2*n-1）！end do：对于n，从1到nmax do d（n） := 2^(2*n个-1)*欧米茄（n）*2^(2*n个-1)) end do：对于n从2到nmax do Zc（n-1） ：=d（n-1）*2/（（2*n-1）*（n-1））端do:c（1） := denom（Zc（1））：对于从1到nmax-1的n，do c（n+1） ：=lcm（c（n）*（n+1）*（2*n+3）/2， 分母（Zc（n+1））：p（n+1 := c（n）end do：对于n，从2到nmax do一a1级（n） ：=p（n）*2^（2*n-3）/（3*阶乘（2*n-1））od:seq(一a1级（n） ，n=2..nmax）；
#结束第一个程序
重新启动; nmax（最大值）:=21; jn公司:=nmax（最大值）+1: 感应电动机n最大值1 := nmax（最大值）+1：对于来自的n1 到 nmax（最大值） 做 对于 我 从 2 0 到感应电动机 n最大值1 进行cfn2(我, 1):=n个, 0 结束 做) : 对于 j个 从 = 1 到 jn公司 做 : cfn2（立方英尺/平方英寸）(1,j个n个, n个) := (双阶乘(2*n个-1 结束 做))^2 日：用于j个 n个 从2 1 到jn公司 n最大值1 为…做我 k个 从2 1 到感应电动机 n个-1 进行cfn2(我,j个n个, k个) :=(2*n个-1)^2*cfn2（立方英尺/平方英寸）(我,j个n个-1, k个-1） +cfn2(我n个-1,j个- k个) 日: 日: 对于 n个 从 1)*(2*j个-三)^2 结束 做 结束 到 n最大值1 做: 增量（n-1） ：=总和（（1-2^（2*k个千1-1） ）*（-1）^（n+1）*（-bernoulli（2*k个千1)/(2*k个千1))*(-1)^(k个千1+n） *cfn2（n-k个+1，n-千1),k个千1=1..n）/（2*4^（n-1）*（2*n-1）！）；兰巴达（-2， n） ：=总和（2*（1-2^（2*k个千1-1） ）*（-伯努利（2*k个千1)/ (2*k个千1))*(-1)^(k个千1+n） *cfn2（n+-1-k个,,n个-千1),k个 千1=1..n）/阶乘（2*n-2）末端do:Lcgz（2） := 1/12:f（2） := 1/12:n从3到nmax（最大值） n最大值1 做Lcgz（n） := 兰巴达（-2， n-1）/（（2*n-2）*（2*n-3））：f（n） ：=Lcgz（n）-（（2*n-3）/（2*n-2））*f（n-1）end do：对于n从1到nmax（最大值） n最大值1 do b（n） := denom（Lcgz（n+1））end do：对于从1到的nnmax（最大值） n最大值1 do b（n） := 2*n*denom（Delta（n-1））/2^（2*n）end do:p（2） := b（1）：n从2到nmax（最大值） n最大值1 do p（n+1） ：=lcm（p（n）*（2*n）*， b（n））结束do：对于n，从2到nmax（最大值） n最大值1 做一a2类（n） := p（n）/（6*阶乘（2*n-2））od:seq(一a2类（n） ，n=2。。nmax（最大值）n最大值1);
#结束第二个程序
