LA=44000
TA=353
LB=42800
TB=373
R=8.314

def dichotomie(xmin,xmax,f,delta):
    while xmax-xmin > delta:
        m=(xmin+xmax)/2
        if f(xmin)*f(m)<0:
            xmax=m
        else:
            xmin=m
    return(m)

def Trosee(x,Teb,Lvap):
    return(Teb/(1-(R*Teb*log(x)/Lvap)))
def rosee_eau(x):
    return(Trosee(x,LB,TB))
def rosee_tol(x):
    return(Trosee(x,TA,LA))

def diff(x):
    return(rosee_eau(x)-rosee_tol(x))
def H(x):
    xH=dicho(diff,0,1,0.01)
    TH=rosee_eau(xH)
    return(xH,TH)
A=[rosee_eau(i*0.01) for i in range(101)]
B=[rosee_tol(i*0.01) for i in range(101)]