Ridders' Method: f(x) = 0
Brief Description of the Method
Given the initial bracketing interval [T0, T2], containing one solution, and such that f(T0) and f(T2) have opposite signs, compute the midpoint T1 = (T0 + T2)/2. Then define the new function h(z) = f(z)/sqrt(f(z)^2 + f(T0)*f(T2)). Find the parameter that ensures h(T1) = (h(T0) + h(T2))/2. Then define T3 as the x-intercept of the line passing through (T0, h(T0)) and (T2, h(T2)). Use T3 as one of the endpoints of the next interval bracketing the solution. The other endpoint is T1 if f(T2) * f(T3) < 0. Otherwise, choose either T0 or T2 based on the requirement that the sign of f(z) at the chosen point must be opposite to the sign of f(T3). Continue until the desired accuracy is reached.