from math import *
def seuil():
    N=0
    U=3
    while ... :
        N=N+1
        U=...
    return ...