>>> seuil(u,1,5,1000)
(1336, 6)