def get_asymptotic_parameters(numax, dnu, teff, d01_mass_offset=-0.073, d01_mass_slope=0.044, d01_offset=-0.063, d02_mass_offset=0.138, d02_mass_slope=-0.014, d02_offset=0.035, d03_slope=0.282, d03_offset=0.16, numax_sun=3150., dnu_sun=134.9, teff_sun=5777):
    """
    Compute relevant asymptotic parameters from `numax`, `dnu`, and `teff`.

    This function provides estimates of stellar mass and radius, as well as the 
    small frequency spacings d02, d01, d03.

    NOTE: The asymptotic parameters d02, d01, d03, are reliable only for RG 
    stars. When using this function for SG and MS, make sure that proper 
    corrections are applied.

    Parameters
    ----------
    numax : float
        The frequency of maximum oscillation power in microHz.
    dnu : float
        The large frequency spacing in microHz.
    teff : float
        Efffective temperature of the star in Kelvin.
    d01_mass_offset : float, default: -0.073
    d01_mass_slope : float, default: 0.044
    d01_offset : float, default: -0.063
    d02_mass_offset : float, default: 0.138
    d02_mass_slope : float, default: -0.014, 
    d02_offset : float, default: 0.035
    d03_slope : float, default: 0.282
    d03_offset : float, defualt: 0.16
    numax_sun : float, default: 3150.
        Solar reference value for frequency of maximum oscillations (`numax`).
    dnu_sun : float, default: 134.9 
        Solar reference value for large frequency spacing (`dnu`).
    teff_sun : float, default: 5777
        Solar reference vaule for effective temperature (`teff`).

    Returns
    -------
    asymptotic_parameters : dict
        A dictionary of the computed asymptotic values. The keys are `'mass'`, 
        `'radius'`, `'d01'`, `'d02'`, and `'d03'`. 
    """
    # Raw mass estimate from scaling relation
    mass = (numax/numax_sun)**3 * (dnu/dnu_sun)**-4 * (teff/teff_sun)**1.5

    # Raw radius estimate from scaling relation
    radius = (numax/numax_sun) * (dnu/dnu_sun)**-2 * (teff/teff_sun)**0.5

    # Scaling relation Dnu - d02 (Corsaro et al. 2012b)
    slope = d02_mass_offset + d02_mass_slope*mass
    d02 = slope*dnu + d02_offset

    # Scaling relation Dnu - d01 (Corsaro et al. 2012b)
    slope = d01_mass_offset + d01_mass_slope*mass
    d01 = slope*dnu + d01_offset

    # Scaling relation Dnu - d03 (Huber et al. 2010)     [1,20] microHz
    d03 = d03_slope*dnu + d03_offset

    asymptotic_parameters = {'d02':     d02,    
                             'd01':     d01,   
                             'd03':     d03,   
                             'mass':    mass,  
                             'radius':  radius,}

    return asymptotic_parameters