功能: 求算术几何平均(Arithmetic-Geometric Mean)

格式: AGM(a, b)

a : 为数值或者矩阵变量
b : 为数值或者矩阵变量

说明:
1. a,b里的元素值不能小于0,否则返回非数字
2. 本函数本质求如下表达式
$$
\begin{aligned}
\mathbf{AGM}(a,b) &= \dfrac{\pi}{2}\int_0^\infty \dfrac{1}{\sqrt{(t^2+a^2)(t^2+b^2)}}dt
\\
\\
a_1 &=a,b_1=b
\\
&\begin{cases}
a_i &= \dfrac{a_{i-1}+b_{i-1}}{2}
\\
b_i &=\sqrt{a_{i-1}b_{i-1}}
\end{cases},i=2,3,4,...,n 
\\
\mathbf{AGM}(a,b) &= a_n
\end{aligned}
$$

参考:{ACM<https://mathworld.wolfram.com/Arithmetic-GeometricMean.html>}

例子:
a = AGM(4,3)
a =
[ 3.48202767635957 ]