import matplotlib.pyplot as plt
import numpy as np
import numpy.random as rd
from scipy.integrate import odeint
from math import *

fr=         # fréquence de résonance
C= 1e-6
u_fr=   # incertitude sur la fréquence de résonance à estimer à l'aide de l'oscilloscope à completer
 # simulation de la variabilité de L
N= 10000
fr_MC=rd.normal(fr,u_fr,N) #on créer des valeurs de fr réparties selon une distribution normale d'écart type u_fr autour de fr trouvée expérimentalement
# plt.hist(fr_MC) à décommenter pour afficher l'histogramme
# plt.show()     à décommenter pour afficher l'histogramme

L=    # à completer avec la formule qui donne l'expression de L en fonction de C et fr_MC
# Calcul de L et de son incertitude
L1=L.mean()
u_L1=L.std()
print("L'inductance vaut: "+str(L1)+" H +-" +str(u_L1)+"H")