写一个计算概率分布函数的代码

以下是一个计算概率分布函数的Python代码示例：

import numpy as np
import matplotlib.pyplot as plt

# 定义概率密度函数
def pdf(x):
    return np.exp(-x**2/2)/np.sqrt(2*np.pi)

# 定义概率分布函数
def cdf(x):
    return (1 + erf(x/np.sqrt(2))) / 2

# 生成数据
x = np.linspace(-5, 5, 1000)
y_pdf = pdf(x)
y_cdf = cdf(x)

# 绘制图形
fig, ax = plt.subplots(1, 2, figsize=(10, 4))
ax[0].plot(x, y_pdf)
ax[0].set_title('Probability Density Function')
ax[0].set_xlabel('x')
ax[0].set_ylabel('pdf(x)')
ax[1].plot(x, y_cdf)
ax[1].set_title('Cumulative Distribution Function')
ax[1].set_xlabel('x')
ax[1].set_ylabel('cdf(x)')
plt.show()

该代码使用了numpy和matplotlib库来生成数据和绘制图形。函数pdf定义了概率密度函数，函数cdf定义了概率分布函数。在代码的后半部分，使用linspace生成了一组数据x，并用pdf和cdf函数计算了对应的y值。最后，使用matplotlib库绘制了两张图，一张是概率密度函数图，一张是累积分布函数图。