n=30；（*来自Rains和Sloane的算法*）
S4[f_，h，x_]：=f[h，x]^4/24+f[h，x]^2 f[h，x^2]/4+f[h，x]f[h，x^3]/3+f[h，x^2]^2/8+f[h，x^4]/4；
S5[f，h，x_]：=f[h，x]^5/120+f[h、x]^3f[h和x^2]/12+f[h，x]*2f[h；x^3]/6+f[h[h，x]f[h；
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]；
求和[Take[CoefficientList[z^（n+1）+S5[T，h-1，z]z-S5[T，h-2，z]z-（T[h-1，z]-T[h-2，z]）（T[h-1，z]-1），z]，n+1]，{h，1，n/2}]+PadRight[{0，1}，n+1）(*罗伯特·拉塞尔2018年9月15日*）
