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=[0.0345999743]
Ptol=[0.0366083871]
j=k=0
for t in range(300,450):
    T.append(t)

while j<150 and Peau[j]<2:
    Peau.append(psat(T[j],373.15,42800))
    j=j+1
Ta=[T[i] for i in range(len(Peau))]
while k<150 and Ptol[k]<2:
    Ptol.append(psat(T[k],383.75,37800))
    k=k+1
Tb=[T[i] for i in range(len(Ptol))]
plt.plot(Ta,Peau,'r',label="Eau")
plt.plot(Tb,Ptol,'g',label="Toluene")
plt.legend()
plt.show()