用蒙特卡洛方法求解π

好的，下面是一个用Python实现用蒙特卡洛方法求解π的示例代码：

import random

N = 1000000  # 随机数总数
n = 0  # 在单位圆内的随机数个数

for i in range(N):
    x = random.uniform(-1, 1)
    y = random.uniform(-1, 1)
    if x**2 + y**2 <= 1:
        n += 1

pi = 4 * n / N  # 计算π的估计值
print("π的估计值为:", pi)

这段代码生成了N个均匀分布在正方形区域内的随机数，其中x和y的取值范围都是[-1,1]。然后判断每个点是否落在单位圆内（即判断$x^2+y^2$是否小于等于1），将在圆内的点的数量n累加起来。最后，根据蒙特卡洛方法的公式计算π的估计值，并输出结果。

需要注意的是，随机数的个数N可以根据需要调整，生成的随机数越多，计算结果越精确。