正态分布检验

def solve():
    from sklearn.datasets import load_iris
    import matplotlib.pyplot as plt
    import numpy as np
    from scipy import stats
    iris = load_iris()
    sepal_length = iris.data[:, ____]
    stat, p_value = stats.normaltest(____)
    print(____)
    print(____)
    print(____)
    plt.figure(figsize=(____, ____))
    plt.subplot(____, ____, ____)
    plt.hist(____)
    x = np.linspace(sepal_length.min(), sepal_length.max(), ____)
    plt.plot(x, stats.norm.pdf(x, np.mean(____), np.std(____)), '____', linewidth=____)
    plt.title(____)
    plt.subplot(____, ____, ____)
    stats.probplot(____)
    plt.title(____)
    plt.tight_layout()
    plt.show()