mata

real scalar agm(real scalar a, real scalar b) {
	real scalar c
	do {
		c=0.5*(a+b)
		b=sqrt(a*b)
		a=c
	} while (a-b>1e-15*a)
	return(0.5*(a+b))
}

agm(1,1/sqrt(2))
end