蒙特卡洛模拟

def solve():
    import numpy as np
    import statsmodels.api as sm
    import matplotlib.pyplot as plt
    
    np.random.seed(____)
    betas = []
    for _ in range(____):
        x = np.random.randn(____)
        y = ____ + ____*x + np.random.randn(____)
        X = sm.add_constant(____)
        b = sm.OLS(____).fit().params[____]
        betas.append(____)
    
    plt.figure(figsize=(____,____))
    plt.hist(____)
    plt.axvline(____)
    plt.axvline(x=np.mean(____), color='____', linestyle='____', linewidth=____, label=f'____')
    plt.title(____)
    plt.xlabel(____)
    plt.ylabel(____)
    plt.legend()
    plt.grid(____)
    plt.tight_layout()
    plt.show()
    
    print(round(np.mean(____), ____))