def plot_expected_limit_bands(axis, bands):
    """Add the expected limit bands to a plot

    Parameters
    ----------

    axis : `matplotlib.Axes`
        The axes we are adding the bands to

    bands : dict
        Dictionary with the bands

    """
    masses = bands['masses']
    axis.fill_between(masses, bands['q02'], bands['q97'], color='yellow')
    axis.fill_between(masses, bands['q16'], bands['q84'], color='green')
    axis.plot(masses, bands['median'], color='gray')
