def suite_frac_pt_initial(p,q,p0,q0):
    '''comme suite_frac mais en partant d'un autre point initial p0/q0 '''
    u = Fraction(p0,q0)
    print('Valeur initiale : '+str(float(u)))
    for i in range(100):  # i varie de 0 à 99
        u = (q+1) * u - p
        if (i+1) % 10 == 0 :  # l'indice mathématique est décalé de 1
            approx = float(u)
            print('Le terme de rang '+str(i+1)+' vaut : '+str(approx))
            print('  écart avec la valeur théorique : ',approx-p/q)