erfen.m

来自「非线性方程组解法源代码」· M 代码 · 共 27 行

function erfen
%二分法解方程
syms x y ;
y = input('请输入迭代函数f(x)：');
v = input('请输入迭代条件[a,b,e]:');
a = v(1);
b = v(2);
e = v(3);
while abs(b-a)>e
    c = 0.5 * (a + b);
    x = c;
    fc = eval(y);
    x = a;
    fa = eval(y);
    if fc==0
        c
    elseif fa==0
        a
    elseif fa*fc>0
        a = c;
    else
        b = c;
    end
end
a
b
0.5 * (a + b)