用python算圆周率π
1.准备第三方库pip
打开cmd 输入代码:pip install requests ,随后就会成功
因为小编已经安装好了,所以就不把图截出来了
2.利用马青公式求π
3.用python语言编写出求圆周率到任意位的程序如下:
from math import *
from tqdm import tqdm
from time import *
total,s,n,t=0.0,1,1.0,1.0 clock()
while(fabs(t)>=1e-6):
total+=t n+=2 s=-s t=s/n k=total*4 print("π值是{:.10f}
运行时间为{:.4f}秒".format(k,clock())) for i in tqdm(range(101)): print("\r{:3}%".format(i),end="") sleep((clock())/100)#用执行程序的总时间来算出进度条间隔的时间
4.在cmd中运行如图所示:
在IDLE中运行执行如下代码:
#- * -coding: utf - 8 - * - from __future__ import division #################### 导入时间模块 import time ############### 计算当前时间 time1 = time.time() ################ 算法根据马青公式计算圆周率#################### number = int(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 range(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')
最后显示结果如下:
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持亿速云。
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。