質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.50%
Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

2回答

1649閲覧

PYTHON 計算 ニュートン法

student_jojo

総合スコア9

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2018/07/10 05:04

前提・実現したいこと

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/ツールのバージョンなど)

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答2

0

ベストアンサー

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

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

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

投稿2018/07/10 05:25

mather

総合スコア6753

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

student_jojo

2018/07/10 06:04

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

2018/07/11 00:13

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

0

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

投稿2018/07/10 05:17

WathMorks

総合スコア1582

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問