xl=0;
xu=2;
tolerance=0.0001;
relative_error=3;
xr=(xl+xu)/2
while(relative_error>tolerance)
fxl=xl^3+4*xl^2-10;
fxr=xr^3+4*xr^2-10;
if(fxl*fxr<0)
       xu=xr;
else
       xl=xr;
end
xrold=xr;
xr=(xl+xu)/2
relative_error=abs(xr-xrold)/xr
end
root=xr