%xsin(x) =1;
clear all;
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.*f(b)-b*f(a))./(f(b)-f(a));
    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