matlab二分法求根

s1=1;s2=2;
answer=0;
f=inline('x^3+4*x^2-10','x')
%别忘记计算需要迭代的次数n
%这个n加一减一是啥意思
n=fix((log(s2-s1)-log(eps))/log(2)-1)+1;
%fix()函数是取整数

for j=1:n
    if((f(s1+s2)/2)*f(s2)<0)
        s1=(s1+s2)/2;
        answer=s1;
    elseif ((f(s1+s2)/2)*f(s2)>0)
        s2=(s1+s2)/2;
        answer=s2;
    else
        if(f(s1)==0)
            answer=s1;
            break;
        elseif(f(s2)==0)
            answer=s2;
            break;
        end
    end
end

answer
n