写一个程序，从一维高斯分布N(μ,σ^2)中产生n个随机数，绘制数据的直方图和概率密度函数，分析说明直方图和概率密度函数的关系。

这是一个很有意思的问题！可以使用Python的NumPy和Matplotlib库来实现。

首先，需要使用NumPy库中的random模块生成n个服从高斯分布的随机数。代码如下：

import numpy as np

n = 1000  # 产生1000个随机数
mu, sigma = 0, 0.1 # 均值μ为0，标准差σ为0.1

s = np.random.normal(mu, sigma, n) # 产生n个随机数，服从N(μ,σ^2)分布

接下来，使用Matplotlib库中的hist函数绘制数据的直方图，使用NumPy库中的linspace函数生成一组横坐标，然后使用Matplotlib库中的plot函数绘制高斯分布的概率密度函数。代码如下：

import matplotlib.pyplot as plt

count, bins, ignored = plt.hist(s, 30, density=True)  # 产生直方图

x = np.linspace(mu - 3 * sigma, mu + 3 * sigma, 100)
pdf = 1 / (sigma * np.sqrt(2 * np.pi)) * np.exp(-(x - mu) ** 2 / (2 * sigma ** 2))  # 产生概率密度函数
plt.plot(x, pdf, linewidth=2, color='r')  # 绘制概率密度函数

plt.show()  # 显示图像

直方图和概率密度函数的关系是：直方图展示了数据的频率分布情况，而概率密度函数则是描述了这一分布的数学模型。同时，从图像可以看出，直方图和概率密度函数是很相似的，因为它们都反映了随机数在高斯分布中的分布情况。