位数=96;F[x_x]:=1-Exp[-x]/Sqrt[Pi*x]-Erf[Sqrt[x]];清除[f,g];
f[m_]:=f[m]=2个整数[(1-(1-f[x])^-1),{x,0,m},工作精度->数字+10];f[m=100];f【m=2米】;打印[“m=”,m];而[realdights[f[m],10,位数+5][[1]]!=实数位数[f[m/2],10,位数+5][[1]],m=2m;打印[“m=,m]];
g[m_]:=g[m]=(8/3)n整体[(1-(1-F[x])^-1)*x,{x,0,m},工作精度->数字+10];克[m=100];g[米=2米];打印[“m=”,m];而[realdights[g[m],10,位数+5][[1]]!=实数[g[m/2],10,位数+5][[1]],m=2m;打印[“m=,m]];
加入[{0},实数位数[g[m]-f[m]^2,10,数字][[1]]]
|