import numpy

def compute_pi(N):
    sum = 0.0
    # Create a vector of n values between 0 and 1
    xs = numpy.linspace(0.0, 1.0, N)
    # compute the values of sqrt(1 - x**2)
    vs = numpy.sqrt(1.0 - xs**2)
    sum += vs.sum()
    return 4.0 * sum / N

print compute_pi(100000000)