def sfrac(p,q,r,s,Niter):
    '''calcule avec le module Fractions les Niter premiers termes
    de la suite et retourne le dernier, en valeur approchée'''
    u = Fraction(p,q)
    v = Fraction(r,s)
    for i in range(Niter): 
        w = v
        v = 2003 - 6002/v + 4000/(u*v)
        u = w
    return(float(v))