みなさん、こんにちは！
オウミ技研、ハード部です！！

以前に、『なぜ円周率は「３．１４・・・」なのか？』という投稿記事を掲載しましたが、
予想に反して、毎日、多数のアクセスを頂いています（幅広い世代の疑問だからでしょうか？）。

読んでいない方は ⇒ こちらから！

この意外にも？ヒットした投稿の続編として、

　『疑似乱数で円周率の近似値を求める』

これを今回の投稿テーマとしたいと思います。

早速、本題に入りますが、
『疑似乱数』をどのように使って『円周率の近似値』を算出するのか？
先ずは、その手順を説明したいと思います。

【手順①】
　・ｘｙグラフの原点を中心に半径１の円を描きます。

【手順②】
　・ｘ、ｙ各々に『０～１の間の疑似乱数』を生成します。
　・乱数が生成されたｘ、ｙ座標をグラフ上にプロットします（赤枠内にプロットされます）。
　・上記を既定回数、繰り返します。

【手順③】
　・『グラフの原点』から『プロットされた各々の座標』までの距離を以下の式から算出します。

　　☆ 原点からプロットされた座標までの距離＝√（ｘ²＋ｙ²）
　・『√（ｘ²＋ｙ²）』が『１（＝円の半径）』の場合、その点は扇の円弧上にプロットされます。
　・『√（ｘ²＋ｙ²）』が『１』より小さい場合、その点は扇の円弧より内側にプロットされます。
　・『√（ｘ²＋ｙ²）』が『１』より大きい場合、その点は扇の円弧より外側にプロットされます。

【手順④】
　・赤枠の正方形と斜線の扇の面積を算出します。

　　☆ 正方形の面積＝１＊１＝１
　　☆ 扇の面積＝πｒ²＊１／４＝π／４

　　☆ 『正方形の面積』と『扇の面積』の比率＝１：π／４
　　　 ・・・①
　　☆ 『正方形内にプロットされた座標の個数』と『扇内（扇の外周も含む）にプロットされた座標の個数』の比率
　　　 ・・・②

　・上記①、②がイコールとなることから、

　　☆ １：π／４＝『ｘ、ｙ座標の疑似乱数の総数』：『扇内にプロットされたｘ、ｙ座標の疑似乱数の個数』
　　☆ ∴ π＝『扇内にプロットされたｘ、ｙ座標の疑似乱数の個数』＊４／『ｘ、ｙ座標の疑似乱数の総数』

　・なお、『扇内にプロットされた疑似乱数』をカウントする際の条件式は『√（ｘ²＋ｙ²）≦１』となります。

以上の『手順①～④』をもって、疑似乱数による円周率の近似値の算出方法が理解できたと思います。
次の段階としては、
実際に、生成した疑似乱数を散布図にプロット表示した上で、円周率の近似値を算出したいと思いますが・・・
次回の投稿まで、お待ち下さい！！

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
　★ ハード部のことをもっと知りたい ⇒ ハード部専用ブログへ！
　★ オウミ技研のＩｏｔに対する取組みを知りたい ⇒ ブログ「中小企業へのＩｏｔ導入に向けて１／２」
　 　　　　　　　　　　　　　　　　　　　　　 　⇒ ブログ「中小企業向け『簡易Ｉｏｔ』１／２／３／４」
　★ ハードウェア・ソフトウェア設計・製造に関するお悩みがございましたら、ぜひお気軽に お問い合わせ を！
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝