#
# pi.py
# (c) Neil Gershenfeld 2/10/08
# calculation of pi by a scalar sum
# pi = 3.14159265358979323846 
#

from numpy import *

NPTS =  2000000
NLOOP = 1
a = 0.5
b = 0.75
c = 0.25
for loop in range (NLOOP):
   i = arange(1,(NPTS+1),dtype=float64)
   x = sum(0.5/((i-0.75)*(i-.25)))
print "pi.py: NPTS, nloop = %d, %d"%(NPTS,NLOOP)
print "pi = %.20f"%x