Такой лайфхачек придумался. Если записать ряд Тейлора функции в виде
(где суммирование идёт по всем целым значениям), то не нужно будет беспокоиться об исключениях, если только предположить
(Тем более, что так оно и есть).
Проиллюстрирую идею на полиномах Лежандра.
Функция переходит в
(стрелка заменяет громоздкий знак суммы), её производные и прочие комбинации в
,
,
,
и в итоге получается следующие рекуррентные соотношения
Поскольку мы оставили факториалы в знаменателе, далее легко выписываются все исключения
Потом, конечно, всё это можно свернуть