#!/usr/bin/env python from __future__ import division from math import sqrt RADIUS = 1000 def quarter_circle_method(): dots = 0 for x in xrange(RADIUS): for y in xrange(RADIUS): if sqrt( x**2 + y**2 ) <= RADIUS: dots += 1 Pi = (4*dots) / (RADIUS**2) return Pi if __name__ == "__main__": print quarter_circle_method()