def METHOD_NAME(areas, ctoas, seed):
    # S=pi*a*b
    # C=2pib + 4(a-b)
    # R = C/S
    # x = RS
    # y = S/pi
    # long = np.sqrt(y-np.pi*y/2+x/4)
    # short = y/longs
    x = ctoas * areas
    y = areas / np.pi
    longs = np.sqrt(y - np.pi * y / 2 + x / 4)
    shorts = y / longs
    axes = list(zip(longs, shorts))
    return axes