def fit_lmfit_schmidt(x, y, yerr, init=None):
    (m, _) = lmfit_powerlaw(x, y, yerr, init=init)
    return m.parameters
def fit_lmfit_schmidt(x, y, yerr, init=None):
    (m, _) = lmfit_powerlaw(x, y, yerr, init=init)
    return m.parameters