Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- Clear[s, sum];
- nn = 12;
- A = Table[Table[If[Mod[n, k] == 0, 1/n, 0], {k, 1, nn}], {n, 1, nn}];
- B = Table[
- Table[If[Mod[k, n] == 0, MoebiusMu[n]*n/k^s, 0], {k, 1, nn}], {n,
- 1, nn}];
- sum = Total[Total[A.B]];
- Im[s /. Table[FindRoot[sum == 0, {s, ZetaZero[n]}], {n, 1, 11}]]
- N[Im[ZetaZero[Range[11]]]]
- % - %%
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement