approximate_pi.py

import numpy as np
import random
from tqdm import tqdm


def test_pi(x, y):
    return int(x ** 2 + y ** 2 <= 1) * 4

def test_it(n):
    return np.mean([
        test_pi(random.uniform(0, 1), random.uniform(0, 1))
        for _ in tqdm(range(n))
    ])

if __name__ == '__main__':
    print(test_it(10000000))