import numpy as np              #bibliothèque numpy renommée np

data=[39900, 39950, 39980, 40000, 40030, 40080]  #  !!!!!liste des 6 frequences à remplacer
N=len(data)                                      #nombre d'observations
moy=np.average(data)                             #moyenne des valeurs de la liste data
sx=np.std(data, ddof=1)                          #écart-type de la liste data avec N-1 au dénominateur grâce à ddof = 1
print('moyenne arithmétique = ', moy)
print('écart type= ', sx)
print('incertitude-type u  = ', sx/np.sqrt(N))


print("moyenne = {:.3e} ".format(moy)+"Hz")
print("incertitude-type  u = {:.1e}".format(sx/np.sqrt(N))+"Hz")