import matplotlib.pyplot  as plt
import numpy as np
from scipy.integrate import odeint
Ne=400;tau=0.01;Te=5*tau/Ne # Déclaration des constantes

def f(y,t):
    return -y/tau # Définition de la fonction f
yini=2.  # Condition initiale
t_ech=np.array([n*Te for n in range(Ne)])  # tableau de valeurs des instants t

sol=odeint(f,yini,t_ech)
plt.plot(t_ech,sol)
plt.axis([0,5*tau,0,yini])
plt.show()