Python实现计算圆周率π的值到任意位的方法示例
作者:开心果汁 时间:2021-09-08 16:47:09
本文实例讲述了Python实现计算圆周率π的值到任意位的方法。分享给大家供大家参考,具体如下:
一、需求分析
输入想要计算到小数点后的位数,计算圆周率π的值。
二、算法:马青公式
π/4=4arctan1/5-arctan1/239
这个公式由英国天文学教授约翰·马青于1706年发现。他利用这个公式计算到了100位的圆周率。马青公式每计算一项可以得到1.4位的十进制精度。因为它的计算过程中被乘数和被除数都不大于长整数,所以可以很容易地在计算机上编程实现。
三、python语言编写出求圆周率到任意位的程序如下:
# -*- coding: utf-8 -*-
from __future__ import division
####################导入时间模块
import time
###############计算当前时间
time1=time.time()
################算法根据马青公式计算圆周率####################
number = int(raw_input('请输入想要计算到小数点后的位数n:'))
# 多计算10位,防止尾数取舍的影响
number1 = number+10
# 算到小数点后number1位
b = 10**number1
# 求含4/5的首项
x1 = b*4//5
# 求含1/239的首项
x2 = b// -239
# 求第一大项
he = x1+x2
#设置下面循环的终点,即共计算n项
number *= 2
#循环初值=3,末值2n,步长=2
for i in xrange(3,number,2):
# 求每个含1/5的项及符号
x1 //= -25
# 求每个含1/239的项及符号
x2 //= -57121
# 求两项之和
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'总共耗时:' + str(time2 - time1) + 's'
运行结果:
请输入想要计算到小数点后的位数n:20
3.14159265358979323846
总共耗时:9.77699995041s
请输入想要计算到小数点后的位数n:50
3.14159265358979323846264338327950288419716939937510
总共耗时:2.30099987984s
运行截图如下:
PS:这里再为大家推荐几款计算工具供大家参考使用:
在线一元函数(方程)求解计算工具:
http://tools.jb51.net/jisuanqi/equ_jisuanqi
科学计算器在线使用_高级计算器在线计算:
http://tools.jb51.net/jisuanqi/jsqkexue
在线计算器_标准计算器:
http://tools.jb51.net/jisuanqi/jsq
希望本文所述对大家Python程序设计有所帮助。
来源:https://blog.csdn.net/u013421629/article/details/72640062
标签:Python,计算,圆周率
![](/images/zang.png)
![](/images/jiucuo.png)
猜你喜欢
解决Python传递中文参数的问题
2021-04-10 09:00:07
Ext.js源码研读总结
2009-03-04 12:06:00
Python如何使用神经网络进行简单文本分类
2022-05-10 00:45:26
![](https://img.aspxhome.com/file/2023/5/67755_0s.jpg)
Django使用HTTP协议向服务器传参方式小结
2023-03-08 11:31:20
![](https://img.aspxhome.com/file/2023/6/71336_0s.png)
golang中的空slice案例
2023-09-02 12:26:36
![](https://img.aspxhome.com/file/2023/3/101043_0s.png)
MySQL的存储过程写法和Cursor的使用
2008-12-03 15:55:00
php中Array2xml类实现数组转化成XML实例
2023-07-14 21:48:13
JavaScript中Infinity(无穷数)的使用和注意事项
2023-08-21 00:59:13
Silverlight与Flash的技术比较
2009-04-19 18:45:00
Oracle 9i轻松取得建表和索引的DDL语句
2010-07-16 13:09:00
如何获得上一个月份是几月?
2009-11-23 20:38:00
简单的asp采集代码教程
2011-04-18 10:39:00
Python测试开源工具splinter安装与使用教程
2022-07-22 01:56:36
![](https://img.aspxhome.com/file/2023/3/65083_0s.png)
python调用pyaudio使用麦克风录制wav声音文件的教程
2023-12-16 13:32:15
解决Django中修改js css文件但浏览器无法及时与之改变的问题
2023-10-20 05:54:08
![](https://img.aspxhome.com/file/2023/0/63530_0s.jpg)
[译]Javascript风格要素(二)
2008-02-29 12:51:00
如何解决ASP脚本运行超时的错误
2007-10-14 11:38:00
客户体验
2008-12-11 13:58:00
python 中yaml文件用法大全
2021-11-30 08:02:29
![](https://img.aspxhome.com/file/2023/4/65964_0s.png)
sql查询点滴记录
2012-01-05 19:07:49