---

def max_diff(xc, yc, xf, yf):
    m = 0.0
    for xi, yi in zip(xc, yc):
        idx = np.argmin(np.abs(xf - xi))
        m = max(m, abs(yi - yf[idx]))
    return m

eps1 = max_diff(x1c, y1c, x1f, y1f)
eps2 = max_diff(x2c, y2c, x2f, y2f)

# Дополнительно: расхождения между системами
diff_coarse = max(abs(ya - yb) for ya, yb in zip(y1c, y2c))
diff_fine   = max(abs(ya - yb) for ya, yb in zip(y1f, y2f))