xl = 3.0; xu = 8.0; eps = 0.001;  maxit = 50;
root1 = bisection( @myfunc3 , xl, xu, eps, maxit );
xl = 1.0; xu = 3.0; eps = 0.001;  maxit = 50;
root2 = bisection( @myfunc3 , xl, xu, eps, maxit );

fprintf('*** root1 = %8.2f\n', root1 );
fprintf('*** root2 = %8.2f\n', root2 );