MatsGranvik

Generalized Harmonic numbers perturbed by an epsilon divided by 2 to the power s

Aug 23rd, 2023 (edited)
123
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.52 KB | None | 0 0
  1. (* start *)
  2. (* Mathematica 8.0.1 *)
  3. n = 20;(*set n=30 for better precision*)
  4. epsilon = 1/10^10;
  5. s = 1 + 14*
  6. I;(*change the number s=1+14*I to get another zero closest to it*)
  7. \
  8.  
  9. s + 1/n +
  10. 1/(1 - Sum[(-1)^(k - 1)*
  11. Binomial[n - 1,
  12. k - 1]/(HarmonicNumber[10^10000, s + k/n] +
  13. epsilon/2^(s + k/n)), {k, 1, n}]/
  14. Sum[(-1)^(k - 1)*
  15. Binomial[n - 1,
  16. k - 1]/(HarmonicNumber[10^10000, s + k/n + 1/n] +
  17. epsilon/2^(s + k/n + 1/n)), {k, 1, n}]);
  18. N[%, 15]
  19. (* end *)
Advertisement
Add Comment
Please, Sign In to add comment