检验功效与β错误

def solve():
    from scipy.stats import norm
    import numpy as np
    import matplotlib.pyplot as plt
    mu0, mu1, sigma, n, alpha = ____, ____, ____, ____, ____
    se = sigma/np.sqrt(____)
    z = norm.ppf(____-alpha/____)
    power = round(____-norm.cdf((____)/se)+norm.cdf((____)/se), ____)
    print(____)
    x = np.linspace(____)
    plt.plot(x, norm.pdf(____), '____', label=f'____')
    plt.plot(x, norm.pdf(____), '____', label=f'____')
    plt.axvline(____)
    plt.axvline(____)
    plt.fill_between(x, norm.pdf(____), where=(____)|(____), alpha=____, color='____', label=f'____')
    plt.xlabel(____); plt.ylabel(____); plt.title(____); plt.legend()
    plt.tight_layout(); plt.show()