﻿from math import *

def D(n):
    y=4+40.2/(0.005*exp(n)+1)
    return y

def Seuil(S):
    n=0
    while D(n)>S:
        n=n+1
    return n

