f=@(x) sin(pi*x)

%-1*(x<0)+ 1*(x>=0)
%x.^3-5*x ;
for m=1:20
a=0.5;
b=2*m+1.5;
N=30;

for n=1:N
c=(a+b)/2 ;
if f(c)*f(a) >0
    a=c;
else
    b=c;
end
end

c=(a+b)/2

end