def Seuil(M): x=0 while 0.0025*x**2-0.2*x+2>M: x=x+0.01 return x