Python 作业( 有进度条的圆周率计算 )

 

这里附上圆周率的计算公式:

 

import time

def Countpi(N = 100):
    pi = 0.0
    k = 0
    start = time.perf_counter()
    for i in range(N):
        pi += 1 / pow(16, i) * (
            4/(8*i + 1) - 2/(8*i 
            + 4) - 1/(8*i + 5) - 1/(8*i + 6))
        if i % 2 == 0:
            k += 1
        Progress(i, start, pi, k)
    return pi

def Progress(i, start, pi, k, L = 50):
    a = '*' * k
    b = '.' * (L - k)
    c = (k/L) * 100
    dur = time.perf_counter() - start
    print('
{:^3.0f}%[{}->{}]计算时间为: {:.2f}s'.format(c, a, b, dur), end='')
    time.sleep(0.1)

def main():
    print('圆周率计算开始'.center(50, '-'))
    pi = Countpi()
    print('
'+'圆周率计算结束'.center(50, '-'))
    print('圆周率计算得: {}'.format(pi))

main()

 结果就是这样啦:

　　其实还有一种蒙特卡罗方法可以计算圆周率, 它是通过在这个正方形内部，随机产生10000个点（即10000个坐标对 (x, y)），计算它们与中心点的距离，从而判断是否落在圆的内部。如果这些点均匀分布，那么圆内的点应该占到所有点的 π/4，因此将这个比值乘以4，就是π的值。

下面就直接附上这种方法的源码(不过没有进度条) :

from random import random
from time import perf_counter
DARTS = 1000*1000
hits = 0.0
start = perf_counter()
for i in range(1, DARTS+1):
    x, y = random(), random()
    dist = pow(x ** 2 + y ** 2, 0.5)
    if dist <= 1.0:
        hits = hits + 1
pi = 4 * (hits/DARTS)
print("圆周率值是: {}".format(pi))
print("运行时间是: {:.5f}s".format(perf_counter() - start))