|
发表于 2018-6-20 07:29:35
|
显示全部楼层
下面是公式 pi = 16 * arctan(1/5) - 4 * arctan(1 / 239) 的Python计算程序:
# Calculate PI, 2018-6-20
# pi = 16 * arctan(1/5) - 4 * arctan(1 / 239)
P5_2 = 1.0 / (5.0 * 5.0)
P239_2 = 1.0 / (239.0 * 239.0)
t5_2 = 1.0 / 5.0
t239_2 = 1.0 / 239.0
s5 = t5_2
s239 = t239_2
for i in range(2, 20):
t5_2 = -t5_2 * P5_2
s5 = s5 + t5_2 / ( 2 * i - 1.0)
t239_2 = -t239_2 * P239_2
s239 = s239 + t239_2 / ( 2 * i - 1.0)
pi_ = 16.0 * s5 - 4.0 * s239
print(pi_)
运行结果:
3.141592653589794
只用20次循环,就精确到小数点后14位。对比你那个是20万次循环,精确到小数点后5位
|
|