import sys

def test(n=ZZ(10)):
    bnd = ZZ(10)**n
    alpha = QQ.random_element(bnd,bnd)
    c = continued_fraction(alpha)
    v = convergents(c)
    if alpha != v[-1]:
        print "ERROR: found by contfrac.sage"
        print "alpha = %s"%alpha
        sys.exit(ZZ(1))

k = ZZ(1)
while True:
    test()
    k = k + ZZ(1)
    if k%ZZ(100) == ZZ(0):
        print k