S[m_,n_,S]:=(t=1;总和=0;difs=1;
当[Abs[difs]>10^(-数字-5)||difs==0时,
difs=(莫比乌斯亩[t]/t)*
对数[如果[s*t==1,DirichletL[m,n,s*t],
求和[Zeta[s*t,j/m]*DirichletCharacter[m,n,j]^t,{j,1,m}]/
m^(s*t)]];总和=总和+difs;t++];总额);
P[m_,n_,s]:=1/EulerPhi[m]*
求和[Conjugate[DirichletCharacter[m,r,n]]*S[m,r,S],{r,1,
欧拉菲[m]}]+
总和[如果[GCD[p,m]>1&&Mod[p,m]==n,1/p^s,0],{p,1,m}];
Z[m_,n_,s]:=(w=1;总和=0;difz=1;
当[Abs[difz]>10^(-数字-5)时,difz=P[m,n,s*w]/w;
sumz=sumz+difz;打印临时[w];w++];实验[sumz]);
(*使用Vaclav Kotesovec的函数ZA301430型. *)
$MaxExtraPrecision(最大额外精度)=1000100; 数字=121;50; (*调整 作为 需要. *)
数字化[c(c)_] :=真数字[Chop[N[Z轴[5,3c(c),2],数字+10]],10,位数-1][[1]] (*之后_瓦茨拉夫 科特索维奇_ *)]];
数字化[Z[5,3,2]]
|