"""
Created on Sun Sep  8 15:22:11 2024

@author: starons
"""
"""
Estimation de la valeur de pi par simulation de Monte-Carlo
"""
# On importe les bibliothèques utiles
# -----------------------------------
import numpy as np
from numpy import random

# Simulation de Monte-Carlo
# -------------------------
N = 10000
x=random.uniform(0,1,N)
y=random.uniform(0,1,N)
r=np.sqrt(x**2+y**2)    # calcul du module

# comptage du nombre de points de module inférieur à 1
# ----------------------------------------------------
n=0
for i in range (N):
    if r[i]<1:
        n=n+1
        
# affichage du résultat 
# ---------------------
print("Une valeur approchée de pi est :",4*n/N)