from matplotlib import pyplot as plt
from math import *
a=input("nombre de point")
n=int(a)
k1=1.3
k2=2*10**4
k5=34
k6=3*10**3
k7=2*10**-2
A=10**-2
f=1.1
X=[0 for i in range (10000*n)]
Y=[0 for i in range (10000*n-1)]
Z=[0 for i in range (10000*n-1)]
Y.insert(0,10**-3)
Z.insert(0,10**-3)

for i in range (10000*n-1):
    X[i+1]=X[i]+(1/n)*(k1*A*Y[i]-k2*X[i]*Y[i]+k5*A*X[i]-2*k6*((X[i])**2))
    Y[i+1]=Y[i]+(1/n)*(-k1*A*Y[i]-k2*X[i]*Y[i]+f*k7*Z[i])
    Z[i+1]=Z[i]+(1/n)*(k5*A*X[i]-k7*Z[i])

T=[0 for i in range(10000*n)]
for i in range(10000*n-1):
    T[i+1]=T[i]+1/n

plt.plot(T,X,label="X")
plt.plot(T,Y,label="Y")
plt.plot(T,Z,label="Z")
plt.legend()
plt.show()