import matplotlib.pyplot as plt
from math import *
R=8.314
def psat(T,TEb,DeltaH):

 return(exp(-(DeltaH/R)*((1/T)-(1/TEb))))

T=[]
Peau=[]
Ptol=[]

for t in range(300,450):
    T.append(t)
for u in range(300,450):
    Peau.append(psat(T[u-300],373.15,42800))
for f in range(300,450):
    Ptol.append(psat(T[f-300],383.75,37800))

plt.plot(T,Peau,'b',label="Eau")
plt.plot(T,Ptol,'g',label="Toluène")
plt.show()
plt.legend()