xl=3.0;
xu=4.0;
xr=(xl+xu)/2.0;
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=(xl+xu)/2.0 end
xl=3.0;
xu=4.0;
xr=(xl+xu)/2.0;
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=(xl+xu)/2.0 end