newton.m

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

function [x,n]=newton(x0,e,et,N)
%x=newton(e,et,n)
n=0;
while fdiff(x0)~=0|n<=N
    x1=x0-f(x0)/fdiff(x0);
    if abs(x1-x0)<e|abs(f(x1))<et
        x=x1;
        return
    end
    x0=x1;
    n=n+1;
end