遅れたのでビュフォンの針でお詫び

本日の営業は終了いたしました

って、今日の投稿忘れてた！！！！！

お詫びのしるしとしてビュフォンの針から円周率をたたき出すプログラムを作成したので、これでご勘弁を(´；ω；`)ｳｩｩ

ちょっと重いので時間はかかりますがそこそこの精度が出ます

import math
import random
p = 100000000
x,t=0,0
for i in range(p):
    x =random.random()-math.sin(random.random()*math.pi)
    if x<=0:
        t+=1
print(p/t*2)

実行結果

一回目：3.1415927146213622
二回目：3.141645912691964
三回目：3.1416860345065514
四回目：3.1414104829401857
五回目：3.1414468486356038
六回目：3.1417024191092917
七回目：3.1409505813671807
八回目：3.1414401379418933
九回目：3.1414525724852482
十回目：3.1414790702334034
十一回目：3.1415133156576087
十二回目：3.1418585749230137
十三回目：3.1413903020694947

さて、以上の結果を分析～

平均
3.141505305
不偏分散
4.70717E-08

このくらいの精度なので実用には向いていませんが、
運が良ければ一回目のようにうまくいきます。