Неберущийся интеграл, как бы аппроксимировать его

так этот интеграл и не возьмётся, ты сам то подставь туда -inf, и получишь +inf.

Можно метод Симпсона попробовать. Но мне кажется, что интеграл расходится.

Ну ты график построй и посмотри, функция вполне себе ограничена.

достаточно бредовый вариант: взять математику, которая норм считает, и нагенерить таблицу значений для твоего диапазона параметров.

Функция выглядит нормально, думаю, что любое базовое правило подойдет. Хотя расчёт экспонент может занять какое-то время.

Переписал бы ты его в другом виде для наглядности, например, так:

2^(-e^(-x/a)) e^(-x/b), b>a>0

с точностью около процента?

Интегрирование методом Монте-Карло © «хватит всем».

Монте-карло от -inf до +inf? Не, тут надо какого-нибудь чувака со стат-физики, это они умеют прикидывать интегралы по фотографии.

А ты уверен, что он сходится?

2^(-e^(-x/a)) e^(-x/b), b>a>0

Доколе? Когда на ЛОРе появятся нормальные формулы вместо этого месива? %)

Монте-карло от -inf до +inf?

Можно заменой переменных перейти от интеграла с бесконечными пределами к интегралу с конечными пределами.

P.S. Ещё можно глянуть Maple/method: _d01amc — метод для вычисления интегралов на бесконечных интервалах ©.

Если ты просто хочешь посчитать без формального доказательства, возми численный от -2^n до 2^n, а потом несколько значаний в 2^m, m ∈ (n.. 100n) и прикинь значение. По виду, у тебя он будет сходиться очень быстро.

Вот когда появятся из коробки, тогда и буду на них писать. А так можешь эту строчку в вольфрам загнать, это нормальная для него формула.

да, тупанул.

t = e^(-x/a)

t'= e^(-x/a)*(-1/a)=-t/a

2^(-t)*(t^(a/b))*t'/t' dx =

-1*2^(-t)*(t^(a/b-1))/a dt

n = a/b-1 в диапазоне (-1,0)

t в диапазоне {0,+inf}

$e^(-log(2)*t)*t^n dt = Г(n+1)/(log(2)^(n+1)) = Г(a/b)/(log(2)^(a/b))

по формуле на второй странице седьмая снизу

https://radfiz.org.ua/files/NOT%20SORTED/zemskoff_IVT_s6_20140125_yse/%EB%D7%...

вот-с.

Спасибо. Точно, гамма получается, только a·log(2)^(-a/b) · Γ(a/b).