Hoje é o dia do Pi. Isso mesmo, 3/14 em notação americana.
O número pi (3,141592….) é um número obtido da relação entre o diâmetro de um círculo e seu perímetro . Nem precisava dizer isto pois os leitores do Zeletron sabem isto.
É um número fascinante: irracional e transcendente.
Aqui vai um programa Python para calcular PI
#!/usr/bin/env python from __future__ import with_statement import decimal def pi_gauss_legendre(): D = decimal.Decimal with decimal.localcontext() as ctx: ctx.prec += 2 pi_fun = lambda a,b,t: (a + b)**2 / (4 * t) a, b, t, p = 1, 1/D(2).sqrt(), 1/D(4), 1 pi, piold = pi_fun(a, b, t), None while pi != piold: piold = pi an = (a + b) / 2 b = (a * b).sqrt() t -= p * (a - an)**2 a, p = an, 2*p pi = pi_fun(a, b, t) return +pi decimal.getcontext().prec = 100 print pi_gauss_legendre() |
2 comentários em “Dia do Pi”