f(x) = x^2 - m,
の根 √m をニュートン法で求める。
x=a。から始め、
 a_{k+1} = a_k - f(a_k)/f '(a_k) = (1/2)(a_k + m/a_k),
に従って進む。
 a_k = (√m) coth((2^k)θ),
ただし cothθ = a。/(√m),