MatsGranvik

A divided by B and X divided by Y

Dec 20th, 2020
208
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.73 KB | None | 0 0
  1. Clear[f, A, B, n, s, a, b, x, m];
  2. f[x_] = Zeta[x];
  3. A[n_, s_] =
  4. Sum[(-1)^(k - 1)*Binomial[n - 1, k - 1]/f[s + k/n - 1/n], {k, 1, n}];
  5. B[n_, s_] =
  6. Sum[(-1)^(k - 1)*Binomial[n - 1, k - 1]/f[s + k/n], {k, 1, n}];
  7. X[n_, s_] =
  8. Sum[(-1)^(k - 1)*Binomial[n - 1, k - 1]/f[-s + k/n - 1/n], {k, 1,
  9. n}];
  10. Y[n_, s_] =
  11. Sum[(-1)^(k - 1)*Binomial[n - 1, k - 1]/f[-s + k/n], {k, 1, n}];
  12. n = 120;
  13. s = N[-5 + 2*I, 200];
  14. Block[{$MaxExtraPrecision = 600}, N[A[n, s]/B[n, s], 20]]
  15. Block[{$MaxExtraPrecision = 600}, N[X[n, s]/Y[n, s], 20]]
  16.  
  17.  
  18. Clear[xy, s, aa, bb, ab, AB, XY];
  19. aa = 4;
  20. bb = 2;
  21. ab = AB /. Solve[1/(1 - AB) + s == -aa, AB][[1]][[1]];
  22. xy = XY /. Solve[1/(1 - XY) - s == -bb, XY][[1]][[1]];
  23. s = N[-5 + 2*I, 20];
  24. ab
  25. xy
Advertisement
Add Comment
Please, Sign In to add comment