python ニュートン法による√xの近似値を求める関数の定義

Question

前提・実現したいこと

ニュートン法を用いて、x の平方根√x の値（近似値）を求める関数を定義したいです。
以下、与えられた条件です。
・初項　a1 = x
・漸化式　an+1 = an - (an^2 - x)/2an

なお、今回はa10 の値を求めるコードを作成したいです。

自分でコードを作成し、実行したところ
計算結果が nan となりました。
どう修正すればよいか教えていただきたいです。

該当のソースコード

def sqrt(x):
    a = x
    n = 1

    while n < 10:
      x = x - (x**2 - a) / 2*x
      n += 1
    return x

Answer 1

正解を書くのもどうかと思うので、以下の指摘にとどめます。

• 漸化式　an+1 = an - (an^2 - x)/2an

と

• x = x - (x**2 - a) / 2*x

を良く見比べてください。

なお、Pythonでは他の多くの言語と同じように、

A / B * C は、(A / B) * C　ですので注意してください。

Comments

[noob_mugy]

ご回答ありがとうございます。

ご指摘通りに修正したところ、正しく計算できました。
ありがとうございました。