## TP 32 Mesure de la célérité du son dans l’air

import numpy as np
import random as rd
def NS(a):
    return('{:.2E}'.format(a))
fr=39950
Dx=0.084
uDx=0.002
Dt=240e-6
uDt=10e-6
def Vs3(x,t):
    return x/t
tirages=[]
for i in range(100000):
    tirages.append(Vs3(rd.normalvariate(Dx,uDx),rd.normalvariate(Dt,uDt)))
Vs,uVs=np.average(tirages),np.std(tirages)
print('Vs3 = ',NS(Vs),'uVs3 = ',NS(uVs))