猫史档案馆

from decimal import Decimal, getcontext
import math
getcontext().prec = 500
pi=Decimal (4)

print('莱布尼茨级数:')
i=Decimal(3)
f=Decimal(4)
t=Decimal(2)
while(i<10000000):
    pi-=f/i
    i+=t
    pi+=f/i
    i+=t
    if(i%1000000==3):
        print(i//4,'pi=',pi)

print('马青公式:')
pi=f*math.acos(0.2)-math.acos(1/239)
print('pi=',pi)

print('Python中的数学库提供的圆周率：',math.pi)