global f_calls;
global bisect_num;

f_calls = 0;

xtol = input ( 'Enter interval tolerance XTOL: ' );
a = 0.0;
b = 1.0;
[ x, a, b ] = bisection5 ( xtol, a, b, @cosxx2 );

fprintf ( 'Estimated root is %20.16g\n', x );
fprintf ( 'Number of f calls =    %d\n', f_calls );
fprintf ( 'Number of bisections = %d\n', bisect_num );
