Bisection Method
Matlab
Code:
LEFT_LIMIT = input('GIVE LEFT
LIMIT :');
RIGHT_LIMIT = input('GIVE RIGHT
LIMIT :');
MIDDLE_VALUE=0.0;
DIFFERECE_LIMIT= 0.0;
DOA=0.001;
x=0.0;
ymiddle= 0.0;
yleft=0.0;
yright=0.0;
counter=0.0;
a = LEFT_LIMIT:.8:RIGHT_LIMIT;
b = a.^2-7*a+10;
plot(a,b,'-rs'),xlabel('x axis'),ylabel('computed
function value'),legend('f(x)'),grid;
title('MATLAB
SESSIONAL');
for counter=1:1:1200
x = LEFT_LIMIT;
yleft = x^2-7*x+10;
x = RIGHT_LIMIT;
yright = x^2-7*x+10;
if (yleft*yright)>0
fprintf('INTERVAL IS NOT CORRECT');
break;
end
DIFFERECE_LIMIT =
abs(RIGHT_LIMIT-LEFT_LIMIT );
if DIFFERECE_LIMIT<DOA
fprintf('APPROXIMATE ANSWER : ');
disp(RIGHT_LIMIT);
break;
else
MIDDLE_VALUE = (RIGHT_LIMIT+LEFT_LIMIT)/2;
x = MIDDLE_VALUE;
ymiddle = x^2-7*x+10;
if (ymiddle*yleft)>0
LEFT_LIMIT = MIDDLE_VALUE;
else
RIGHT_LIMIT = MIDDLE_VALUE;
end
end
end