蒙特卡洛求π

def solve():
    import numpy as np
    import matplotlib.pyplot as plt
    np.random.seed(____)
    n = ____
    x = np.random.uniform(____)
    y = np.random.uniform(____)
    inside = x**____ + y**____ <= ____
    pi_estimate = ____ * np.sum(____) / n
    print(____)
    print(____)
    plt.figure(figsize=(____, ____))
    plt.scatter(____)
    plt.scatter(____)
    theta = np.linspace(____, ____*np.pi, ____)
    plt.plot(np.cos(____), np.sin(____), '____', linewidth=____)
    plt.axis(____)
    plt.title(____)
    plt.show()