test.py 文件源码

def limite_central():
    N=5000
    epsilon = 5e-2
    x1 = zeros(N)
    mu = 0.5
    for i in range(N):
        x1[i] = sum(rand(i+1)) / (i+1)

    plt.hist(x1, bins=100, range=(mu - epsilon, mu + epsilon), normed=True)
    formatter = FuncFormatter(to_percent)
    plt.gca().yaxis.set_major_formatter(formatter)

#limite_central()