# -*- coding: utf-8 -*-
import matplotlib.pyplot as plt
import numpy as np

# Valeurs sur lesquelles on va calculer
x = np.linspace(0, 10 * np.pi, 1000)
# application d'une fonction à chaque point
y1 = np.cos(x) # np.cos s'applique sur des tableaux
y2 = np.exp(-x / 10.0) * y1

# y1 en fonction de x, y2 en fonction de x
plt.plot(x, y1)
plt.plot(x, y2)

# Titres sur les axes
plt.title('Oscillateurs harmoniques libre et amorti')
plt.xlabel('x')
plt.ylabel('y = f(x)')

# Tracé de l'axe y=0 et grille
plt.axhline(y=0, color='k')
plt.grid()

plt.show()