def u(n): return -0.2*n**2+1.2*n+1000 def seuil(u,debut,s): u_n=u(debut) compteur=debut while u_n>s: compteur+=1 u_n=u(compteur) return u_n,compteur