Le point est-il dans le quart de disque ?

def dans_le_disque(x, y):
    return x*x + y*y <= 1

print(dans_le_disque(0.5, 0.5))   # True  : 0.5 < 1
print(dans_le_disque(0.9, 0.9))   # False : 1.62 > 1
print(dans_le_disque(0, 1))       # True  : 1 <= 1, le bord compte

On teste $x^2+y^2\le 1$ plutôt que $\sqrt{x^2+y^2}\le 1$ : c'est équivalent et cela évite une racine carrée à chaque point — sur un million de tirages, l'économie compte.