Использование scipy.integrate.quad для выполнения 3D-интеграла ⇐ Python

[Anonymous]

Мотивация вопроса

Я пытаюсь проинтегрировать функцию f(x,y,z) по всему пространству.

Я пробовал использовать scipy.integrate.tplquad и scipy.integrate.nquad для интегрирования, но оба метода возвращают интеграл как 0 (когда интеграл должен быть конечным). Это связано с тем, что по мере увеличения объема интегрирования область, в которой подынтегральная функция не равна нулю, выбирается все меньше и меньше. Интеграл «пропускает» эту область пространства. Однако scipy.integrate.quad, похоже, способен справиться с интегралами от [-infinity, infinity], выполняя замену переменных...

Вопрос

Можно ли использовать scipy.integrate.quad 3 раза для выполнения тройного интеграла. Код, который я имею в виду, будет выглядеть примерно так:

Код: Выделить всё

x_integral = quad(f, -np.inf, np.inf)
y_integral = quad(x_integral, -np.inf, np.inf)
z_integral = quad(y_integral, -np.inf, np.inf)

where f is the function f(x, y, z), x_integral should integrate from x = [- infinity, infinity], y_integral should integrate from y = [- infinity, infinity], and z_integral should integrate from z = [- infinity, infinity]. I am aware that quad wants to return a float, and so does not like integrating a function f(x, y, z) over x to return a function of y and z (as the x_integral = ... line from the code above is attempting to do). Is there a way of implementing the code above?

Thanks


Источник: https://stackoverflow.com/questions/513 ... d-integral