import random inside = 0 nsamples = 120000 for i in range(nsamples): x = random.random(); y = random.random(); if (x*x)+(y*y)<1: inside += 1 pi = (4.0 * inside)/nsamples print "Computed value of pi is",pi