erfenfa.m

来自「这是非线性方程（组）的数值解法」· M 代码 · 共 27 行

function [x,n,w]=Erfenfa(f,a,b,eps)
% x为根，n为迭代次数，w为先验绝对误差界
if sign(f(a))*sign(f(b))>0
    disp('不满足二分法条件！')
    return;
end
if nargin==3
    eps=1e-5;
end
a1=a;
b1=b;
n=1;
while abs((a-b)/2)>eps
    x=(a+b)/2;
    if abs(f(x))<eps
        return;
    end
    if sign(f(a))*sign(f(x))<0
        b=(a+b)/2;
    else
        a=(a+b)/2;
    end
    n=n+1;
end
x=(a+b)/2;
w=(b1-a1)/(2^n);
end