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