Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- (* Sequence A282896 in the OEIS *)
- (*Comment 1:*)Monitor[
- Flatten[Position[
- Table[Floor[
- Im[ZetaZero[n]]/(2*Pi)*Log[Im[ZetaZero[n]]/(2*Pi*E)] + 11/8] -
- n + 1, {n, 1, 110}], 1]], n]
- (*Comment 2:*)
- Monitor[Flatten[
- Position[Table[
- Sign[Im[Zeta[1/2 + I*2*Pi*E*Exp[LambertW[(n - 11/8)/E]]]]], {n, 1,
- 110}], -1]], n]
- (*From Mats Granvik,Jun 17 2017:(Start)*)
- (*Comment 3:*)Flatten[
- Position[Table[
- Floor[2*(RiemannSiegelTheta[Im[ZetaZero[n]]]/Pi -
- Floor[RiemannSiegelTheta[Im[ZetaZero[n]]]/Pi])], {n, 1, 110}],
- 1]]
- (*Comment 3:*)
- Flatten[Position[
- Floor[2*FractionalPart[
- N[RiemannSiegelTheta[Im[ZetaZero[Range[110]]]]/Pi, 30]]], 1]]
- (*Comment 4:*)
- Clear[nn, n, k, t, FrancaLeclair, NumberOfZetaZeros];
- nn = 110;
- FrancaLeclair[n_] = 2*Pi*Exp[1]*Exp[ProductLog[(n - 11/8)/Exp[1]]];
- NumberOfZetaZeros[t_] =
- RiemannSiegelTheta[t]/Pi + Im[Log[Zeta[1/2 + I*t]] + I*Pi]/Pi;
- Monitor[b =
- N[Table[-(1 +
- 2*Sum[(NumberOfZetaZeros[FrancaLeclair[k + 1]] -
- 1) - (NumberOfZetaZeros[FrancaLeclair[k]] - 1) - 1, {k,
- 1, n}]), {n, 0, nn - 1}]], n];
- (*Notice that sequence b is integer in itself despite the use of the \
- Round function below:*)
- Flatten[Position[Round[b], 1]]
- (*End*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement