モンテカルロ法による円周率算出 (2020年10月)

0~1.0までの間の浮動小数のランダム値を二つ得て、1.0 x 1.0の正方形内に仮想的にプロット。 x*x+y*y=1の方程式により、x*x+y*y<=1なら円内としてカウント(C)。 *4してそれを回数(N)で割ると円周率の近似値になります。 近似解π=4C/N。Nが多いほど円周率に…