def METHOD_NAME(self, magnitudes, delta_m):
    '''If a single magnitude spacing is input then create the bins
    :param numpy.ndarray magnitudes:
        Vector of earthquake magnitudes
    :param float delta_m:
        Magnitude bin width
    :returns: Vector of magnitude bin edges (numpy.ndarray)
    '''
    min_mag = np.min(magnitudes)
    max_mag = np.max(magnitudes)
    if (max_mag - min_mag) < delta_m:
        raise ValueError('Bin width greater than magnitude range!')
    mag_bins = np.arange(np.floor(min_mag), np.ceil(max_mag), delta_m)
    # Check to see if there are magnitudes in lower and upper bins
    is_mag = np.logical_and(mag_bins - max_mag < delta_m,
                            min_mag - mag_bins < delta_m)
    mag_bins = mag_bins[is_mag]
    return mag_bins