Метод Монте-Карло, число π

Напишите программу, которая при помощи метода Монте-Карло определяет приближённое значение числа π.

import random

n = 10**6
S = 4
k = 0

for _ in range(n):
    x = random.uniform(-1, 1)
    y = random.uniform(-1, 1)
    if x**2 + y**2 <= 1:
        k += 1
        
print((k/n) * S)