In [26]:
%matplotlib inline
In [20]:
import numpy as np
from numpy.random import randint
from itertools import count
import matplotlib.pyplot as plt
In [36]:
probability = []
ntrials = 100000
people_in_group = range(2,70)
for number_of_people in people_in_group:
cnt_common_birthdays = 0
for i in range(ntrials):
random_birthdays = randint(365, size=number_of_people)
if len(set(random_birthdays)) < number_of_people:
cnt_common_birthdays += 1
probability.append(cnt_common_birthdays*100/float(ntrials))
In [37]:
plt.plot(people_in_group, probability)
Out[37]:
In [ ]: