Computing PI using Monte Carlo method

最新推荐文章于 2021-10-20 07:01:32 发布

原创最新推荐文章于 2021-10-20 07:01:32 发布 · 547 阅读

0 ·

CC 4.0 BY-SA版权

Python 专栏收录该内容

16 篇文章

订阅专栏

本文介绍了一种利用Python编程实现蒙特卡洛方法来估算圆周率π的算法。通过生成随机点并计算它们与单位圆的关系，我们可以估计π的值。程序包括点的随机生成、距离单位圆中心的计算以及最终π值的计算输出。

import math
import random

class Point():
    def __init__(self, x = 0, y = 0):
        self.x = x
        self.y = y

    def set_rand_range(self, rand_start, rand_end):
        self.rand_start = rand_start
        self.rand_end = rand_end

    def get_random(self):
        r = random.uniform(self.rand_start, self.rand_end)
        return r

    def set_x(self):
        self.x = self.get_random()

    def get_x(self):
        return self.x

    def set_y(self):
        self.y = self.get_random()

    def get_y(self):
        return self.y

def get_randpoint():
    p = Point()
    p.set_rand_range(0, 1)
    p.set_x()
    p.set_y()
    return p.get_x(), p.get_y()

def compute_pi(n):
    x, y = 0, 0
    count = 0
    
    for i in range(n):
        x, y = get_randpoint()
        if math.sqrt(x * x + y * y) <= 1:
            count += 1

    pi = 4.0 * count / n
    print pi
        
        

def main():
    L = [100, 1000, 10000, 100000, 1000000]
    for i in range(len(L)):
        compute_pi(L[i])

if __name__ == '__main__':
    main()
    print math.pi