def model(x, theta):
    '\n        theta = (beta, kappa)\n        '
    return np.log(schmidt_law(x, theta))
theta = (beta, kappa)
def model(x, theta):
    '\n        \n        '
    return np.log(schmidt_law(x, theta))