展开全部
#coding=utf-8
'''
Created on 2014-11-04
@author: Neo
'''
import sys
import math
from decimal import *
def bbp(n):
pi=Decimal(0)
k=0
while k
pi+=(Decimal(1)/(16**k))*((Decimal(4)/(8*k+1))-(Decimal(2)/(8*k+4))-(Decimal(1)/(8*k+5))-(Decimal(1)/(8*k+6)))
k+=1
return pi
def main(argv):
if len(argv) !e69da5e6ba903231313335323631343130323136353331333335343331=2:
exit('Usage: BaileyBorweinPlouffe.py ')
getcontext().prec=(int(sys.argv[1]))
my_pi=bbp(int(sys.argv[2]))
accuracy=100*(Decimal(math.pi)-my_pi)/my_pi
print "Pi is approximately "+str(my_pi)
print "Accuracy with math.pi: "+str(accuracy)
if __name__=="__main__":
main(sys.argv[1:])
result:
d:\workspace\PyDemo>python test.py 10 10
Pi is approximately 3.141592653
Accuracy with math.pi: 1.877369797E-8
d:\workspace\PyDemo>python test.py 25 25
Pi is approximately 3.141592653589793238462644
Accuracy with math.pi: -3.898171852150198570978563E-15
d:\workspace\PyDemo>python test.py 40 40
Pi is approximately 3.141592653589793238462643383279502884195
Accuracy with math.pi: -3.898171832519375446564294863483061824361E-15
d:\workspace\PyDemo>