@staticmethod
def evaluate(x, amplitude, x_break, alpha_1, alpha_2):
    'One dimensional broken power law model function'
    alpha = np.where((x < x_break), alpha_1, alpha_2)
    xx = (x / x_break)
    return (amplitude * (xx ** (- alpha)))
One dimensional broken power law model function
@staticmethod
def evaluate(x, amplitude, x_break, alpha_1, alpha_2):
    
    alpha = np.where((x < x_break), alpha_1, alpha_2)
    xx = (x / x_break)
    return (amplitude * (xx ** (- alpha)))