def METHOD_NAME(
    distance_categories_to_highest_spectrum,
    constant,
    coefficients,
    population,
    score,
    real_numbers,
    output_unmet_demand,
    vector_base_map=False,
    vector_methods=None,
    vector_column_prefix=None,
):
    """
    Parameters
    ----------
    distance_categories_to_highest_spectrum:
        Map of distance categories to areas of highest recreation spectrum
    constant:
        Constant for the mobility distance function
    coefficients :
        Set of coefficients for distance > 4 category, as
        (the example) presented in the following table:
        |----------+---------+---------|
        | Distance | Kappa   | Alpha   |
        |----------+---------+---------|
        | >4       | 0.06930 | 0.00057 |
        |----------+---------+---------|
    population:
        A map with the distribution of the demand per distance category and
        predefined geometric boundaries (see `r.stats.zonal` deriving the
        'demand' map ).
    score :
        A score value for the mobility function
    real_numbers :
        If real_numbers is False (which is the default), the mapcalc expression
        used to build the mobility function, will derive integer values.
    vector_base_map:
        A 'base' vector map to update with columns for unmet demand statistics
    vector_methods:
        Descriptive statistics for the `v.rast.stats` call--see update_vector()
        function
    vector_column_prefix:
        Prefix for 'unmet demand' columns added in the 'base' vector map
    Returns
    -------
    This function does not return anything
    """
    unmet_demand_expression = build_unmet_demand_expression(
        distance=distance_categories_to_highest_spectrum,
        constant=constant,
        coefficients=coefficients,
        population=population,
        score=score,
        real_numbers=real_numbers,
    )
    # suitability=suitability)  # Not used.
    # Maybe it can, though, after successfully testing its
    # integration to build_distance_function().
    msg = "*** Unmet demand function: {f}"
    grass.debug(_(msg.format(f=unmet_demand_expression)))
    unmet_demand_equation = EQUATION.format(
        result=output_unmet_demand, expression=unmet_demand_expression
    )
    r.mapcalc(unmet_demand_equation, overwrite=True)
    if vector_base_map:
        update_vector(
            vector=vector_base_map,
            raster=output_unmet_demand,
            methods=METHODS,
            column_prefix=COLUMN_PREFIX_UNMET,
        )