n=50;(*雷恩和斯隆的算法*)
S2[f,h,x]:=f[h,x]^2/2+f[h,x^2]/2;
S3[f,h,x^3]:=f[h,x]^3/6+f[h,x]f[h,x^2]/2+f[h,x^3]/3;
T[-1,z_9]:=1;T[h,z}:=T[h,z]=表[z^k,{k,0,n}]。取[CoefficientList[z^(n+1)+1+S2[T,h-1,z]z,z],n+1];
Sum[取[CoefficientList[z^(n+1)+S3[T,h-1,z]z-S3[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[{1,1},n+1](*罗伯特A.罗素2018年9月15日*)
|