
import math

g = 9.8
def w(L):
    return math.sqrt(g / L)

def k(L, b):
    denom = math.sqrt(g) + b * math.sqrt(L)
    return g / (denom ** 2)
    

# Let's say our first pendulum is 1 meter long, w is 3.13 radians per second
# and with 15 more pendulums we want to reach about 9.13, so 6 / 15 is 0.4
L = 1
print L, w(L)
for i in range(15):
    L = L * k(L, 0.4)
    print "L: %s, w: %s" % (L, w(L))

