import numpy as np
import matplotlib.pyplot as plt

rng = np.random.default_rng()

N = 101
p = 0.75
x = rng.binomial(1, p, size = N) # Bernoulli(p = 0.75)

np.mean(x == 1)

0.7524752475247525

ndx = np.arange(1, N + 1)
cm = np.cumsum(x == 1) / ndx
plt.plot(ndx, cm)