Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- (*Mathematica start*)
- nn = 200;
- f[t_] = D[RiemannSiegelTheta[t], t];
- g2 = Plot[
- Sign[RiemannSiegelZ[t]]*
- Abs[Zeta[1/2 + I*t]*
- Total[Table[
- Total[MoebiusMu[Divisors[n]]/Divisors[n]^(1/2 + I*t - 1)]/
- n, {n, 1, nn}]]]/(f[t] + HarmonicNumber[nn]), {t, 0, 60},
- PlotStyle -> Thickness[0.004], ImageSize -> Large,
- PlotRange -> {-4, 4}]
- (*end*)
- (*Mathematica analytic continuation start*)
- f[t_] = D[RiemannSiegelTheta[t], t];
- epsilon = 1/40;
- Monitor[Plot[
- Sign[RiemannSiegelZ[t]]*
- Abs[Zeta[1/2 + I*t]/Zeta[1/2 + I*t + (1 + epsilon) - 1]*
- Zeta[(1 + epsilon)]]/(f[t] + 1/epsilon + EulerGamma), {t, 0,
- 60}, PlotStyle -> Thickness[0.004], ImageSize -> Large,
- PlotRange -> {-4, 4}], n]
- (*end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement