Estimation d'une aire par la méthode de Monte-Carlo

[_funambule_]

Si quelqu'un peut m'aider à mieux comprendre l'exercice, qui de base fait partie d'un dm sur les probabilités... 

 

[Black Jack]

Bonjour,

 

from random import*
def point():
    x=uniform(0,5)
    y=uniform(0,2)
    if y<=0.072*x**2+0.2:
      p=1
    else:
      p=0
    return p
    
def echantillon(n):
    c=0
    for i in range(n):
        c = c+point()
    return c   
      
somme = 0
somme = echantillon(500000)
print(somme/500000)

En lançant quelques fois le programme, la valeur sortie est :

0.399216 ; 0.400404 ; 0.399668 ; 0.400524 ...

L'aire A vaut environ 0,4 * (5 * 2) = 4 unités d'aire.

************

Mais ce n'est utile que si tu fais l'effort de comprendre ce qui a été fait et pourquoi.