import matplotlib.pyplot as plt
import numpy as np

# Données issues de la droite de pente 1.75
inv_a = np.linspace(0, 50, 100)
L = 1.75 * inv_a

plt.figure(figsize=(7, 5.5))
plt.plot(inv_a, L, color='purple', linewidth=2, label=r'$L = 1{,}75 \cdot \frac{1}{a}$')

# Points remarquables lisibles sur le quadrillage
plt.scatter([0, 20, 40], [0, 35, 70], color='red', zorder=5, label='Points de lecture')

# Titres et étiquettes
plt.title("Modélisation linéaire : L en fonction de 1/a", fontsize=12, fontweight='bold', pad=15)
plt.xlabel(r'$\frac{1}{a} \ (\mathrm{mm}^{-1})$', fontsize=11)
plt.ylabel(r'$L \ (\mathrm{mm})$', fontsize=11)

# Configuration précise du pas des axes (Pas de 10 en Y)
plt.xticks(np.arange(0, 61, 10))
plt.yticks(np.arange(0, 101, 10))

plt.xlim(0, 60)
plt.ylim(0, 100)
plt.grid(True, which='both', linestyle='--', linewidth=0.5, color='gray')
plt.legend(loc='upper left')
plt.show()