import random

def compute_pi(N):
    x = 0.0
    sum = 0.0

    while x < 1.0:
        sum += (1.0 - x*x) ** 0.5
        x += 1.0/N

    return 4.0 * sum / N

print compute_pi(100000000)
