powered by
ecd.uniroot(f, lower, upper, use.mpfr = FALSE, tol = .Machine$double.eps^0.25, maxiter = 1000)