xl=3.0;
xu=6.0;
fxl=exp(xl)-10*xl;
fxu=exp(xu)-10*xu;
xr=xu-fxu*(xl-xu)/(fxl-fxu);
fxr=exp(xr)-10*xr;
while abs(fxr)>1.0e-6
fxl=exp(xl)-10*xl;
fxu=exp(xu)-10*xu;
fxr=exp(xr)-10*xr;
if fxr*fxl < 0 xu= xr; elseif fxr*fxu < 0 xl= xr; end xr=xu-fxu*(xl-xu)/(fxl-fxu) end