Как можно численно интегрировать Exp [-Abs [x1] -Abs [x2]] / Abs [x1-x2] для x1 и x2 из -Infinity в Infinity?
Какую стратегию можно адаптировать, Монте-Карло? Алгоритм Вегаса (я пробовал реализации в Python и Mathematica) может работать, но я не могу найти пакет, который позволяет мне исключить набор точек, что делает NIntegrate Mathematica (здесь: Exclusion-> x1 == x2). Последний, однако, дает разные ответы в зависимости от стратегии (без сообщений об ошибках и с интегральными ошибками, доходящими только до 1,1, используя https://mathematica.stackexchange.com/questions/75426/obtaining-an-nintegrate-error-estimate). Я также go касаюсь пространства, используя x-> Tan [alpha], чтобы сжать пределы интегрирования до + -pi / 2.
Спасибо.
Плотность График функции вдоль x1 и x2 (в касательном пространстве)