trootnewt1.m

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

function trootnewt1(x0,eps)
k=1
xk=x0-(2*exp(-x0)-sin(x0))/((-2)*exp(-x0)-cos(x0))
while (xk-x0)>eps
    k=k+1
    p=xk
    xk=xk-(2*exp(-xk)-sin(xk))/((-2)*exp(-xk)-cos(xk))
    x0=p
end
r=xk