Advertisement
FSZ

NewMathema Mfied

FSZ
Feb 24th, 2020
144
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.83 KB | None | 0 0
  1. Block[{$MaxExtraPrecision = 10000}, epsloc = {}; emptwob = {};
  2. For[v = 2, v < 106,
  3. v++ For[j = 4, j < 5, j++,
  4. For[h = 1, h < j - 1, h++,
  5. If[h == (j - 1)/2, Continue,
  6. For[f = 1, f < v + 1, f++,
  7. g = Floor[(1/
  8. 2*(((2*h*(j^(v) - 1)/(j - 1)) + 1/2)^(1/f) + 1/2))];
  9. d = g + Sqrt[g^2 - 1]; If[d > j, x = (d^f + (1/d^f))/2;
  10. y = h*(j^v - 1)/(j - 1);
  11. If[x == y, kappaa = Log[j, d]; z = 1;
  12. q = Denominator[Convergents[kappaa, z][[z]]];
  13. If[6*10^(15) < q,
  14. epslo = Abs[-Log[j, (2*h)/(j - 1)]*q -
  15. Round[-Log[j, (2*h)/(j - 1)]*q]] -
  16. 10^(15)*Abs[kappaa*q - Round[kappaa*q]];
  17. If[epslo > 0, AppendTo[epsloc, epslo]; AppendTo[emptwob, q],
  18. While[epslo <= 0, z++;
  19. q = Denominator[Convergents[kappaa, z][[z]]];
  20. epslo =
  21. Abs[-Log[j, (2*h)/(j - 1)]*q -
  22. Round[-Log[j, (2*h)/(j - 1)]*q]] -
  23. 10^(15)*Abs[kappaa*q - Round[kappaa*q]]];
  24. AppendTo[epsloc, epslo]; AppendTo[emptwob, q]],
  25. While[q <= 6*10^(15), z++;
  26. q = Denominator[Convergents[kappaa, z][[z]]] ];
  27.  
  28. epslo = Abs[-Log[j, (2*h)/(j - 1)]*q -
  29. Round[-Log[j, (2*h)/(j - 1)]*q]] -
  30. 10^(15)*Abs[kappaa*q - Round[kappaa*q]];
  31. If[epslo > 0, AppendTo[epsloc, epslo]; AppendTo[emptwob, q],
  32. While[epslo <= 0, z++;
  33. q = Denominator[Convergents[kappaa, z][[z]]];
  34. epslo =
  35. Abs[-Log[j, (2*h)/(j - 1)]*q -
  36. Round[-Log[j, (2*h)/(j - 1)]*q]] -
  37. 10^(15)*Abs[kappaa*q - Round[kappaa*q]]];
  38. AppendTo[epsloc, epslo]; AppendTo[emptwob, q]] ],
  39. Continue], Continue]]]]]]; Print[Max[emptwob]];
  40. Print[Min[epsloc]]]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement