def multi_phase_options(printout=True):
    """Retrieve all the options for multiple scattering radiation"""
    if printout:
        print(
            "Can also set delta_eddington=True/False in approx['delta_eddington']"
            )
    return ['N=2', 'N=1']
