import numpy as np

n = 100000
r = 1  #radius of the circle
inside = 0

for k in range (0, n):
    # Coordinates in first quadrant:
    [x_coor, y_coor] = np.random.random(2)
    dist = np.sqrt(x_coor**2 + y_coor**2)
    
    if (dist < r):
        inside = inside + 1

print('Pi is esitmated to be:', f'{ 4*inside/n :.4f}')
# Inside is mult. by 4 as we only consider first quadrant