two.m

非线性方程求根 包括二分法、牛顿迭代法的matlab程序。

function [x,n]=two(a,b,e)
%x=two(f,a,b,e)
if f(a)*f(b)>=0
    return
end
n=0;
while abs(a-b)>e
    x=(a+b)/2;
    y=f(x);
    if abs(y)<e
        return
    end
    if f(a)*f(x)<0
        b=x;
    end
    if f(x)*f(b)<0
        a=x;
    end
    n=n+1;
end
n