def get_number_density(m, mfunc, dm=None):
    """
    Performs integration of Halo Mass Function.
    
    Args:
        m (np.ndarray): Array of Exponent of Halo Masses in M_sun.
        mfunc (np.ndarray): Array of Corresponding Halo Mass Functions
        dm (float): Exponent of Mass Resolution for integration.
                    [i.e., default dM = 10^7 to integrate dn/dln(M)]
    """

    if dm == None:
        # Assumes Linearly-Spaced Data
        dm = m[1] - m[0]
        
    return sum(mfunc/m) * 10**dm