bisection.m

来自「Numerical analysis （Gauss、Bisection、jaco」· M 代码 · 共 19 行

function s=Bisection(fx,xa,xb,n)
%bisection does n bisections to approximation
%a root of fx
x=xa;
fa=eval(fx);
x=xb;
fb=eval(fx);
for i=1:n
    xc=(xa+xb)/2;
    x=xc;
    fc=eval(fx);
    x=[i,xa,xb,xc,fc];
    disp(x);
    if fc*fa<0
        xb=xc;
    else 
        xa=xc;
    end % of if/else
end % of for loop