function fxroot

% -- Get possible interval with the root
xb = brackPlot('myfx',0,5);

% -- Use bisection to refine the interval
r = bisect('myfx',xb);

% -- Add the value returned by bisect for visual confirmation
hold('on')
fr = myfx(r);
plot(r,fr,'*')
hold('off')

fprintf('\nRoot found at x = %8.5f  where f(x) = %12.3e\n',r,fr)

end