MATLAB实现二分法

syms x;

f=input('(输入函数形式)fx=');

a=input('（输入二分法下限)a=');

b=input('（输入二分法上限)b=');

tic;

c=(a+b)/2;

k=1;

c=double(c);

fprintf('k=%d x=%f\n',k,c);

while abs(c-a)>0.0001&&k<30;

if subs(f,x,a)*subs(f,x,c)<0;

b=c;

elseif subs(f,x,c)*subs(f,x,b)<0; a=c;

else

a=c;b=c;

end

c=((a+b)/2);

c=double(c);

k=k+1;

fprintf('k=%d x=%f\n',k,c);

end;

toc;