% test_bisect.m % Test the implementation of the bisection method in bisect.m g = @(x) cos(x) - x; a = 0; b = pi/2; r = bisect( g, a, b, 1.0e-19 )