【研究】 筑波大がスパコンで７２時間かけ樹立した円周率の桁数世界記録、仏の技術者が２７万円のＰＣで１３１日かけ破る★５

筑波大学がπ計算で使ってたプログラムのベースになってるガウス＝ルジャンドルのアルゴリズム
これは200年も前に考案されたアルゴリズム。単純だけど、わずか3ループしただけでも
倍精度の精度限界までπを求められる驚速アルゴリズム
昔の天才数学者すげー
//──────────────────────────────
//Windows環境なら以下をメモ帳にコピペして、拡張子を"js"に変更して実行
//ガウス＝ルジャンドルのアルゴリズム

var loop = 3;//ループ回数　これ以上増やしても無駄だから1や2に減らしてみてもいい
var a = 1;
var b = Math.sqrt(2) / 2;
var t = 0.25;
var p = 1;
var x = 0;
var y;
var pie;
var i;

for (i = 0; i < loop; i++){
x = (a + b) / 2;
y = Math.sqrt(a * b);
t -= p * (a - x) * (a - x);
p = 2 * p;
a = x;
b = y;
}

pie = (a + b) * (a + b) / (4 * t);

var WSHShell = new ActiveXObject("WScript.Shell");
WSHShell.Popup(pie, 0, "ガウス＝ルジャンドル", 0);
//──────────────────────────────