MatsGranvik

Sqrt[MoebiusMu[k]*k]

Aug 8th, 2022 (edited)
203
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.88 KB | None | 0 0
  1. (*start*)
  2. nn = 200;
  3. s = N[8, 20];
  4. (Sum[Sum[If[Mod[n, k] == 0, (MoebiusMu[k]*k)^(1/2)/n^s, 0], {k, 1,
  5. nn}], {n, 1, nn}])^2
  6. (A = Table[
  7. Table[If[Mod[n, k] == 0, (MoebiusMu[k]*k)^(1/2)/n^s, 0], {k, 1,
  8. nn}], {n, 1, nn}]);
  9. N[Total[Total[A.Transpose[A]]], 20]
  10. Zeta[s]*Zeta[s]/Zeta[2*s - 1]
  11. (*end*)
  12.  
  13. Out[446]= 1.00814049998960894694 + 0.01168365002895327553 I
  14.  
  15. Out[448]= 1.0081404999896092120
  16.  
  17. Out[449]= 1.008140499989609234
  18.  
  19.  
  20. (*start*)
  21. nn = 200;
  22. s = N[8, 20];
  23. (Sum[Sum[If[Mod[n, k] == 0, (MoebiusMu[k])^(1/2)/n^s, 0], {k, 1,
  24. nn}], {n, 1, nn}])^2
  25. (A = Table[
  26. Table[If[Mod[n, k] == 0, (MoebiusMu[k])^(1/2)/n^s, 0], {k, 1,
  27. nn}], {n, 1, nn}]);
  28. N[Total[Total[A.Transpose[A]]], 20]
  29. Zeta[s]*Zeta[s]/Zeta[2*s]
  30. (*end*)
  31.  
  32. Out[656]= 1.00815593032900188265 + 0.008189189928767649031 I
  33.  
  34. Out[658]= 1.0081559303290019133
  35.  
  36. Out[659]= 1.008155930329001935
Advertisement
Add Comment
Please, Sign In to add comment