求圆周率Π

Posted 2022-04-21 chenziqing

蒙特卡洛方法求Π

import random
count = 0
for i in range(1000000):
    x, y = random.random(), random.random()
    dist = pow(x ** 2 + y ** 2, 0.5)
    if dist < 1:
        count += 1
print((count / 1000000)* 4)

pi = 0
k = 0
while True:

    pi +=  (1/(16**k))*           (4/(8*k+1) - 2/(8*k+4) - 1/(8*k+5) - 1/(8*k+6))

    print(pi)
    k += 1

作　　者：豆瓣酱瓣豆
出　　处：https://www.cnblogs.com/chenziqing/
查看其它博文请点击：https://www.cnblogs.com/chenziqing/
声援博主：如果您觉得文章对您有帮助，可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力！

微信: