新建 文本文档 (3).txt

来自「matlab中GUI界面来实现运筹学中常用的几种最优算法」· 文本 代码 · 共 18 行

function x=bisect(fun,a,b,ess)
fa=feval(fun,a);
fb=feval(fun,b);
if fa*fb>0
    error('函数两端点的值应异号');
end
x=(a+b)/2
while (b-a)>ess
    fx=feval(fun,x);
    if fa*fx<0
        b=x;
        fb=fx;
    else
        a=x;
        fa=fx;
    end
    x=(a+b)/2
end