二分法解方程MATLAB

clear
a=input('a=');
b=input('b=');
e=input('c=');
c=(a+b)/2
n=log2((b-a)/e)-1
if f(a)*f(b)>0
    Solution='May be No Solution'
else
    if f(c)==0
        Solution=c
    else
        for i=1:1:n
            c=(a+b)/2
            if f(c)<0
                a=c
            else
                b=c
            end
        end
        Solution=c
    end
end