How to use the impedance.models.circuits.fitting.rmse function in impedance

ts[k-1] = 10**(np.log10(t_min) +
                               ((k-1)/(M-1))*np.log10(t_max/t_min))

        return ts

    R0 = min(np.real(Z))
    if c is not None:
        M = 0
        mu = 1
        while mu > c and M <= max_M:
            M += 1
            ts = get_tc_distribution(f, M)
            p_values, mu = fitLinKK(f, ts, M, Z)

            if M % 10 == 0:
                print(M, mu, rmse(eval_linKK(p_values, R0, ts, f), Z))
    else:
        M = max_M
        ts = get_tc_distribution(f, M)
        p_values, mu = fitLinKK(f, ts, M, Z)

    return M, mu, eval_linKK(p_values, R0, ts, f), \
        residuals_linKK(p_values, R0, ts, Z, f, residuals='real'), \
        residuals_linKK(p_values, R0, ts, Z, f, residuals='imag')

((k-1)/(M-1))*np.log10(t_max/t_min))

        #ts /= 2*np.pi

        return ts

    if c is not None:
        M = 0
        mu = 1
        while mu > c and M <= max_M:
            M += 1
            ts = get_tc_distribution(f, M)
            p_values, mu = fitLinKK(f, ts, M, Z)

            if M % 10 == 0:
                print(M, mu, rmse(eval_linKK(p_values, ts, f), Z))
    else:
        M = max_M
        ts = get_tc_distribution(f, M)
        p_values, mu = fitLinKK(f, ts, M, Z)

    return M, mu, eval_linKK(p_values, ts, f), \
        residuals_linKK(p_values, ts, Z, f, residuals='real'), \
        residuals_linKK(p_values, ts, Z, f, residuals='imag')