Question

matlab实现二分法

function [x]=erfen(f,a,b,eps)
c=(a+b)/2;
fa=feval(f,a);
fb=feval(f,b);
fc=feval(f,c);
while abs(b-a)>eps;
if fb*fa==0;
return;
else if fb*fc<0
a=c;
else b=c;
end;
end;
end;
x=c;

麻烦大神看眼哪里出问题了

Answer 1

c=(a+b)/2;fa=feval(f,a);fb=feval(f,b);fc=feval(f,c); 这些应该在循环内 不然你循环里面都是搞空事啊