def seuil_var(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(a,a,p,q,20)