PYTHON 計算 ニュートン法

Question

前提・実現したいこと

PYTHONを使用して方程式を解いています。
ニュートン法で解きたいです。
ブログ等を参考にし、自分なりにコードを書いてみたのですが
以下のエラーメッセージが発生しました。

発生している問題・エラーメッセージ

Traceback (most recent call last):
  File "C:\Users\jojo2\Desktop\�R�clab\PYTHON\�s�X�g����~�ʒu�[�������v�Z.py", line 25, in <module>
    if abs(x2 - x) < 0.1:
  File "C:\Users\jojo2\AppData\Local\Programs\Python\Python36\lib\site-packages\sympy\core\relational.py", line 195, in __nonzero__
    raise TypeError("cannot determine truth value of Relational")
TypeError: cannot determine truth value of Relational

該当のソースコード

from sympy import *
x = Symbol('x')

p10 = 500000
p20 = 60000
v10 = 0.00444
v20 = 0.00188
r1 = 1.4
r2 = 1.67
s = 0.00096

g = (p10/p20) * (v10/v20) * ((r2 - 1)/(r1 - 1)) * (((1 + (s*x / v10)) ** (1 - r1)) - 1) + (((1-(s*x / v20))** (1-r2)) -1)
f = diff(g)
print (f)

#　微分完了
#　ここからニュートン法

x = 1.8

while True:

    x2 =  x - g / f

    if abs(x2 - x) < 0.1:
        break

    x = x2

print(x)

試したこと

誤差の指定が悪かったのかと、誤差を変更してみたのですが同様に失敗し、
なぜエラーが発生するのか色々と検索をかけてみましたが
エラー発生原因も理解できておりません...

補足情報（FW/ツールのバージョンなど）

Answer 1

f, g は x に関する式ですので、 g / f も x に関する式です。 x の値は代入されていません。
そのため、そもそも絶対値などの計算ができない状況です。

f, g に x の値を代入して計算するようにしましょう。

なお、 x = Symbol('x') と宣言しているのですから x = 1.8 と上書きしてしまうのはやめましょう。
x0, x1 など別の名前を使うほうが良いです。

Comments

[mather]

参考： https://note.nkmk.me/python-sympy-factorization-solve-equation/

[student_jojo]

mather様
ご回答、参照のご提示ありがとうございます。
解法を掴めたような気がしますので、頑張って完成させてみます。

[student_jojo]

解決しました！
本当にありがとうございました！
次はグラフ化にチャレンジしようと思います！

Answer 2

f,gに現在の値xを代入する.