二分法解方程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