小編給大家分享一下如何使用Python實(shí)現(xiàn)計(jì)算圓周率π的值到任意位，希望大家閱讀完這篇文章之后都有所收獲，下面讓我們一起去探討吧！
具體如下：
一、需求分析
輸入想要計(jì)算到小數(shù)點(diǎn)后的位數(shù)，計(jì)算圓周率π的值。
二、算法：馬青公式
π/4=4arctan1/5-arctan1/239
這個(gè)公式由英國(guó)天文學(xué)教授約翰·馬青于1706年發(fā)現(xiàn)。他利用這個(gè)公式計(jì)算到了100位的圓周率。馬青公式每計(jì)算一項(xiàng)可以得到1.4位的十進(jìn)制精度。因?yàn)樗挠?jì)算過(guò)程中被乘數(shù)和被除數(shù)都不大于長(zhǎng)整數(shù)，所以可以很容易地在計(jì)算機(jī)上編程實(shí)現(xiàn)。
三、python語(yǔ)言編寫(xiě)出求圓周率到任意位的程序如下：
重慶分公司
重慶分公司
# -*- 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'運(yùn)行結(jié)果：
請(qǐng)輸入想要計(jì)算到小數(shù)點(diǎn)后的位數(shù)n:20
3.14159265358979323846
總共耗時(shí)：9.77699995041s
請(qǐng)輸入想要計(jì)算到小數(shù)點(diǎn)后的位數(shù)n:50
3.14159265358979323846264338327950288419716939937510
總共耗時(shí)：2.30099987984s
運(yùn)行截圖如下：
看完了這篇文章，相信你對(duì)“如何使用Python實(shí)現(xiàn)計(jì)算圓周率π的值到任意位”有了一定的了解，如果想了解更多相關(guān)知識(shí)，歡迎關(guān)注創(chuàng)新互聯(lián)行業(yè)資訊頻道，感謝各位的閱讀！
