程式碼如下：

#-*- coding:utf-8 -*-
# Python實現正態分佈
# 繪製正態分佈概率密度函式
import numpy as np
import matplotlib.pyplot as plt
import math
u = 0   # 均值μ
u01 = -2
sig = math.sqrt(0.2)  # 標準差δ

x = np.linspace(u - 3*sig, u + 3*sig, 50)
y_sig = np.exp(-(x - u) ** 2 /(2* sig **2))/(math.sqrt(2*math.pi)*sig)
print x
print "="*20
print y_sig
plt.plot(x, y_sig, "r-", linewidth=2)
plt.grid(True)
plt.show()