%xsin(x) =1
clear all;
hold off;
clc;
f = inline('x*sin(x)-1');
fplot(f,[0,2]);
hold on;
grid on;
a=0;
b=2;
plot(a,0,'.g');
plot(b,0,'.r');
for i =1:20;
    c=(a+b)/2;
    if f(a).*f(c)>0
        a=c;
    else
        b=c;
    end
disp([a,b,c]);
plot(a,0,'.g');
plot(b,0,'.r');
end