import numpy as np

N = 100000                              # nombre de tirages
D = np.random.normal(80.0, 0.5, N)      # valeur mesurée, incertitude, nb de tirages
d = np.random.normal(40.0, 0.5, N)

f = (D**2 - d**2) / (4 * D)             # la formule du TP, appliquée aux N tirages

print("f' =", f.mean(), "cm")
print("u(f') =", f.std(ddof=1), "cm")