double Newton(cut c, double eps) {
    double x0, x1 = (c.first + c.second) / 2.0;
    do {
        x0 = x1;
        x1 = x0 - f_(x0) / f__(x0);
    } while (abs(x1 - x0) > eps);
    return x1;
}