Advertisement
Guest User

Untitled

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