from math import *
TB=373
LB=42000
R=8.314
C=75
Qpoint=-1500
LA=37800
n0=1
TH=357

def Trosee(x,Teb,Lvap):
    return(Teb/(1-(R*Teb*log(x))/Lvap))

def rosee_eau(x):
    return(Troseex,TB,LB)

def TdeXb(x):
    return(rosee_eau(x))

def XbdeT(T):
    return(exp(-(LB/R)*((1/T)-(1/TB))))

def dt1(dT):
    return(C*dT/Qpoint)

def refroidissement(xb):
    T1=TdeXb(xb)
    Ti=T1+20
    Tf=TH-10
    t=0
    T=Ti
    tt=[0]
    TT=[Ti]
    dT=(T1-Ti)/100
    for i in range(100):
        T+=dt
        t+=dt1(dT)
        tt.append(t)
        TT.append(T)