import numpy as np
import matplotlib.pyplot as pltN = 10000
rng = np.random.default_rng()
x = rng.normal(size = N)
xbar = np.mean(x)
np.var(x)1.0116761420756855ndx = np.arange(N) + 1
cm = np.cumsum(x) / ndx
plt.plot(ndx, cm)
cv = np.cumsum((x - xbar) ** 2) / ndx
plt.plot(ndx, cv)
np.mean(x >= 0) # P[X >= 0] = P[X in R^+]0.4934cp = np.cumsum(x >= 0) / ndx
plt.plot(ndx, cp)