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)
ここまで書いて手が止まってしまいました。
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。