前に「C言語による実用アルゴリズム入門」からshort intの配列を使って、高精度の数値計算をやって見てました。
勿論それはそれで良いんでしょうが、Rubyのマニュアルにπの計算が載ってました。それによれば、BigDecimalを
使って、小数点以下を計算してました。
何とかできないものか、四苦八苦してみました。上の図の後半部分の数行で何とかなりましたよ。特別なテクニックも
要りませんでした。驚きです。( ;∀;)
仕上げはループ回数です。やみくもにループしないで、収束したら終わるですね。10の-1001乗で1000桁まで
誤差が出ないようです。