from math import sqrt

def arithmeticMean (a,b):
    c=(a+b)/2
    return c

def geometricMean (a,b):
    c=sqrt(a*b)
    return c

def harmonicMean (a,b):
    c=1/(arithmeticMean(1/a,1/b))
    return c

def archimedes(iterations):
    outer = 6/sqrt(3)
    inner = 3
    sides = 6
    for i in range(iterations):
        outer = harmonicMean(inner, outer)
        inner = geometricMean(inner, outer)
        sides = sides * 2
        print(sides, inner, outer)