Question

Matlab 入门级问题

function [c,err,yc]=bisect(f,a,b,delta)
ya=feval(f,a);
yb=feval(f,b);
if ya*yb>0,return,end
max1=1+round((log(b-a)-log(delta))./log(2));
for k=1:max1
c=(a+b)/2;
yc=feval(f,c);
if yc==0
a=c;
b=c;
elseif yb*yc>0
b=c;
yb=yc;
else
a=c;
ya=yc;
end
if b-a<delta,return,end
end
c=(a+b)/2;
err=abs(b-a);
yc-feval(f,c);
一个二分法的程序，几乎是跟着书本打的，但输入
[a,b,c]=bisect('f',0,2,.0001)
时出现
??? Output argument "err" (and maybe others) not assigned during call to "E:\mathlab\work\bisect.m (bisect)".

Error in ==> bisect at 2
ya=feval(f,a);

求破。
function y=f(x)
y=x.*sin(x)-1;

Answer 1

feval(f,a);
这里面这个f是什么函数，把代码也贴上来吧。

ok，搞定了：
原因是这行代码：
if b-a<delta,return,end
如果进入这个条件语句，那么return，之后err就得不到了，但是你还要输出来，即输出一个并不存在的数。
你把return改为break就行了。
我已经测试通过了。