buffon.py

# coding: utf-8

import random
import math

def fall():
    rad = random.random() * math.radians(360)
    c = random.random() +  math.fabs(math.cos(rad))
    if c > 1:
        return 1
    else:
        return 0


N = 1000000
cnt = 0
all = 0
while True:
    all += N
    for i in xrange(N):
        cnt += fall()
    print(cnt)
    print(all)
    print(float(cnt)/all)