中心极限定理模拟

def solve():
    import numpy as np
    import matplotlib.pyplot as plt
    from scipy import stats
    np.random.seed(____)
    means = [np.mean(np.random.uniform(____,____,____)) for _ in range(____)]
    print(____)
    x = np.linspace(____,____,____)
    plt.hist(____)
    plt.plot(x, stats.norm.pdf(x,____,____/np.sqrt(____)), '____', lw=____)
    plt.title(____); plt.tight_layout(); plt.show()