Bisection method using mathematica