#J.COURTIN 01/2023         MONTE CARLO  pour BESSEL

from random import gauss




#formule de Bessel
def focale(D,e):
    return





#Génération des mesuree aléatoire fictive
def mesureAleatoire(D,e,DD,De,N=100):


    #??

    return #??




###### MAIN ##########

D =   #paramètres de distance
e =    #mesurés D et e

DD =    #incertitude mesurée
De =

N=1000     # nombre de mesures simulées


#Simulation :
mes_D, mes_e, ma_focale = mesureAleatoire(D,e,DD,De,N)

print("D = ", mean(mes_D), " +/- ", Sigma(mes_D) )
print("e = ", mean(mes_e), " +/- ", Sigma(mes_e) )


print("focale : f' = ", mean(ma_focale), " +/- ", Sigma(ma_focale) )