Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- epsloc=[]
- emptwob=[]
- for v in range(2,3):
- for j in range(4,5):
- for h in range(1,j-1):
- if h==(j-1)/2:
- continue
- for f in range(1,v+1):
- g=floor((1/2*(((2*h*(j^(v)-1)/(j-1))+1/2)^(1/f)+1/2)))
- d=g+sqrt(g^2-1)
- if d>j:
- x=(d^f+(1/d^f))/2
- y=h*(j^v-1)/(j-1)
- if x==y:
- kappaa=log(d,j)
- z=0
- q=denominator(continued_fraction(kappaa).convergent(z))
- if 6*10^(15)<q:
- 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))
- if epslo>0:
- epsloc.append(epslo)
- emptwob.append(q)
- else:
- while q<=6*10^(15):
- z+=1
- q=denominator(continued_fraction(kappaa).convergent(z))
- if 6*10^(15)<q:
- 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))
- if epslo>0:
- epsloc.append(epslo)
- emptwob.append(q)
- else:
- while epslo<=0:
- z+=1
- q=denominator(continued_fraction(kappaa).convergent(z))
- 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))
- if epslo>0:
- epsloc.append(epslo)
- emptwob.append(q)
- else:
- continue
- else:
- continue
- print(max(emptwob))
- print(N(min(epsloc)))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement