python 蒙特·卡罗法计算圆周率

【题目描述】

蒙特·卡罗方法是一种通过概率来得到问题近似解的方法，在很多领域都有重要的应用，其中就包括圆周率近似值的计问题。假设有一块边长为2的正方形木板，上面画一个单位圆，然后随意往木板上扔飞镖，落点坐标(x,y)必然在木板上（更多的时候是落在单位圆内），如果扔的次数足够多，那么落在单位圆内的次数除以总次数再乘以4，这个数字会无限逼近圆周率的值。这就是蒙特·卡罗发明的用于计算圆周率近似值的方法。编写程序，模拟蒙特·卡罗计算圆周率近似值的方法，输入掷飞镖次数，然后输出圆周率近似值。

【源代码程序】

import random

 

 

def monte_carlo_pi(throws):

    inside_circle = 0

    # 模拟掷飞镖

    for _ in range(throws):

        # 在正方形内随机生成一个点的坐标

        x = random.uniform(-1, 1)

        y = random.uniform(-1, 1)

        # 如果点在单位圆内，则计数加1

        if x ** 2 + y ** 2 <= 1:

            inside_circle += 1

    # 使用蒙特·卡罗方法计算圆周率的近似值

    pi_approx = (inside_circle / throws) * 4

    return pi_approx

 

 

# 输入掷飞镖次数

throws = int(input())

 

# 计算圆周率近似值

approx_pi = monte_carlo_pi(throws)

 

# 输出结果

print(approx_pi)

【运行测试】