python實(shí)現(xiàn)計(jì)算π的值到任意位的方法

這篇文章給大家介紹python實(shí)現(xiàn)計(jì)算π的值到任意位的方法，內(nèi)容非常詳細(xì)，感興趣的小伙伴們可以參考借鑒，希望對(duì)大家能有所幫助。

具體代碼如下：

# -*- coding: utf-8 -*-
from __future__ import division
####################導(dǎo)入時(shí)間模塊
import time
###############計(jì)算當(dāng)前時(shí)間
time1=time.time()
################算法根據(jù)馬青公式計(jì)算圓周率####################
number = int(raw_input('請(qǐng)輸入想要計(jì)算到小數(shù)點(diǎn)后的位數(shù)n:'))
# 多計(jì)算10位，防止尾數(shù)取舍的影響
number1 = number+10
# 算到小數(shù)點(diǎn)后number1位
b = 10**number1
# 求含4/5的首項(xiàng)
x1 = b*4//5
# 求含1/239的首項(xiàng)
x2 = b// -239
# 求第一大項(xiàng)
he = x1+x2
#設(shè)置下面循環(huán)的終點(diǎn)，即共計(jì)算n項(xiàng)
number *= 2
#循環(huán)初值=3，末值2n,步長(zhǎng)=2
for i in xrange(3,number,2):
  # 求每個(gè)含1/5的項(xiàng)及符號(hào)
  x1 //= -25
  # 求每個(gè)含1/239的項(xiàng)及符號(hào)
  x2 //= -57121
  # 求兩項(xiàng)之和
  x = (x1+x2) // i
  # 求總和
  he += x
# 求出π
pai = he*4
#舍掉后十位
pai //= 10**10
############ 輸出圓周率π的值
paistring=str(pai)
result=paistring[0]+str('.')+paistring[1:len(paistring)]
print result
time2=time.time()
print u'總共耗時(shí)：' + str(time2 - time1) + 's'

關(guān)于python實(shí)現(xiàn)計(jì)算π的值到任意位的方法就分享到這里了，希望以上內(nèi)容可以對(duì)大家有一定的幫助，可以學(xué)到更多知識(shí)。如果覺得文章不錯(cuò)，可以把它分享出去讓更多的人看到。