圆周率的计算（高斯-勒让德的算术几何平均法）

我们将验证筑波大学的高桥等人用于计算出精确度达2兆5769亿8037万位圆周率的AGM法（高斯-勒让德算法，即算术几何平均法）。
以a0=1、b0=1/√2、t0=1/4、p0=1为初始值，
an+1 = (an + bn) / 2
bn+1 = √(anbn)
tn+1 = tn - pn * (an+1 - an)^2
pn+1 = 2*pn
据此，
π≒(an + bn)^2 / (4*tn)
进行计算。