scipy.special.agm

scipy.special.agm(a, b)

Arithmetic, Geometric Mean

Start with a_0=a and b_0=b and iteratively compute

a_{n+1} = (a_n+b_n)/2 b_{n+1} = sqrt(a_n*b_n)

until a_n=b_n. The result is agm(a,b)

agm(a,b)=agm(b,a) agm(a,a) = a min(a,b) < agm(a,b) < max(a,b)