Bisection method in matlab

HI I wanna graph the bisection method with the function that I have but Idk how to do it. I also want to Iterate until the relative approximate error falls below 0.01% or the number of iterations exceeds 100. this is what I have so far but for some

this is the code

clc

clear

lc=3; lp=3; w=160;

T= 700;

f=@(d) (w*lc*lp/(d*sqrt(lp^2-d^2)))-T;

xl=input('enter the value for Xl');

xu=input('enter the value for Xu');

xm = (xu + xl) / 2;

if f(xl)*f(xu)>0

xmnew=xl

xunew=xu

%calculate the error

error=abs((xmnew-xm)/xmnew)

else if f(xl)*f(xu)<0

xmnew=xu

xlnew=xl

error=abs((xmnew-xm)/xmnew);

else

fprintf('xm is the root')

end

end

[Asad (Mehrzad) Khoddam]

There is no loop for the iteration

[Dai Nguyen]

I don't really know how to create the loop for the iteration, can you help me?

[Dai Nguyen]

is it for i=1:xm in this case

Answer 1

[Asad (Mehrzad) Khoddam]

clc
clear
lc=3; lp=3; w=160;
T= 700;
f=@(d) (w*lc*lp/(d*sqrt(lp^2-d^2)))-T;
xl=input('enter the value for Xl ');
xu=input('enter the value for Xu ');
for i=1:100
    xm = (xu + xl) / 2;
    error=abs(xm-xu)/xm;
    if error< 0.01
        break;
    end
    if f(xl)*f(xm)>0
        xl = xm;
    elseif f(xl)*f(xu)<0
        xu = xm;
    else
        fprintf('%f is the root\n', xm)
        break;
    end
    
end