Robust ray intersection with interval arithmetic

For more general types of equations f (t) = 0, the problem of finding roots can be divided into two steps. First the roots must be isolated by finding intervals [ti , ti+1] which are known to contain one and only one root of the function. Secondly, the intervals about each root can be refined by reducing the size of the isolating interval until the root is located as accurately as possible using machine arithmetic.