1import matplotlib.pyplot as plt
2import numpy as np
3import scipy.stats as stats
4import math
5
6mu = 0
7variance = 1
8sigma = math.sqrt(variance)
9x = np.linspace(mu - 3 * sigma, mu + 3 * sigma, 100)
10plt.plot(x, stats.norm.pdf(x, mu, sigma))
11plt.show()