import matplotlib.pyplot as plt
import numpy.random as rd
import numpy as np

def MonteCarlo(alphaM,alphaU,N):
    generateur = rd.default_rng()
    return generateur.normal(alphaM,alphaU,N)

a=MonteCarlo(0.9,0.1,100000)

f=9.81/np.cos(a)+np.tan(a)
plt.hist(f,bins=100)
plt.show()

print(np.mean(f),np.std(f))