#
# pi = 3.14159265358979323846
#

from numarray import *

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