python求pi的方法


Posted in Python onOctober 08, 2014

本文实例讲述了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程序设计有所帮助。

Python 相关文章推荐
Python 时间处理datetime实例
Sep 06 Python
Python实现抓取网页生成Excel文件的方法示例
Aug 05 Python
Python实现批量压缩图片
Jan 25 Python
利用scrapy将爬到的数据保存到mysql(防止重复)
Mar 31 Python
python制作mysql数据迁移脚本
Jan 01 Python
python 通过麦克风录音 生成wav文件的方法
Jan 09 Python
django框架防止XSS注入的方法分析
Jun 21 Python
python之拟合的实现
Jul 19 Python
Python+OpenCV实现旋转文本校正方式
Jan 09 Python
解决paramiko执行命令超时的问题
Apr 16 Python
python 通过exifread读取照片信息
Dec 24 Python
python中子类与父类的关系基础知识点
Feb 02 Python
python实现简单的TCP代理服务器
Oct 08 #Python
python操作CouchDB的方法
Oct 08 #Python
python基于queue和threading实现多线程下载实例
Oct 08 #Python
python实现封装得到virustotal扫描结果
Oct 05 #Python
python解析xml文件操作实例
Oct 05 #Python
python写xml文件的操作实例
Oct 05 #Python
python实现上传样本到virustotal并查询扫描信息的方法
Oct 05 #Python
You might like
php运行报错Call to undefined function curl_init()的最新解决方法
2016/11/20 PHP
php中this关键字用法分析
2016/12/07 PHP
几行代码轻松搞定jquery实现flash8类似的连接效果
2007/05/03 Javascript
jQuery EasyUI中对表格进行编辑的实现代码
2010/06/10 Javascript
JS HTML5 音乐天气播放器(Ajax获取天气信息)
2013/05/26 Javascript
浅析基于WEB前端页面的页面内容搜索的实现思路
2014/06/10 Javascript
判断在css加载完毕后执行后续代码示例
2014/09/03 Javascript
介绍JavaScript中Math.abs()方法的使用
2015/06/14 Javascript
浅析jQuery Ajax请求参数和返回数据的处理
2016/02/24 Javascript
使用NodeJs 开发微信公众号(三)微信事件交互实例
2016/03/02 NodeJs
js模拟支付宝密码输入框
2017/04/11 Javascript
JS如何设置元素样式的方法示例
2017/08/28 Javascript
基于openlayers4实现点的扩散效果
2020/08/17 Javascript
Vue项目中添加锁屏功能实现思路
2018/06/29 Javascript
Element Table的row-class-name无效与动态高亮显示选中行背景色
2018/11/30 Javascript
微信小程序实现简单跑马灯效果
2020/05/26 Javascript
Angular6新特性之Angular Material
2018/12/28 Javascript
vue 项目@change多个参数传值多个事件的操作
2021/01/29 Vue.js
在Python中使用matplotlib模块绘制数据图的示例
2015/05/04 Python
浅谈python抛出异常、自定义异常, 传递异常
2016/06/20 Python
Python 模块EasyGui详细介绍
2017/02/19 Python
python解析基于xml格式的日志文件
2017/02/25 Python
python实现求解列表中元素的排列和组合问题
2018/03/15 Python
python实现zabbix发送短信脚本
2018/09/17 Python
JSON文件及Python对JSON文件的读写操作
2018/10/07 Python
python爬虫 基于requests模块的get请求实现详解
2019/08/20 Python
教你使用Canvas处理图片的方法
2017/11/28 HTML / CSS
面向中国市场的在线海淘美妆零售网站:Beauty House美丽屋
2021/03/02 全球购物
鲜果饮品店创业计划书
2014/01/21 职场文书
岗位说明书怎么写
2014/07/30 职场文书
商务经理岗位职责
2014/07/30 职场文书
2014领导班子四风剖析对照检查材料思想汇报
2014/09/20 职场文书
钢铁是怎样炼成的读书笔记
2015/06/29 职场文书
初中语文教学反思范文
2016/03/03 职场文书
聊聊JS ES6中的解构
2021/04/29 Javascript
Springboot如何同时装配两个相同类型数据库
2021/11/17 Java/Android