Мотивация: У меня есть многомерный интеграл, который для полноты я воспроизвел ниже. Он исходит из вычисления второго вириального коэффициента при существенной анизотропии:
Здесь W - функция всех переменных. Это известная функция, для которой я могу определить функцию python для.
Вопрос программирования. Как мне получить scipy
для интеграции этого выражения? Я думал о соединении двух тройных квадов (scipy.integrate.tplquad
), но я беспокоюсь о производительности и точности. Существует ли более высокий размерный интегратор в scipy
, который может обрабатывать произвольное число вложенных интегралов? Если нет, каков наилучший способ сделать это?