Вычислим следующий несобственный интеграл:

\int_0^\infty e^{x^2}\mathrm{d} x

Данный интеграл называется интеграл Гаусса. Он играет очень большую роль в математике, особенно в теории вероятности.

Я тут приведу классический способ его вычисления. Его можно найти во многих учебниках по дифференциальному и интегральному исчислению в разделах кратных интегралах. Я, например, его взял из учебника Н. С. Пискунова “Дифференциального и интегрального исчисления для втузов” Москва издательство “Наука” 1985г. Но проанализировав это решение, мы выведем необычную аппроксимацию интеграла:

\int_{-a}^a e^{x^2}\mathrm{d} x

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

Для этого умножим исходный несобственный интеграл на интеграл:

\int_{-\infty}^{\infty}{e^{-y^2}\mathrm{d}y}

Мы получили следующий двойной несобственный интеграл

\int_{-\infty}^{\infty}{e^{-y^2}\mathrm{d}y}\cdot\int_{-\infty}^{\infty}{e^{-x^2}\mathrm{d}x}

Определенный интеграл от одной и той же функции с одинаковыми верхними и нижними пределами всегда равны, независимо какой буквой мы обозначим переменную. Следовательно, данное произведение будет равно квадрату исходного интеграла.

Внеся подынтегральное выражения левого интеграла под правый, получим следующий двойной интеграл:

\int_{-\infty}^\infty\left(\int_{-\infty}^\infty \mathrm{e}^{x^2}\cdot\mathrm{e^{x^2}}\mathrm{d}x\right)\mathrm{d} y=\int_{-\infty}^\infty\left(\int_{-\infty}^\infty\mathrm{e}^{-x^2-y^2}\mathrm{d}x\right)\mathrm{d}y

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

\int_{-\infty}^{\infty}\left(\int_{-\infty}^{\infty}{e^{-x^2-y^2}\mathrm{d}x}\right)\mathrm{d}y=\int_{0}^{2\cdot\pi}\left(\int_{0}^{\infty}{e^{-\rho^2}\cdot\rho\cdot \mathrm{d}\rho}\right)\mathrm{d}\theta

Так как внутренний интеграл не зависит от θ, то их можно представить как произведение двух независимых интегралов. Внутренний несобственный интеграл представим как произведением двух независимых интегралов. Внутренний несобственный интеграл представим как предел при R стремящейся к бесконечности от интеграла с пределами от нуля до R:   

\int_{0}^{2\cdot\pi}\left(\int_{0}^{\infty}{e^{-\rho^2}\cdot\rho d\rho}\right)d\theta=\int_{0}^{2\pi}d\theta\cdot\lim_{R\rightarrow\infty}{\int_{0}^{R}{e^{-\rho^2}p\mathrm{d}\rho}}=2\cdot\pi\cdot\lim_{R\rightarrow\infty}{\frac{1-e^{-R^2}}{2}}=\pi

Так как выше было показано, что данный интеграл является квадратом исходного интеграла, то имеем:

I=\int_{-\infty}^\infty\mathrm{e}^{x^2}\mathrm{d}x=\sqrt{\pi}

Мы нашли значение интеграла Пуассона. Но давайте посмотрим на выражение, которое мы получили в пределе. Оно монотонно возрастает с R и имеет определенное значение при R=∞. При этом мы показали, что данный придел является квадратом исходного интеграла. Следовательно:

\lim_{a\rightarrow\infty}{\int_{-a}^{a}e^{-x^2}dx}=\sqrt{\pi\cdot\left(1-e^{-2\cdot a^2}\right)}

Исследуем следующее выражение:

\left|\int_{-a}^{a}e^{-x^2}dx-\sqrt{\pi\cdot\left(1-e^{-2\cdot a^2}\right)}\right|

Фактически данное выражение является абсолютной ошибки предполагаемой предложенной формулы аппроксимации и точного значения интеграла Гаусса при разных пределах

Ниже приводится график абсолютной ошибки, причем ось ординат имеет логарифмический масштаб:

График ошибки нашей аппроксимации
График ошибки нашей аппроксимации

Из графика видно, что максимальное значение данной функции в районе 0.7 и составляет примерно 0.200. Ниже приводится таблица точного интеграла, интеграла, вычисленного по предложенной формуле для аппроксимации и ее ошибки, при разных пределах интегрирования.

аргумент

Значение интеграла

Значение интеграла, вычисленное по аппроксимирующей формуле

абсолютная погрешность

0.7

1.201371069

1.400899019

0.199527950

1.0

1.493648146

1.648157857

0.154509708

1.5

1.712376674

1.762581264

0.050204590

2.0

1.764162592

1.772156530

0.007993947

2.5

1.771732394

1.772450582

0.000718154

3.0

1.772414425

1.772453837

0.000039118

На основание этих данных можно оценить интеграл:

\int_{-a}^a\mathrm{e}^{-x^2}

с помощью выражения:

\sqrt{\pi\cdot\left(1-\mathrm{e}^{-a^2}\right)}

При a>2 с достаточно высокой точностью можно вычислить данный интеграл.

Предложенная формула не является строгим следствием приведённого выше рассуждения, однако численный анализ показывает её высокую точность в указанном диапазоне. Для точных вычислений рекомендуется использовать функцию ошибок, но данная аппроксимация удобна для быстрых инженерных оценок.