Approximation d’une intégrale par la méthode de Monte-Carlo. On s’intéresse à l’aire comprise sous la courbe représentative de la fonction f(x)=x**2 et entre les droites d’équation x=0 et x=1.
from math import * from matplotlib.pyplot import * from random import * def f(x): return x**2 def graph(N): inside=0 X=[i/100 for i in range(101)] Y=[f(x) for x in X] plot(X,Y) for i in range(N): x=random() y=random() if y<f(x): scatter(x,y,color="blue") inside+=1 else: scatter(x,y,color="red") show() return inside/N