2
0
Fork 0
mirror of https://github.com/MartinThoma/LaTeX-examples.git synced 2025-04-26 06:48:04 +02:00
LaTeX-examples/tikz/birthday-paradox/calculate.py
2012-08-09 07:04:38 +02:00

14 lines
328 B
Python

#!/usr/bin/python
# -*- coding: utf-8 -*-
from math import factorial
from gmpy import bincoef
f = open('data.csv', 'w')
f.write('People\tprobability\n')
def prob(people):
return 1.0 - float(factorial(people)*bincoef(365,people))/(365**people)
for people in xrange(60+1):
f.write("%i\t%f\n" % (people, prob(people)))