'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)))
