用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')
最后显示结果如下:
以上就是本文的全部内容,希望对大家的学习有所帮助。