08/05/06 13:56:32
ニュートン法を使った平方根計算を作りたいのですが、例えば
Nの平方根を求める。初期値 a0、許容誤差 epsとして
def my_sqrt(n, a0, eps):
x = a0
new_x = 0.5 * ( x + n/x)
if abs(new_x -x) < eps:
# print 'In if abs statement:', new_x
return new_x
else:
my_sqrt(n, new_x, eps)
>>> my_sqrt(5.0, 2.0, 1.e-5)
# 5の平方根を初期値2から求める。誤差1.e-5になったら終了。
としても、returnで値が帰ってこないんですがなぜでしょうか？
#でコメントアウトしたprint文は(コメントを取ると)出力されるので、
再帰計算中に、許容誤差範囲に収束しないという事ではなさそうです。