trootsqrt.m

数值分析各种公式的解法

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