from math import*
from random import*
from matplotlib.pyplot import*

# Calcul d'une approximation de Pi par la méthode de Monte-Carlo

n=int(input("Combien d'essais dans l'échantillon?"))
r=0

for i in range(n):
    x=random() #Tirage uniforme d'un point du carré unité 
    y=random()
    if x**2+y**2<1: #Vérification de sa position par rapport au cercle de rayon 1
        plot(x,y,'ro') # à l'intérieur
        r=r+1 # on le comptabilise
    else:
        plot(x,y,'bo') #à l'extérieur


p=(r/n) # Calcul de l'approximation de l'aire obtenue

axis([0,1,0,1]) # affichage des résultats de l'expérience
text(0.2,0.2, 'Aire='+str(float(p)),size=30)
show()