Вычислял параметрический интеграл


с помощью квадратурной формулы Лагерра. Фиксировал параметр p изменялся от 0.1 до 10.
Для отрезка [0.1 , 1.0) параметр изменялся с шагом 0.1, для [1 , 10] с шагом 0.5.
Чтобы вычисления были в общем виде для всех пределов сделал замену:

и затем

. Интеграл привел к новому виду:

Записал алгоритм на языке С++. Построил график

Рисунок из ответов:

Видно, что функция должна иметь максимум, мое решение этот максимум не дает.
Задача 420 (с.119) из учебного пособия:
А. И. Алексеев Сборник задач по классической электродинамике, Москва, 1977.
Ответ к задаче на с. 262