這篇文章給大家介紹python實(shí)現(xiàn)計(jì)算π的值到任意位的方法,內(nèi)容非常詳細(xì),感興趣的小伙伴們可以參考借鑒,希望對大家能有所幫助。
具體代碼如下:
# -*- 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('請輸入想要計(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,步長=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)容可以對大家有一定的幫助,可以學(xué)到更多知識(shí)。如果覺得文章不錯(cuò),可以把它分享出去讓更多的人看到。
分享名稱:python實(shí)現(xiàn)計(jì)算π的值到任意位的方法-創(chuàng)新互聯(lián)
網(wǎng)址分享:http://www.rwnh.cn/article6/djpjog.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供網(wǎng)站改版、商城網(wǎng)站、網(wǎng)站收錄、網(wǎng)站制作、品牌網(wǎng)站設(shè)計(jì)、關(guān)鍵詞優(yōu)化
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場,如需處理請聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來源: 創(chuàng)新互聯(lián)
猜你還喜歡下面的內(nèi)容
移動(dòng)網(wǎng)站建設(shè)知識(shí)