xianjie.m

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

function [x,n]=Xianjie(f,x0,x1,eps)
% xianjie迭代法
if nargin ==3
    eps=1e-5;
end
n=0;
x2=x1-f(x1)*(x1-x0)/(f(x1)-f(x0));
while abs(x1-x0)>eps && abs(f(x1))>eps
    x0=x1;
    x1=x2;
    x2=x1-f(x1)*(x1-x0)/(f(x1)-f(x0));
    n=n+1;
end
x=x2;
end