圆周率计算（算术几何平均）
使用算术几何平均（AGM公式）计算圆周率π。
当上次和本次的计算结果相等时，程序即结束。增加运算位数可提高π的精度。
1976年，萨拉明和布伦特发现，使用高斯公式(1809)可以使圆周率以极快的速度实现二次收敛。
这种收敛性表现为：随着循环次数的增加，精度位数会以5位、10位、20位为单位成倍增长。
1976年，萨拉明和布伦特发现，使用高斯公式(1809)可以使圆周率以极快的速度实现二次收敛。
这种收敛性表现为：随着循环次数的增加，精度位数会以5位、10位、20位为单位成倍增长。
| \(\normalsize Gauss-Legendre\ method\ 1809\\ \hspace{80px}\pi={\large\frac{2\hspace{1px}{\rm AGM}^2\left(1,{\large\frac{1}{\sqrt{2}}}\right)}{1-{\large\displaystyle \sum_{\small k=0}^{\small \infty}}2^kc_k^2}}\\\\ The\ Square\ {\rm AGM}\\ \hspace{80px} by\ \ Salamin\ \&\ Brent,\ 1976\\ (1)\ a_0=1,\hspace{20px}b_0={\large\frac{1}{\sqrt{2}}},\hspace{20px}t_0={\large\frac{1}{4}}\\ (2)\ a_{n+1}={\large\frac{1}{2}}(a_n+b_n),\hspace{20px}b_{n+1}=\sqrt{a_nb_n},\\ \hspace{50px}t_{n+1}=t_n-2^n(a_n-a_{n+1})^2\\ (3)\ \pi=\displaystyle \lim_{n\to\infty}{\large\frac{(a_n+b_n)^2}{4t_n}}\\ \normalsize The\ Quartic\ {\rm AGM}\ by\ J.M.\ Borwein\\ \hspace{120px} \&\ P.B.\ Borwein,\ 1985\\ (1)\ y_0=\sqrt{2}-1,\hspace{20px}a_0=6-4\sqrt{2}\\ (2)\ y_{n}={\large\frac{1-\sqrt[4]{1-y_{n-1}^4}}{1+\sqrt[4]{1-y_{n-1}^4}}},\\ \hspace{10px}a_{n}=(1+y_{n})^4a_{n-1}-2^{2n+1}y_{n}(1+y_{n}+y_{n}^2)\\ (3)\ \pi={\large\displaystyle \lim_{\small n \to\infty}\frac{1}{a_n}}\\ \) |
目前还没有关于此计算的用户评价。