def maxVelocity(l1, l2, v2):
    """
    This function takes perihelion, aphelion distances and minimum velocity
    to calculate the maximum velocity at the aphelion in AU/year.
    This equation is taken from Lab4 handout.
    """
    return (l2*v2)/l1