>>> seuil(u,1,2000,250)
(69.7698357481176, 10)