def seuil(f,M,u0):
    n=0
    u=u0
    while u<=M:
        n=n+1
        u=f(u)
    return n