def agm(x, y, n):
    a = (x + y)/2
    g = (x * y)**0.5

    for x in range(1, n):

        a, g = (a + g)/2, (a * g)**0.5

    return [a, g]