我试图绘制一组数据点的概率密度函数，这些数据点是从正态分布中抽样得到的，如下所示：

mu, sigma = 0, 1
N = 100000
s = np.random.normal(mu, sigma, N)
def normal_dist(x , mean , sd):
    prob_density = (np.pi*sd) * np.exp(-0.5*((x-mean)/sd)**2)
    return prob_density

pdf = normal_dist(s, mu, sigma)
plt.plot(s,pdf , color = 'red')
plt.xlabel('Data points')
plt.ylabel('Probability Density')
plt.show()

但这是我得到的图片。我有兴趣得到一行代表pdf的行。我哪里做错了？