В одной
статье (формула

) есть такой ряд:


Также даются соотношения для вычисления внутренних сумм, в результате использования которых у меня получается:


где

-
символ Pochhammer'а, а

- полиномы, выражения для которых приведены в приложении статьи.
Однако, остаётся непонятным, как всё-таки далее суммировать это выражение по

. Я пробовал скормить исходное выражение Wolfram Mathematica, в общем случае не дождался. Когда подставляю конкретные значения

, после упрощения (без которого иногда получаются выражения с гипергеометрическими функциями) получаются числа вроде

но главная проблема - даже конкретные значения вычисляются
очень медленно (несколько минут на одно значение, а надо их тысячи).
Пробовал также вычислить численно, с помощью NSum - не дождался. Пробовал вынести конечные суммы за знак бесконечного суммирования - оказалось, что такая сумма, выражающаяся через гипергеометрическую функцию

, требует вычисления с огромной точностью этой функции (ибо коэффициенты внешних сумм оказываются очень большими), что, опять-таки, слишком медленно - не дождался.
С другой стороны, авторы говорят (стр. 20):
Цитата:
when

raises, the formula for calculating the coefficients

becomes much complicated
(здесь

- максимальное значение для

), что явно указывает на то, что у них есть метод получения аналитических решений. Я пробовал написать авторам, но электронный адрес, указанный в статье, уже невалиден (письмо вернулось).
Поэтому вопрос: как вообще подходить к аналитическому суммированию этого ряда?