I am in need of a good interval bisection routine to help me find the roots of ''difficult'' equations. I am following numerial recipes text in which they seem to do it in two steps. FIrst determine a number of intervals where the curve crosses x axis and then using an interval bisection program refine and output the root within the interval.

Is there anything already written in Mathematica to do this? I am having trouble with Newtons Method sometimes.
Are there any texts or articles on the subject that you can refer me to?