from fractions import *

def suite_frac(p,q,r,s,N):
    u = Fraction(p,q)
    v = Fraction(r,s)
    for i in range(N): 
        w = v
        v = 2003 - 6002/v + 4000/(u*v)
        u = w
        if (i+1) % 10 == 0 :  
            approx = float(v)
            print('Le terme de rang '+str(i+1)+' vaut : '+str(v))
            print('  valeur approchée : '+str(approx))