trootxian.m

来自「数值分析各种公式的解法」· M 代码 · 共 14 行

 function trootxian(x0,x1,eps)
k=1
xk=x1-f1(x1)*(x1-x0)/(f1(x1)-f1(x0))
x0=x1;
while abs(xk-x0)>eps
    k=k+1
    q=xk;
    xk=xk-f1(xk)*(xk-x0)/(f1(xk)-f1(x0))
    x0=q;
end
r=xk
function f1=f1(x)
f1=x*exp(x)-1;
%f=0