Advertisement
MatsGranvik

Characteristic Polynomial alternating series Riemann

Aug 1st, 2021
155
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.93 KB | None | 0 0
  1. (*start*)
  2. number = 10;
  3. log1 = 0;
  4. log2 = Round[Log[2]*number]/number;
  5. log3 = Round[Log[3]*number]/number;
  6. log4 = Round[Log[4]*number]/number;
  7. Clear[s];
  8. x = s /. FindRoot[
  9. 1 - 1/(E^(log2))^s + 1/(E^(log3))^s - 1/(E^(log4))^s == 0, {s,
  10. 1 - I*20}, WorkingPrecision -> 100]
  11. 1 - 1/(E^(Log[2]))^s + 1/(E^(Log[3]))^s - 1/(E^(Log[4]))^s
  12. s = x;
  13. "(1) alternating series"
  14. 1 - 1/(E^(Log[2]))^s + 1/(E^(Log[3]))^s - 1/(E^(Log[4]))^s
  15. 1 - 1/(E^(log2))^s + 1/(E^(log3))^s - 1/(E^(log4))^s
  16. Clear[x]
  17. 1 - x^Round[Log[2]*number] + x^Round[Log[3]*number] -
  18. x^Round[Log[4]*number]
  19. "(2) alternating polynomial"
  20. X = Exp[-s/number]
  21. 1 - X^Round[Log[2]*number] + X^Round[Log[3]*number] -
  22. X^Round[Log[4]*number]
  23. (*end*)
  24. h = {log1, log2, log3, log4}*number
  25. "(3) Characteristic polynomial"
  26. Expand[CharacteristicPolynomial[
  27. Table[Table[(-I)^(n + k - 0)*x^((h[[n]] + 1)/2)*
  28. x^((h[[k]] + 1)/2), {k, 1, 4}], {n, 1, 4}], x]/x^4 - 1]
  29. (*end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement