本文實(shí)例講述了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)=2for 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 resulttime2=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)行截圖如下:


PS:這里再為大家推薦幾款計(jì)算工具供大家參考使用:
在線一元函數(shù)(方程)求解計(jì)算工具:
http://tools.jb51.net/jisuanqi/equ_jisuanqi
科學(xué)計(jì)算器在線使用_高級(jí)計(jì)算器在線計(jì)算:
http://tools.jb51.net/jisuanqi/jsqkexue
在線計(jì)算器_標(biāo)準(zhǔn)計(jì)算器:
http://tools.jb51.net/jisuanqi/jsq
更多關(guān)于Python相關(guān)內(nèi)容感興趣的讀者可查看本站專題:《Python數(shù)學(xué)運(yùn)算技巧總結(jié)》、《Python數(shù)據(jù)結(jié)構(gòu)與算法教程》、《Python函數(shù)使用技巧總結(jié)》、《Python字符串操作技巧匯總》、《Python入門(mén)與進(jìn)階經(jīng)典教程》及《Python文件與目錄操作技巧匯總》
新聞熱點(diǎn)
疑難解答
圖片精選