def cvgce(Niter,Npts,eps):
    '''divise l'intervalle en Npts régulièrement espacés, et en chaque point du maillage
    de I^2 ainsi obtenu, calcule lim de u_n avec Niter itérations ; 
    le point du maillage est rangé dans une liste qui est fonction de la limite obtenue,
    identifiée à eps près à 1, 2, 2000 ou autre'''
    lx1, ly1, lx2, ly2, lx2000, ly2000, lxa, lya = [], [], [], [], [], [], [], []
    for i in range(Npts) :
        pi, qi = Npts - 1 + 8*i, 2*(Npts - 1)
        for j in range(Npts) :
            rj, sj = Npts - 1 + 8*j, 2*(Npts - 1)
            zij = sfrac(pi,qi,rj,sj,Niter)
            if abs(zij - 1) < eps :
                lx1.append(pi/qi)
                ly1.append(rj/sj)
            elif abs(zij - 2) < eps :
                lx2.append(pi/qi)
                ly2.append(rj/sj)
            elif abs(zij - 2000) < eps :
                lx2000.append(pi/qi)
                ly2000.append(rj/sj)
            else :
                lxa.append(pi/qi)
                lya.append(rj/sj)
    return(lx1, ly1, lx2, ly2, lx2000, ly2000, lxa, lya)