xxxx

a = 1.0;
b = 2.0;
epos = 0.00010;

f2=@(x) exp(x)/x-3;
f=@(x) x*x*x-x*x-3;
xk1 = 1;

for j=0:20
    xk=(a+b)/2;

    if (abs(xk-xk1)<0.00010)
        break;
    elseif ((f(a)*f(xk)) < 0)
        b = xk;
    else
        a = xk;
    end
    xk1=xk;
end
    xk
    j