中心极限定理验证

def solve():
    import matplotlib.pyplot as plt
    import numpy as np
    from scipy import stats
    np.random.seed(____)
    population_mean = ____
    sample_size = ____
    n_samples = ____
    sample_means = []
    for _ in range(____):
        sample = np.random.exponential(____)
        sample_means.append(np.mean(____))
    sample_means = np.array(____)
    theoretical_std = population_mean / np.sqrt(____)
    print(____)
    print(____)
    print(____)
    fig, (____) = plt.subplots(____, ____, figsize=(____, ____))
    x = np.random.exponential(____)
    ax1.hist(____)
    ax1.set_title(____)
    ax1.set_xlabel(____)
    ax2.hist(____)
    x_norm = np.linspace(sample_means.min(), sample_means.max(), ____)
    ax2.plot(x_norm, stats.norm.pdf(____), '____', linewidth=____, label='____')
    ax2.set_title(____)
    ax2.set_xlabel(____)
    ax2.legend()
    plt.tight_layout()
    plt.show()