Мне нужна помощь в вычислении Pi. Я пытаюсь написать программу python, которая будет вычислять цифры Pi до X. Я попробовал несколько из списка рассылки python, и он должен замедляться для моего использования. Я прочитал о Gauss-Legendre Algorithm, и я попытался перенести его на Python без успеха.
Я читаю из здесь, и я был бы признателен за любые данные о том, где я ошибаюсь!
Выводится: 0.163991276262
from __future__ import division
import math
def square(x):return x*x
a = 1
b = 1/math.sqrt(2)
t = 1/4
x = 1
for i in range(1000):
y = a
a = (a+b)/2
b = math.sqrt(b*y)
t = t - x * square((y-a))
x = 2* x
pi = (square((a+b)))/4*t
print pi
raw_input()