import numpy as np from scipy import stats mu=10 sig=1 sigma = np.sqrt(sig) x = np.random.normal(mu, sigma, 50) xa = np.average(x) print(xa)