Advertisement
FSZ

NewSMath

FSZ
Feb 24th, 2020
161
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.60 KB | None | 0 0
  1. epsloc=[]
  2. emptwob=[]
  3. for v in range(2,3):
  4. for j in range(4,5):
  5. for h in range(1,j-1):
  6. if h==(j-1)/2:
  7. continue
  8. for f in range(1,v+1):
  9. g=floor((1/2*(((2*h*(j^(v)-1)/(j-1))+1/2)^(1/f)+1/2)))
  10. d=g+sqrt(g^2-1)
  11. if d>j:
  12. x=(d^f+(1/d^f))/2
  13. y=h*(j^v-1)/(j-1)
  14. if x==y:
  15. kappaa=log(d,j)
  16. z=0
  17. q=denominator(continued_fraction(kappaa).convergent(z))
  18.  
  19. if 6*10^(15)<q:
  20. epslo=abs(-log((2*h)/(j-1),j)*q-round(-log((2*h)/(j-1),j)*q))-10^(15)*abs(kappaa*q-round(kappaa*q))
  21. if epslo>0:
  22. epsloc.append(epslo)
  23. emptwob.append(q)
  24. else:
  25. while q<=6*10^(15):
  26. z+=1
  27. q=denominator(continued_fraction(kappaa).convergent(z))
  28. if 6*10^(15)<q:
  29. epslo=abs(-log((2*h)/(j-1),j)*q-round(-log((2*h)/(j-1),j)*q))-10^(15)*abs(kappaa*q-round(kappaa*q))
  30. if epslo>0:
  31. epsloc.append(epslo)
  32. emptwob.append(q)
  33. else:
  34. while epslo<=0:
  35. z+=1
  36. q=denominator(continued_fraction(kappaa).convergent(z))
  37. epslo=abs(-log((2*h)/(j-1),j)*q-round(-log((2*h)/(j-1),j)*q))-10^(15)*abs(kappaa*q-round(kappaa*q))
  38. if epslo>0:
  39. epsloc.append(epslo)
  40. emptwob.append(q)
  41.  
  42.  
  43.  
  44. else:
  45. continue
  46.  
  47. else:
  48. continue
  49. print(max(emptwob))
  50. print(N(min(epsloc)))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement