def METHOD_NAME(self, P,M_min,M_max,b):
    """
    Returns the inverse of the Doubly Truncated Gutenberg Richter distribution
    for the specified probability, magnitude range, and coefficient.
    """
    B = b*np.log(10)
    M = np.log(np.exp(-B*M_min) - P*(np.exp(-B*M_min)-np.exp(-B*M_max)))/(-B)
    return M