#
# pi = 3.14159265358979323846
#

from numarray import *

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