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

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

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

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

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

Q&A

解決済

1回答

5561閲覧

sympyで積分した結果が数式であるものを数値にしたい

yoko28

総合スコア14

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

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

0グッド

0クリップ

投稿2020/08/19 04:56

編集2020/08/19 04:57

Pythonは、初級者です。
質問の文も拙いと思いますがよろしくお願いいたします。

前提・実現したいこと

  • sympyで積分した結果が数式になっているものを数値にしたいです。

ここでは、( sin(x) × tanh(x) )を[-20,21]の範囲で定積分を行います

発生している疑問、試したこと

  • 上記の式をsympyで積分した結果が以下のような数式として出力されます。

(下記のソースコードのprintの出力結果)

Integral(sin(r1)*tanh(r1), (r1, -20, 21)) <class 'sympy.integrals.integrals.Integral'>
  • また、積分の計算サイトで同じ式で計算した結果、1.50479という結果が得られました(現在の設定であれば数値化できる)

計算サイトの結果LINK

  • 計算サイトの結果のように積分の結果を数値にするには、何かに変換する必要があるのでしょうか。

該当のソースコード

python

1 r1 = sympy.Symbol('r1') 2 x1 = -20 3 x2 = 21 4 res = sympy.integrate( sympy.tanh(r1) * sympy.sin(r1) , (r1, x1,x2 )) # <-積分の計算部分 5 print(res) 6 print(type(res))

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

python 3.6

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

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

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

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

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

guest

回答1

0

ベストアンサー

解析解が計算できないので、res は積分式のままになっているのだと思います。
evalf() で評価すると、数値解になります。

python

1import sympy 2 3r1 = sympy.Symbol("r1") 4x1 = -20 5x2 = 21 6res = sympy.integrate(sympy.tanh(r1) * sympy.sin(r1), (r1, x1, x2)) 7print(res.evalf()) # 1.50478609907259

integrate sin(x)*tanh(x), -20 <= x <= 21 - Wolfram|Alpha

投稿2020/08/19 05:36

編集2020/08/19 05:38
tiitoi

総合スコア21956

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

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

yoko28

2020/08/19 05:49

回答ありがとうございます。evalfというのがあるのですね
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.46%

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

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

質問する

関連した質問