def seuil(a,eps,k):
    '''attention ça n'a d'intérêt que si a = 1 ou 2 et eps = 1 ou -1'''
    u = Fraction(a+eps*10**(-k))
    p, q = u.numerator, u.denominator
    return sfrac(p,q,p,q,20)

def test_seuil(a,eps,tol=0.1):
    k = 1
    while abs(seuil(a,eps,k) - a) > tol :
        k += 1
    return k