def METHOD_NAME(areas, ltos, seed):
    # S = pi*a*b
    # R = a/b
    # b = np.sqrt(S/(R*pi))
    # a = np.sqrt(S/(R*pi)) * R
    np.random.seed(seed)
    while len(areas) > len(ltos):
        ltos = np.tile(ltos, 2)
    ltos = ltos[0 : len(areas)]
    shorts = np.sqrt(areas / (ltos * np.pi))
    longs = (shorts * ltos).astype(np.uint16)
    shorts = shorts.astype(np.uint16)
    axes = list(zip(longs, shorts))
    return axes