モンテカルロ法でπを求めます。
モンテカルロ法は、正方形に内接する円を考え、乱数で正方形内に任意の点を打ちます。
円の半径をｒとすると、円の面積はπｒ２で正方形の面積は４ｒ２ですので面積比は π：４ となり、円の内側にある点の数（ｍ）と全ての点の数（ｎ）を数えれば、π/４＝ｍ/ｎ からπが求まることになります。ここではｒ＝１とします。
表のＸ，Ｙは０〜１の乱数です。=RAND() の関数を使います。４行目〜６５５３６行目（エクセルの最大行）まで入力します。
Ｘ２＋Ｙ２を求めます。
円の中かどうかは Ｘ２＋Ｙ２＜ｒ２（＝１） なら円の中ですので、判定結果の列に =IF(E4<1,"",0) を入れて円の外なら０を表示するようにします。（円の外の方が数が少ないので処理が早い）
判定結果の列の上に =COUNT(F4:F65536) で判定結果の列の４行目から６５５３６行目（エクセルの最大行）までにいくつ数字があるか数えます。これが円の外の数です。全体（６５５３３）から引けば円の中の数が求まります。円の中の数を全体で割って４倍するとπが求まります。
さて実行結果は
３．１４４６１４
３．１３７１６８
３．１４４０６５
３．１４９００９
３．１３８６９３
３．１４１３７９
などが得られました。
判定を「＜（未満）」としたのですが、線上はどうなるのかと「＜＝（以下）」の判定でやっても結果は変わりません。線上はそれほどはないのでしょうね。
再計算の処理時間は「＜」のとき７分４０秒なのが「＜＝」とすると１１分２５秒と処理時間が増えてしまいます。判定に時間がかかるのでしょうね。
また、外側を数えるのではなく内側を数えると１１分２４秒と処理時間が増えてしまいます。処理回数が増えるからでしょうね。
ちなみに、ＣＰＵはアスロンＸＰの＋１７００でメモリーは５１２Ｍです。