X = np.random.normal(size=N)
counts, bins = np.histogram(X, bins=50, density=True)
bins = bins[:-1] + (bins[1] - bins[0]) / 2
print(np.trapz(counts, bins))
Как построить функцию массы вероятности в python