Partiendo de cero, suponemos que la DE $$ (x^2+2)y''+xy'-y=0 \tag{1}$$ tiene una solución analítica en una vecindad de cero, dada por $$ y(x) = \sum_{n\geq 0} a_n x^n \tag{2} $$ En cuanto a los coeficientes, $(1)$ se traduce en:
$$ (x^2+2)\sum_{n\geq 2} n(n-1)a_n x^{n-2} + x\sum_{n\geq 1}na_n x^{n-1} - \sum_{n\geq 0}a_n x^n = 0 \tag{3}$$ que equivale a: $$ \sum_{n\geq 2}n(n-1)a_n x^n +\sum_{n\geq 0}2(n+2)(n+1)a_{n+2}x^n + \sum_{n\geq 1} n a_n x^n -\sum_{n\geq 0}a_n x^n = 0 \tag{4} $$ Ahora consideramos cuál es el coeficiente de $x^n$ en la última expresión y obtenemos eso,
por cada $n\geq 2$ , $$ n(n-1) a_n + 2(n+2)(n+1) a_{n+2} + (n-1)a_n = 0\tag{5} $$ por lo tanto: $$ \frac{a_{n+2}}{a_n} = - \frac{(n-1)}{2(n+2)} \tag{6} $$ y por inducción: $$ a_{n+2k} = \frac{(-1)^k}{2^k}\cdot\frac{(n+2k-3)(n+2k-5)\cdot\ldots\cdot(n-1)}{(n+2k)(n+2k-2)\cdot\ldots\cdot(n+2)}\,a_n \tag{7}$$ por lo que una solución de $(1)$ viene dada por una función hipergeométrica.