import matplotlib.pyplot as plt
import matplotlib
from math import *
La=7000
Ta=90
Lb=6000
Tb=77
R=8.314
T=[i for i in range(77,95)]
def f(x,y,t):
    return(exp((-x/R)*((1/y)-(1/t))))

P=[]
L=[]

for u in T:
    P.append((1-f(La,Ta,u))/(f(Lb,Tb,u)-f(La,Ta,u)))
for k in T:
    L.append(f(Lb,Tb,k)*((1-f(La,Ta,k))/(f(Lb,Tb,k)-f(La,Ta,k))))

plt.plot(P,T)
plt.plot(L,T)
matplotlib.pyplot.axis([0,1,70,95])
plt.title("binaire O2/N2 idéal")
plt.show()