代码之家 › 专栏 › 技术社区 › Alexey Kirov

莱布尼茨数行的Python脚本

pi python

Alexey Kirov · 技术社区 · 9 年前

我试图用Python 2.7编写简单的脚本来计算Pi。我对莱布尼茨公式感兴趣:

= (4/1) - (4/3) + (4/5) - (4/7) + (4/9) - (4/11) + (4/13) - (4/15) ...

所以我认为“(4/1)-(4/3)”是X,1-3-5是y,y+2,y+4。。。。

我写了这个脚本:

def pi():
    p = 0.0
    y = 1.0
    x = 4.0/y - 4.0/(y+2.0)
    for i in range(10000):
        p = p + x
        y += 4.0
    print p

pi()

2 回复 | 直到 9 年前

Ukimiku 9 年前

这个工作原理:

def pi():
    p = 0.0
    y = 3.0
    x = 4.0 - 4/y
    for i in range(5000):
        y += 2
        x = x + 4/y
        y += 2
        x = x - 4/y
        print x

pi()

Keiwan Thirler 9 年前

x 和 y 方法

下面是使用实际公式的替代方法:

def pi():
    p = 0.0
    for i in range(10000):
        p += 4 * ((-1)**i) / float(2 * i + 1)

    print p

对于 range(10000) 这将打印

3.14149265359

推荐文章

Google User · Django管理员在`list_display中未显示`creation_date`字段`

9 月前

user29747013 · 如何创建一个新的数据框架,其中包含原始数据框架中列的聚合列?

9 月前

ÎÎÎ½Î· ÎÎ®Î¹Î½Î¿Ï · Python lxml.html语法错误:使用lxml find时XPATH的谓词无效

9 月前

user29715306 · from_users=和chats=电视节目中的差异

9 月前

Redshoe · 当执行numpy.genfromtxt()时,python是否会读取文件的所有行?

9 月前

RASEL MAHMUD · 为什么以及如何在is_even()函数内的IF条件中递归X变量在满足0后递增?[副本]

9 月前

prayner · 更新嵌套字典包含列表中的项

9 月前

Bringo Jr · 我可以在O(n)中解决这个问题吗?

9 月前

Dave · 如何在for循环中修改列表值

9 月前

Shukurullox Komiljonov · 从记录中获得相互和解。使用SQL

9 月前