trootnewt.m

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

function trootnewt(x0,eps)
k=1
xk=4*x0^2-14*x0+14
while abs(xk-x0)>eps
    k=k+1
    p=xk;
    xk=4*xk^2-14*xk+14
    x0=p;
end
r=xk