import matplotlib.pyplot as plt
import matplotlib
from math import *
La=44000
Ta=353
Lb=42800
Tb=373
R=8.314
Y=[0.001*i for i in range(1,1000)]


def f(x,t,m):
    return(x/(1-(R*x*log(m))/t))


P=[]
L=[]

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

U=[]
for e in range (len(Y)):
    U.append(abs(P[e]-L[e]))

j=0
while U[j]>0.1 and j < len(U)-1:
    j=j+1
m=P[j]
d=Y[j]
del P[0:j]
O=list(Y)
del O[0:j]
del L[j:len(L)]
del Y[j:len(Y)]
plt.plot(O,P)
plt.plot(Y,L)
plt.plot([h*0.01 for h in range(100)],[m for k in range (100)])
matplotlib.pyplot.axis([0,1,320,380])
print(m,d)
plt.show()