def seuil(p): x,y=0,0 while y<2: x=x+p y=x**3+2*x**2 return x