Clean up code a bit.
Inline EPSILON and TOLERANCE for better performance in Numerical.
Share code that handles bounds in roots between solveQuadratic() and solveCubic()