print agm(1, 1/sqr(2))
print using("#.#################",agm(1, 1/sqr(2)))
end

function agm(a,g)
    do
        absdiff = abs(a-g)
        an=(a+g)/2
        gn=sqr(a*g)
        a=an
        g=gn
    loop while abs(an-gn)< absdiff
    agm = a
end function