trootnewt.m

数值分析各种公式的解法

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