import numpy as np
import matplotlib.pyplot as plt

# Positions mesurées de l'objet xm et de l'image ym
xm=
ym=
# Calcul de la vergence
def vergence(x,y):
    return ()
ve=                  # valeur de la vergence déduite de l'expérience

# Précisions sur les mesures (Delta : demi-largeur de l'intervalle des valeurs mesurées raisonnablement acceptables)
Deltaxm=
Deltaym=


# Nombre de simulations à effectuer 
N = 

# Calculs de la vergence issus de tirages aléatoires avec une distribution de probabilité uniforme
V = []
for i in range(N):
    X = 
    Y = 
    V.append(vergence(X,Y))

# Représentation des valeurs obtenues dans un histogramme
plt.hist(V,bins='rice')
plt.title('Résultat du tirage aléatoire des vergences')
plt.xlabel("v (m^-1)")
plt.show()

# Calcul et affichage moyenne et écart type
v_moy = np.mean(V)
v_ecartype = np.std(V,ddof=1)
print("incertitude-type sur la vergence u(v) =", v_ecartype,"m^-1 (rappel : ne conserver que 2 CS)")
print("vergence expérimentale v =",ve,"m^-1 (rappel : conserver un nombre de CS cohérent avec la précision de u(v))") 
print("vergence moyenne issue des tirages aléatoires v =",v_moy,"m^-1 (rappel : conserver un nombre de CS cohérent avec la précision de u(v))")