Advertisement
MatsGranvik

Newton Raphson iteration applied to Riemann zeta function

Aug 11th, 2017
248
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.73 KB | None | 0 0
  1. (*Newton Raphson iteration*)
  2. (*start*)
  3. f[t_] = Zeta[1/2 + I*t]/Zeta'[1/2 + I*t];
  4. g[t_] = Im[(1/2 + I*t) - f[t]];
  5. min = 0;
  6. max = 60;
  7. Show[Plot[g[g[g[g[g[g[g[g[g[g[g[t]]]]]]]]]]], {t, min, max}],
  8. Plot[t, {t, min, max}],
  9. Table[Graphics[{PointSize[Large],
  10. Point[{Im[ZetaZero[n]], Im[ZetaZero[n]]}]}], {n, 1, 14}]]
  11. (*end*)
  12.  
  13. (* variant of logarithmic derivative Newton Raphson iteration*)(*start*)
  14. c = 1 + 1/1000;
  15. f[t_] = -1/((Zeta[1/2 + I*t]*Zeta[c])/Zeta[1/2 + I*t + c - 1] -
  16. Zeta[c]);
  17. g[t_] = Im[(1/2 + I*t) - f[t]];
  18. min = 0;
  19. max = 60;
  20. Show[Plot[g[g[g[g[g[g[t]]]]]], {t, min, max}], Plot[t, {t, min, max}],
  21. Table[Graphics[{PointSize[Large],
  22. Point[{Im[ZetaZero[n]], Im[ZetaZero[n]]}]}], {n, 1, 14}]]
  23. (*end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement