Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- (*start*)
- Clear[n, s, c, t, nn, delta, k];
- c = 1;
- delta = 20;
- nn = 60;
- kk = 120;
- Monitor[Table[
- a = Sum[
- 1/N[n]^c -
- Zeta[1/2 + I*t]*
- Total[MoebiusMu[Divisors[n]]/Divisors[n]^(1/2 + I*t - 1)]/
- n^c, {n, 1, k}];
- Export[StringJoin["plot", ToString[k], ".gif"],
- ListLinePlot[(Table[
- delta*RiemannSiegelTheta[t], {t, 0, nn, 1/delta}] +
- Accumulate[Table[Re[a], {t, 0, nn, 1/delta}]])/delta/Pi,
- PlotStyle -> Thickness[0.003], DataRange -> {0, nn},
- ImageSize -> Large]], {k, 1, kk}], k]
- (*end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement