trootsqrt.m

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

function trootsqrt(c,x0,eps)
k=1
xk=(x0+c/x0)/2
while abs(xk-x0)>eps
    k=k+1
    p=xk
    xk=(xk+c/xk)/2
    x0=p;
end
sqrtc=xk