数字=1004;
nmax0=50;(*总和术语的初始编号*)
dnmax=10;(*nmax增量*)
dd=10;(*精度超标*)
清除[PrimeZeta41];
f[s_]:=(1+2^-s)^-1*DirichletBeta[s]Zeta[s]/Zeta[2s];
PrimeZeta41[s_,nmax_]:=PrimeZeta 41[s,nmax]=(1/2)总和[MoebiusMu[2 n+1]*Log[f[(2 n+1)*9]/(2 n+1),{n,0,nmax}]//n[#,digits+dd&;
PrimeZeta41[9,nmax=nmax0];
PrimeZeta41[9,nmax+=dnmax];
当[Abs[PrimeZeta41[9,nmax]-PrimeZeta41[9,nmax-dnmax]]>10^-(数字+dd)时,打印[“nmax=”,nmax];nmax+=dnmax];
PrimeZeta41[9]=PrimeZeta 41[9,nmax];
联接[{0,0,0,10,0},RealDigits[PrimeZeta41[9],10,digits][[1](*Jean-François Alcover公司2011年6月24日,继X.Gourdon和P.Sebah之后,于2021年5月7日更新*)
|