trootait.m

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

function r=trootait(x0,eps)
k=1
xk=ff(x0)
xk1=ff(xk)
x=xk1-(xk1-xk)^2/(xk1-2*xk+x0)
while abs(x-x0)>eps
    k=k+1
    p=x;
    xk=ff(x)
    xk1=ff(xk)
    x=xk1-(xk1-xk)^2/(xk1-2*xk+x)
    x0=p;
end
r=x