|
n=30；（*来自Rains和Sloane的算法*）
S4[f，h，x_]：=f[h，x]^4/24+f[h、x]^2f[h和x^2]/4+f[h，x]f[h；x^3]/3+f[w，x^2]^2/8+f[n，x^4]/4；
T[-1，z_]：=1；T[h_，z_]：=T[h，z]=表[z^k，{k，0，n}]。取[系数表[z^（n+1）+1+S4[T，h-1，z]z，z]，n+1]；
求和[取[系数列表[z^（n+1）+（T[h，z]-T[h-1，z]）^2/2+（T[h，z^2]-T[h-1,z^2]）/2，z]，n+1]，{h，0，n/2}](*罗伯特·拉塞尔2018年9月15日*）
| 