使用蒙特卡罗方法计算圆周率近似值代码

蒙特卡洛方法是一种基于随机抽样的数值积分方法，通过模拟大量点落在单位正方形和其内切圆内的概率来估算圆周率π。下面是一个简单的Python代码示例：

import random

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

# 例如，我们想要计算100万次样本的圆周率估计
n_samples = 1_000_000
pi_estimate = monte_carlo_pi(n_samples)
print(f"Monte Carlo estimate of π with {n_samples} samples: {pi_estimate}")