Apply bisection method (#585)

Apply bisection method (#585)

Repeatedly halve the interval containing the root and select the subinterval where the function changes sign.