import matplotlib.pyplot  as plt
import numpy as np
from scipy.integrate import odeint

# Déclaration des constantes
Ne=400
tau=0.01
zeta=-0.1 # essayer d'autres valeurs de zeta, en changeant la valeur et le signe ; zeta est le coefficient d'amortissement
f0=1000
w0=2*np.pi*f0
Te=5/f0/Ne

def f(y,t_ech):
 # à compléter

yini=np.array([0.1,0.])  # 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[:,0])
plt.show()