import numpy as np

rng = np.random.default_rng()
x = rng.normal(size = 100)

m = np.mean(x)
np.var(x)
np.mean(x > 2)

0.02

np.mean(x), np.mean((x - m) ** 2), np.mean(x > 2)

(0.041916703859957584, 1.077890480291087, 0.02)