NA-lab-bisection

code

clc;

f=@(x) x^2-29;

a=input("Enter lower bound");

b=input("Enter upper bound");

e=input("Enter tolerance");

if(f(a)*f(b)<0)

while(abs(a-b)>e)

c=(a+b)/2;

if(f(a)*f(c)<0)

b=c;

elseif(f(b)*f(c)<0)

a=c;

else

break;

end

fprintf("The root is %0.4f\n",c);

end

else

fprintf("Not possible sign doesnt change in given interval");

end