sep.m

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

function [y,n]=sep(fun,a,b,ess)

if nargin<4
    ess=0.0001
end
if feval(fun,a)*feval(fun,b)<0
    count=1;
    c=(a+b)/2;
    while c>ess
        if feval(fun,a)*feval(fun,c)<0
            b=c;
            c=(a+b)/2;
        elseif feval(fun,c)*feval(fun,b)<0
            a=c;
            c=(a+b)/2;
        else
            y=c;
            ess=100000;
        end
        n=n+1;
    end
    y=c;
elseif feval(fun,a)==0
    y=a;
elseif feval(fun,b)==0
    y=b;
else disp('cannot find the superior resolution')
end