import numpy as np, matplotlib.pyplot as plt
import scipy.stats as sts

z05 = sts.norm.ppf(0.95)
z5 = sts.norm.ppf(0.5)

z05_2 = sts.norm.interval(0.95)

x = np.linspace(-2.5,2.5,200)
N01 = sts.norm.pdf(x)

print(z05, z5, z05_2)
for i in range(10):
    rnds = sts.norm.rvs(size=10)
    plt.scatter(rnds, [i*0.01 for x in rnds])
plt.plot(x,N01)
plt.show()