newton.m

这是非线性方程（组）的数值解法

function [x,n]=Newton(f,x0,eps)
% Newton迭代法
if nargin ==2
    eps=1e-5;
end
g=inline(diff(sym(f))); % f 的导函数
x1=x0-f(x0)/g(x0);
n=1;
while abs(x1-x0)>eps && abs(f(x1))>eps
    x0=x1;
    x1=x0-f(x0)/g(x0);
    n=n+1;
end
x=x1;
end