Python割圆法求圆周率

Python中可以使用Leibniz公式（也称无穷级数法）来估算圆周率π，最常见的是使用Monte Carlo模拟的割圆法。这种方法通过随机投掷点到单位圆内和圆外，然后统计比例来逼近π/4的值。以下是简单的Python代码示例：

import random

def estimate_pi(n):
    inside_circle = 0
    total_points = 0

    for _ in range(n):
        x = random.uniform(0, 1)
        y = random.uniform(0, 1)
        distance_from_origin = (x**2 + y**2) ** 0.5
        if distance_from_origin <= 1:
            inside_circle += 1
        total_points += 1

    pi_estimate = 4 * inside_circle / total_points
    return pi_estimate

# 调用函数并打印结果
estimated_pi = estimate_pi(1000000)  # 指定更多的点会得到更精确的结果
print(f"估计的圆周率π ≈ {estimated_pi}")