@staticmethod
def evaluate(mag, phi_star, m_star, alpha):
    'Schechter luminosity function model function.'
    if (isinstance(mag, Quantity) or isinstance(m_star, Quantity)):
        raise ValueError('mag and m_star must not have units')
    factor = (10 ** (0.4 * (m_star - mag)))
    return ((((0.4 * np.log(10)) * phi_star) * (factor ** (alpha + 1))) * np.exp((- factor)))
Schechter luminosity function model function.
@staticmethod
def evaluate(mag, phi_star, m_star, alpha):
    
    if (isinstance(mag, Quantity) or isinstance(m_star, Quantity)):
        raise ValueError('mag and m_star must not have units')
    factor = (10 ** (0.4 * (m_star - mag)))
    return ((((0.4 * np.log(10)) * phi_star) * (factor ** (alpha + 1))) * np.exp((- factor)))