from math import *
from matplotlib import pyplot as plt
import numpy as np
R=8.314
DrH=101800
DrS=-11.8

def K(T,A,B):
    a=-A+T*B
    a=a/(R*T)
    return(exp(a))

def Y(T,P):
    A=K(T,DrH,DrS)
    Pco=P/(1+A)
    Pco2=P*A/(1+A)
    return(Pco2/Pco)




T=np.linspace(1400,2000,1000)
X=[Y(e,0.5) for e in T]
F=[K(e,DrH,DrS) for e in T]
plt.plot(T,F)
plt.plot(T,X)
plt.show()