🎄teratailクリスマスプレゼントキャンペーン2024🎄』開催中!

\teratail特別グッズやAmazonギフトカード最大2,000円分が当たる!/

詳細はこちら
Python

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

Q&A

1回答

4304閲覧

Pythonで重積分を行いたい

Lily1007

総合スコア10

Python

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

0グッド

1クリップ

投稿2020/12/23 11:05

編集2020/12/23 11:23

Pythonで重積分を行いたい

写真のような重積分を行いたいです。
イメージ説明
(しかし,dx'の積分は∞までではなく、+5まで。C, l'はそれぞれ値1を代入しているとします。)
tau_0 = 10**(np.linspace(-1.5, 5., 100))として,各τ(0)に帯するW'の積分値を出してそれをplotするのが最終目標です。

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

写真の積分が出来ません。まずdtの積分の中にもx'が含まれており,dtの積分後にx'の積分をするからです。ネットで出てきた重積分を行う方法では,被積分関数が1つのみ(例:x^2 + y^2のx,yでの積分)のものしかでませんでした。

この重積分はどのようにすればよいのでしょうか。
よろしくお願いいたします。

補足情報

python

1x = np.linspace(0, 5, 100) 2t = np.linspace(-1, 1., 100) 3tau_0 = [1/16., 1/8., 1/4., 1/2., 1., 2., 4., 8] 4C = 1.0 5l = 1.0 6func1 = np.exp(-(1/(2*math.pi)*(t/1-t**2)-x)**2)/(1+t**2) 7func2 = integrate.simps(func1, t) 8func2 = 1-np.exp(-2*C*tau_0*func2)

ここまで書いて手が止まってしまいました。

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

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

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

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

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

guest

回答1

0

数値積分
の57ページ以降に書かれてるやり方で、どうでしょうか?

投稿2020/12/24 22:58

jbpb0

総合スコア7653

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

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

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだベストアンサーが選ばれていません

会員登録して回答してみよう

アカウントをお持ちの方は

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

ただいまの回答率
85.36%

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

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

質問する

関連した質問