import numpy as np
import matplotlib.pyplot as plt

A1=10E13
VO=10
QE=10E-2 #courbe plus cohérente que avec 10E2
Ea=99.7*10E2
R=8.31
CAE=5*10E2
TE1=290
TE2=300
TE3=310
Drh=-20000
Cp=2200
u=850



def X(T):
    return (A1*VO*(np.exp(-Ea/(R*T)))/((A1*VO*(np.exp(-Ea/(R*T))))+QE))

def Z1(T):
    return((T-TE1)/-CAE/Drh*Cp*u)

def Z2(T):
    return((T-TE2)/-CAE/Drh*Cp*u)

def Z3(T):
    return((T-TE3)/-CAE/Drh*Cp*u)

T=np.linspace(250,400,200)
y=X(T)
z1=Z1(T)
z2=Z2(T)
z3=Z3(T)

plt.title("Evolution de X(T)")
plt.xlabel("T")
plt.ylabel("X(T)")
plt.xlim(250,400)
plt.ylim(-0.1,1.1)
plt.plot(T,y,'b',label="cinétique")
plt.plot(T,z1,'k',label="thermodynamique TE=290K")
plt.plot(T,z2,'k--',label="thermodynamique TE=300K")
plt.plot(T,z3,'k-.',label="thermodynamique TE=310K")
plt.grid()
plt.legend()
plt.show()