Python实现计算圆周率π的值到任意位的方法示例
本文实例讲述了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
运行截图如下:
最新文章
- (转载)android:android.content.res.Resources$NotFoundException: String resource ID #..
- MySQL数据库
- Web程序员常见的5个错误及解决方案
- Android每次运行项目时重新启动一个新的模拟器的解决办法
- C#/JS 获取二维数组组合
- 桂电在linux环境下使用出校器
- android中ADT和SDK的关系(转)
- C# json object互转工具
- 认识DOS
- delphi TOpenDialog
- Linux下find与grep指令的相关用法
- hpu第五届acm比赛
- (Java 多线程系列)java synchronized详解
- KoaHub平台基于Node.js开发的Koa router路由插件代码信息详情
- c#下载文件选择路径控件
- Asp.net core 学习笔记 Fluent Validation
- MySQL创建新用户以及ERROR 1396 (HY000)问题解决
- Oracle sqlplus失去响应解决方法/如何在数据库失去响应时转储状态信息(转)
- summernote 如何设置为只读?
- linux_jdk_mysql_tomcat
热门文章
- P2057 善意的投票 最小割理解
- PC客户端自动化测试工具--pywinauto
- [SDOI2010]代码拍卖会——DP
- C# 通过Process.Start() 打开程序 置顶方法
- @GetMapping、@PostMapping、@PutMapping、@DeleteMapping、@PatchMapping、@RequestMapping详解
- python 省略号 三个点...的含义
- JS 中获取服务器时间的注意点
- learning armbian steps(9) ----- armbian 源码分析(四)
- HDU 5544 Ba Gua Zhen ( 2015 CCPC 南阳 C、DFS+时间戳搜独立回路、线性基 )
- gzip/bzip/xz/tar