xianjie.m

来自「这是非线性方程（组）的数值解法」· M 代码 · 共 15 行

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