from matplotlib import pyplot as plt
from math import *
import time


k10=20
k20=10
R=8.314
Ea1=70000
Ea2=80000
T0=273

deltaT=input("ecart de temperature par rapport a 293K(ref)")
Y=input("nombre de point")
J=input("echelle en abscisse , de 0 a  ...")
K=int(deltaT)
X=int(Y)
L=float(J)
P=X/L
D=int(P)

k1=k10*exp(-((Ea1)/(R))*((1/(T0+K))-1/(T0)))
k2=k20*exp(-((Ea2)/(R))*((1/(T0+K))-1/(T0)))

A=[0 for i in range(X-1)]
B=[0 for i in range(X)]
C=[0 for i in range(X)]
A.insert(0,1)

for i in range(1,X):
    A[i]=A[i-1]+(1/D)*(-k1)*A[i-1]
    B[i]=B[i-1]+(1/D)*(k1*A[i-1]-k2*B[i-1])
    C[i]=C[i-1]+(1/D)*(k2*B[i-1])

V=[i/D for i in range(X)]
plt.plot(V,A)
plt.plot(V,B)
plt.plot(V,C)
plt.show()