MatsGranvik

Floor log polynomial Abs value Histogram close to Abs Exp -zeta zeros

Jun 18th, 2022 (edited)
102
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.47 KB | None | 0 0
  1. (*start*)
  2. Clear[rr, j, x]
  3. c = 300;
  4. Sum[(-1)^(n - 1)*x^Floor[Log[n]*c], {n, 1, 100}]
  5. rr = DeleteCases[roots = N[x /. NSolve[Total[%] == 0, x]],
  6. Indeterminate];
  7. Histogram[Abs[rr]];
  8. Graphics[{Pink,
  9. Line[{{Abs[Exp[-(1/2 + I*1)/c]], 0}, {Abs[Exp[-(1/2 + I*1)/c]],
  10. 800}}]}];
  11. Show[%%, %]
  12. ListPlot[Table[{Re[rr[[j]]], Im[rr[[j]]]}, {j, 1, Length[rr]}],
  13. AspectRatio -> 1]
  14. g1 = ListLinePlot[Sort[N[Abs[rr]]]];
  15. rr = N[Exp[-Table[1/2 + I*t, {t, 1, 1400}]/c]];
  16. g2 = ListLinePlot[N[Abs[rr]], PlotStyle -> Pink];
  17. Show[g1, g2]
  18. (*end*)
  19.  
  20.  
  21. (*version 2, logarithms approximated by von Mangoldt function matrix*)
  22. (*start*)
  23. Clear[rr, j, x]
  24. c = 500;
  25.  
  26. nn = c;
  27. a[n_] := Total[Divisors[n]*MoebiusMu[Divisors[n]]];
  28. TableForm[
  29. b = Table[
  30. Sum[If[n == 1, 0, a[GCD[n, k]]/k], {k, 1, nn}], {n, 1, nn}]];
  31. log = Table[
  32. Sum[If[Mod[n, k] == 0, b[[n/k]], 0], {k, 1, n}], {n, 1, nn}];
  33. Show[ListLinePlot[log], Plot[Log[n], {n, 1, nn}, PlotStyle -> Pink]]
  34.  
  35.  
  36. polynomial = Sum[(-1)^(n - 1)*x^Floor[log[[n]]*c], {n, 1, nn}]
  37. rr = DeleteCases[roots = N[x /. NSolve[Total[%] == 0, x]],
  38. Indeterminate];
  39. ListPlot[Table[{Re[rr[[j]]], Im[rr[[j]]]}, {j, 1, Length[rr]}],
  40. AspectRatio -> 1]
  41. Histogram[Abs[rr]];
  42. Graphics[{Pink,
  43. Line[{{Abs[Exp[-(1/2 + I*1)/c]], 0}, {Abs[Exp[-(1/2 + I*1)/c]],
  44. 2000}}]}];
  45. Show[%%, %]
  46. g1 = ListLinePlot[Sort[N[Abs[rr]]]];
  47. rr = N[Exp[-Table[1/2 + I*t, {t, 1, 3000}]/c]];
  48. g2 = ListLinePlot[N[Abs[rr]], PlotStyle -> Pink];
  49. Show[g1, g2]
  50. (*end*)
Add Comment
Please, Sign In to add comment