円の面積から円周率PIEを出してみる。シンプソン法。

上は台形の公式を使って求めたもの。四季が単純で分かりやすいかも。

　X軸は0から100まで、区分は120個。C言語のサンプルがありました。ネットから拝借しました。3.14まで合ってます。一番驚いたのが、PDFファイルからCOPYできるようになっていたことです。

分割数を250個にあげてみました。3.141で一桁上がりました。

分割数を500に上げても3.141までしか合ってません。

分割数を1000に上げても、3.1415止まり。C言語では小数点以下6桁が限度かも。

　上のC言語のサンプルを参考にRubyでやってみました。for i in rangeでiは単に繰り返しで使ってるだけで、それでもOKのようです。r=100000にして、分割数200000。繰り返しのrange = 1..199999になることに注意ですかね。この数値の組み合わせでエラーになることもありました。追求してません。3.14159265まで合ってますかね。