Advertisement
MatsGranvik

Fast converging roots of Riemann spectrum

Jan 14th, 2018
196
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.31 KB | None | 0 0
  1. Clear[s, sum];
  2. nn = 12;
  3. A = Table[Table[If[Mod[n, k] == 0, 1/n, 0], {k, 1, nn}], {n, 1, nn}];
  4. B = Table[
  5. Table[If[Mod[k, n] == 0, MoebiusMu[n]*n/k^s, 0], {k, 1, nn}], {n,
  6. 1, nn}];
  7. sum = Total[Total[A.B]];
  8.  
  9. Im[s /. Table[FindRoot[sum == 0, {s, ZetaZero[n]}], {n, 1, 11}]]
  10. N[Im[ZetaZero[Range[11]]]]
  11. % - %%
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement