Такой лайфхачек придумался. Если записать ряд Тейлора функции в виде

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

(Тем более, что так оно и есть).
Проиллюстрирую идею на полиномах Лежандра.

Функция переходит в

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

,

,

,

и в итоге получается следующие рекуррентные соотношения
![$$\frac{1}{{k!}}u_{k + 2} = \left[ {\frac{1}{{(k - 2)!}} + \frac{2}{{(k - 1)!}} - \frac{{n(n + 1)}}{{k!}}} \right]u_k $$ $$\frac{1}{{k!}}u_{k + 2} = \left[ {\frac{1}{{(k - 2)!}} + \frac{2}{{(k - 1)!}} - \frac{{n(n + 1)}}{{k!}}} \right]u_k $$](https://dxdy-04.korotkov.co.uk/f/7/e/3/7e3487b387b7b292625e8ea252d69b2982.png)
Поскольку мы оставили факториалы в знаменателе, далее легко выписываются все исключения

Потом, конечно, всё это можно свернуть
![$$\[
u_{k + 2} = \left[ {k(k + 1) - n(n + 1)} \right]u_k \quad \forall k \in \mathbb{N}
\]
$$ $$\[
u_{k + 2} = \left[ {k(k + 1) - n(n + 1)} \right]u_k \quad \forall k \in \mathbb{N}
\]
$$](https://dxdy-03.korotkov.co.uk/f/6/2/6/62617ab2342cd1bef3245d69319ad11782.png)