誕生日が一致する確率（N人中k人以上）をモンテカルロシミュレーションで計算

昨日は近似式で計算してみた。

さて、今日はモンテカルロシミュレーションで確認しましょう。乱数はもちろん、メルセンヌツイスタを使用。N=1000000試行して、確率を求めた。Excel　VBAを使用してますよ。

で比較。

近似式の精度、非常にいいですね。