import numpy as np, matplotlib.pyplot as plt

Ta = 283
U = 10
A = 100
Cp = 1000/3
drxH = -18000
T0 = 293
Fa0 = 6
Ca0 = 3

R = 1.987

A1 = 17.2
A2 = 41.7
E1 = 11600
E2 = 29600

Xa = np.linspace(0,1,1000)

T = ((U*A*Ta)/Fa0 + T0*Cp - Xa*drxH)/((U*A/Fa0) + Cp)
k1 = np.exp(A1 - E1/(R*T))
k2 = np.exp(A2 - E2/(R*T))

r = Ca0*(k1 - Xa*(k1+k2))

print('max r : ',max(r))
print('X : ', Xa[list(r).index(max(r))])
print('T : ', T[list(r).index(max(r))])
plt.plot(Xa, r)
plt.show()

plt.plot(T,r)
plt.show()
