let no_root_out = function
  p -> let
         bound = div_rat (sum_of_coef_abs_val (p), abs_rat (leading_coef (p)))
       in
         (* bound > 0 then the interval is non-empty *)
         cons_open (opp_rat (bound), bound)