Advertisement
MatsGranvik

Animated gif zeta zero count

Jun 23rd, 2018
377
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.55 KB | None | 0 0
  1. (*start*)
  2. Clear[n, s, c, t, nn, delta, k];
  3. c = 1;
  4. delta = 20;
  5. nn = 60;
  6. kk = 120;
  7. Monitor[Table[
  8. a = Sum[
  9. 1/N[n]^c -
  10. Zeta[1/2 + I*t]*
  11. Total[MoebiusMu[Divisors[n]]/Divisors[n]^(1/2 + I*t - 1)]/
  12. n^c, {n, 1, k}];
  13. Export[StringJoin["plot", ToString[k], ".gif"],
  14. ListLinePlot[(Table[
  15. delta*RiemannSiegelTheta[t], {t, 0, nn, 1/delta}] +
  16. Accumulate[Table[Re[a], {t, 0, nn, 1/delta}]])/delta/Pi,
  17. PlotStyle -> Thickness[0.003], DataRange -> {0, nn},
  18. ImageSize -> Large]], {k, 1, kk}], k]
  19. (*end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement