def norm_pdf(x, mu=0, sigma=1): u = (x - mu)/abs(sigma) y = (1/(math.sqrt(2 * math.pi) * abs(sigma))) * math.exp(-u*u/2) return y