蒙特卡洛算法 python

蒙特卡洛方法（Monte Carlo method）是一种以概率统计理论为指导的数值计算方法，它使用随机数（或伪随机数）来解决科学和工程中的许多计算问题。在Python中，可以使用NumPy库来实现蒙特卡洛算法。

下面是一个蒙特卡洛算法的Python代码示例：

import numpy as np

def monte_carlo_pi(n):
    count = 0
    for _ in range(n):
        x = np.random.uniform(0, 1)
        y = np.random.uniform(0, 1)
        if x**2 + y**2 <= 1:
            count += 1
    return 4 * count / n

# 使用蒙特卡洛算法估算圆周率
n = 1000000
pi_estimate = monte_carlo_pi(n)
print("估算的圆周率为:", pi_estimate)

这段代码通过生成n个在[0, 1]范围内的随机点，并计算落在单位圆内的点的比例来估算圆周率。最后输出的结果即为估算的圆周率。