python求pi的方法
作者:shichen2014 发布时间:2023-04-04 06:49:17
标签:python,方法
本文实例讲述了python求pi的方法,是一篇翻译自国外网站的文章,分享给大家供大家参考。
具体实现方法如下:
#_*_ coding=utf-8 *_*
## {{{ http://code.activestate.com/recipes/578130/ (r5)
def pi(places=10):
"""Computes pi to given number of decimal places
参数places表示要返回的pi的小数点后位数
方法:先整体扩大10**8(10的八次方)倍,然后计算完成后再缩小10的八次方倍
"""
# 3 + 3*(1/24) + 3*(1/24)*(9/80) + 3*(1/24)*(9/80)*(25/168)
# The numerators 1, 9, 25, ... are given by (2x + 1) ^ 2
# The denominators 24, 80, 168 are given by (16x^2 -24x + 8)
extra = 8
one = 10 ** (places+extra)
t, c, n, na, d, da = 3*one, 3*one, 1, 0, 0, 24
#这里的n 和d 分别为每一项的分子与分母 ,na 和 da 分别为分子和分分母后一项比前一项增加的数值
#这里的//可不是C++中的注释,而是除的意思
while t > 1:
n, na, d, da = n+na, na+8, d+da, da+32
t = t * n // d
c += t
return c // (10 ** extra)
def picirc(radius, aspect_ratio=5):
"""Display the digit of pi in a circle of given radius
radius:显示的半径
aspect_ratio:调节显示的比率参数
"""
#display_width为各行的显示长度
display_width = int(radius * aspect_ratio + 10)
pi_str = repr(pi(int(2 * radius ** 2 * aspect_ratio)))
pos = 0
#cols为每一行中要显示的数字个数
for i in range(2 * radius):
cols = int(0.5 + aspect_ratio * (radius**2 - (radius-(i+0.5))**2) ** 0.5)
print(pi_str[pos:pos+cols].center(display_width)) #将产生的pi数值生成的文本列表中相应的位数取出来显示在当前行
pos += cols
if __name__ == '__main__':
picirc(16)
## end of http://code.activestate.com/recipes/578130/ }}}
显示的结果为:
31415926535897932384
6264338327950288419716939937510582
0974944592307816406286208998628034825342117
06798214808651328230664709384460955058223172535940
81284811174502841027019385211055596446229489549303819644
288109756659334461284756482337867831652712019091456485669234
6034861045432664821339360726024914127372458700660631558817488152
09209628292540917153643678925903600113305305488204665213841469519415
11609433057270365759591953092186117381932611793105118548074462379962749
5673518857527248912279381830119491298336733624406566430860213949463952247
371907021798609437027705392171762931767523846748184676694051320005681271452
63560827785771342757789609173637178721468440901224953430146549585371050792279
689258923542019956112129021960864034418159813629774771309960518707211349999998
3729780499510597317328160963185950244594553469083026425223082533446850352619311
88171010003137838752886587533208381420617177669147303598253490428755468731159562
86388235378759375195778185778053217122680661300192787661119590921642019893809525
72010654858632788659361533818279682303019520353018529689957736225994138912497217
75283479131515574857242454150695950829533116861727855889075098381754637464939319
2550604009277016711390098488240128583616035637076601047101819429555961989467678
374494482553797747268471040475346462080466842590694912933136770289891521047521
62056966024058038150193511253382430035587640247496473263914199272604269922796
782354781636009341721641219924586315030286182974555706749838505494588586926
9956909272107975093029553211653449872027559602364806654991198818347977535
66369807426542527862551818417574672890977772793800081647060016145249192
17321721477235014144197356854816136115735255213347574184946843852332
3907394143334547762416862518983569485562099219222184272550254256
887671790494601653466804988627232791786085784383827967976681
45410095388378636095068006422512520511739298489608412848
86269456042419652850222106611863067442786220391949
4504712371378696095636437191728746776465757
3962413890865832645995813390478027
59009946576407895126
希望本文所述对大家的Python程序设计有所帮助。
0
投稿
猜你喜欢
- 本文实例讲述了php查询whois信息的方法。分享给大家供大家参考。具体如下:这里使用php通过查询whois信息的网站列表进行查询func
- 这个函数是前几年刚流行小偷程序的时候,偶写来用于小偷程序中截取代码的;可能有些朋友在我以前的代码中看见过了,但没有写用法,现在把调用方法及使
- Python3将数据保存为txt文件的方法,具体内容如下所示:f = open("data/model_Weight.txt&qu
- 在IE比较简单,大家都知道用setHomePage来设置,懒人写法:<a href="#setHomePage"
- 如今我使用 Python 已经很长时间了,但当我回顾之前写的一些代码时,有时候会感到很沮丧。例如,最早使用 Python 时,我写了一个名为
- 使用.net2005自带的SQL-Express连接不上。解决方法:1.网络防火墙阻止数据库连接;2.默认SQL-Express没有启动Sa
- 我发现有的网站利用了SQL SERVER提供的通过EXCHANGE或OUTLOOK收发邮件的扩展存储过程来完成收发和自动处理邮件(这句话太长
- 常用目标检测模型基本都是读取的PASCAL VOC格式的标签,下面代码用于生成VOC格式的代码,根据需要修改即可:from lxml imp
- 前言:前几天上课闲着没事写了一个python敲击木鱼积累功德的小项目,当时纯粹就是写着玩,回顾一下鼠标事件的东西还记不记得,发现这个博客的点
- 阅读上一篇:打造设计你自己的字体 Ⅱ永远都在寻觅字体设计的灵感。夏天过后,我买了一套便宜的书法钢笔,说服自己,它会让我的鸡爬字产生脱胎换骨的
- 使用MySQL,安全问题不能不注意。以下是MySQL提示的23个注意事项:1.如果客户端和服务器端的连接需要跨越并通过不可信任的网络,那么就
- rs.open sql,conn,A,B A: ADOPenforwardonly (=0) 只读,且当前数据记录只能向下移动。 ADOPe
- 腾讯“月捐计划”倡导爱心人士,通过每月小额捐款的形式,长期关注和支持公益项目。并和亿万爱心网友一起,每人每月1份爱,点滴付出,汇成爱海,形成
- 导言在前面的教程里我们学习了DataList提供了一些风格样式的属性.而且我们还学习了如何定义HeadStyle, ItemStyle, A
- 今天记录一下pandas筛选出一个表中满足另一个表中所有条件的数据。例如:list1 结构:名字,ID,颜色,数量,类型。list1 = [
- JAN-1(January) FEB-2(February) MAR-3(March)APR-4(April) MAY-5(Ma
- 元素浮动导致的问题及解决办法大家都应该很熟悉了,举个简单的例子:<style type="text/css">
- 下标所谓下标就是编号,就好比超市中存储柜的编号,通过这个编号就能找到相应的存储空间。Python中字符串,列表,元祖均支持下标索引。例如:#
- 本文通过一个详细的例子,来阐述了在线编辑XML文档数据的方法。由于Netscape对XML的支持比较弱,因此,要实现跨平台的数据交换,数据的
- 实际数据分析中遇到需求,把某个Excel表格按照某一列分为多个sheet,并且要求如果某个key对应的行数较少应该合并到一个sheet中。i