Advertisement
MatsGranvik

Sequence A288640 in OEIS Riemann zeta zeros related

Jun 19th, 2017
210
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.26 KB | None | 0 0
  1. (*Definition:*)
  2. FrancaLeClair[n_] = 2*Pi*Exp[1]*Exp[ProductLog[(n - 11/8)/Exp[1]]];
  3. Monitor[f =
  4. Table[Sign[Im[ZetaZero[n]] - FrancaLeClair[n]], {n, 1, 90}];, n];
  5. (1 + f)/2
  6. (*Comment 1:*)
  7. Monitor[a =
  8. Table[Floor[
  9. Im[ZetaZero[n]]/(2*Pi)*Log[Im[ZetaZero[n]]/(2*Pi*Exp[1])] +
  10. 11/8 - n + 1], {n, 1, 90}];, n];
  11. a
  12. (*Comment 2:*)
  13. Monitor[Table[(1 -
  14. Sign[Im[Zeta[1/2 + I*2*Pi*E*Exp[LambertW[(n - 11/8)/E]]]]])/
  15. 2, {n, 1, 90}], n]
  16. (*Comment 3:*)
  17. Table[Floor[
  18. 2*(RiemannSiegelTheta[Im[ZetaZero[n]]]/Pi -
  19. Floor[RiemannSiegelTheta[Im[ZetaZero[n]]]/Pi])], {n, 1, 90}]
  20. (*Comment 3:*)
  21. (*For n>1:*)
  22. Floor[2*FractionalPart[
  23. N[RiemannSiegelTheta[Im[ZetaZero[Range[90]]]]/Pi, 30]]]
  24. (*Comment 4:*)
  25. Clear[nn, n, k, t, FrancaLeclair, NumberOfZetaZeros];
  26. nn = 89;
  27. FrancaLeclair[n_] = 2*Pi*Exp[1]*Exp[ProductLog[(n - 11/8)/Exp[1]]];
  28. NumberOfZetaZeros[t_] =
  29. RiemannSiegelTheta[t]/Pi + Im[Log[Zeta[1/2 + I*t]] + I*Pi]/Pi;
  30. Monitor[b =
  31. N[Table[-(1 +
  32. 2*Sum[(NumberOfZetaZeros[FrancaLeclair[k + 1]] -
  33. 1) - (NumberOfZetaZeros[FrancaLeclair[k]] - 1) - 1, {k,
  34. 1, n}]), {n, 0, nn}]], n];
  35. (*Notice that sequence b is integer in itself despite the use of the \
  36. Round function below*)
  37. (1 + Round[b])/2
  38. (*End*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement