def METHOD_NAME(fluxes_obs, fluxes_model):
    '''
    Computes the residuals of the input fluxes and best fit model
    Parameters
    ----------
    fluxes_obs: array-like
        Observed fluxes
    fluxes_model: array-like
        Model fluxes
    Returns
    -------
    residuals_mean: float
        The mean of the residuals
    residuals_stdev: float
        The standard deviation of the residuals
    '''
    residuals = fluxes_obs - fluxes_model
    return np.mean(residuals), np.std(residuals