Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- f[y_] := InputForm[Zeta[Rational[1, 2] + Complex[0, 2] Pi (Rational[-11, 8]
- + y)/ProductLog[E^(-1) (Rational[-11, 8] + y)]]]
- v[y_] := InputForm[(Im[Derivative[1][f][y]]^2 + Re[Derivative[1][f][y]]^2)^
- Rational[1, 2]]
- [Kappa][y_] := InputForm[Abs[Im[Derivative[2][f][y]] Re[Derivative[1][f][y]] -
- Im[Derivative[1][f][y]] Re[Derivative[2][f][y]]] (Im[Derivative[1][f][y]]^2 + Re[
- Derivative[1][f][y]]^2)^Rational[-3, 2]]
- d[y_] := InputForm[(Im[f[y]]^2 + Re[f[y]]^2)^Rational[1, 2]]
- r = 2; rr = 15;
- zeros = Rest@ Table[y /. FindRoot[
- 2 Pi (-11/8 + y)/LambertW[(-11/8 + y)/E] == Im[ZetaZero[n]], {y,n}], {n, r, rr}];
- Plot[{v[y], 2 [Pi]/[Kappa][y], [Pi] d[y]}, {y, r, rr},
- PlotRange -> All, PlotTheme -> "Garnet",
- Filling -> {1 -> {2}, 2 -> {3}}, Axes -> False, Frame -> True,
- ImageSize -> 400, GridLines -> {zeros, {}}, Epilog -> {Red, PointSize[Large],
- Point[Transpose@{zeros, Table[v[y], {y, zeros}]}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement