Calculando π con Python

¿Cuántos decimales creéis que se podrían calcular en 1 segundo? ¿Y en 0,024 s? Por supuesto esto depende de la potencia de vuestro ordenador. Si queréis comprarlo aquí os dejo un programa en Python que permite calcular los decimales de π que queramos.

#!/usr/bin/env python
# -⁻- coding: UTF-8 -*-
print "   πππππππππππ"
print "   π π     π"
print "     π     π"
print "     π     π"
print "  Pi Calcutator"
print 

def arccot(x, unity):
sum = xpower = unity // x
n = 3
sign = -1
while 1:
xpower = xpower // (x*x)
term = xpower // n
if not term:
break
sum += sign * term
sign = -sign
n += 2
return sum

def pi(digits):
unity = 10**(digits + 10)
pi = 4 * (4*arccot(5, unity) - arccot(239, unity))
return pi // 10**10

print "Inserte el número de decimales que desea calcular: "
decimales=int(raw_input())
print pi(decimales)

Si además queremos saber el tiempo ejecutamos el programa con el comando time delante. Por ejemplo, de esta forma:

time ./pi.py

Calculando π

Y a vosotros, ¿cuánto os tarda en calcular 1000 decimales?
Filed under: GNU/Linux, Programación, Python, Software libre