python画随机曲线,用理论曲线绘制直方图：随机实现

我需要编写一个程序来生成Cauchy分布的随机实现

{1美元^

具有空位置和单位比例。在

我还需要做一个介于-5和5个格之间的直方图，对于1000个点的随机实现，和理论曲线确保它们具有相同的单位。在

我计算了Cauchy分布的累积分布函数

：

我写了下面的python代码：

在

在from __future__ import division

import scipy

import random

import matplotlib.pyplot as plt

import numpy as np

import math as m

valuesX = []

for q in range(1000):

R = random.random()

x = m.tan(m.pi*(R-0.5)) #Cumulative Function

valuesX.append(x)

z = np.linspace(-10,10,1000)

y = 1/(m.pi*(1+z**2)) #Theoretical Cauchy

plt.plot(y,z)

plt.hist(valuesX, bins = 50, range = [-5,5], normed=True)

我的结果是：

我不知道这是否可以接受，因为我正在根据概率密度函数绘制离散值(随机实现)。我怎么能比较他们两个呢？因为我需要找到上面曲线的分数差，并确定两条曲线之间的整体均方根偏差，作为随机实现大小的函数。在