Problem: Given positive integer $n$, solve the recurrence relation $a_k = (n+k)a_{k-1} - ka_{k-2}, \ k\ge 1$ with $a_0$ and $a_{-1}$.
As Somos suggested in comment, we use the exponential generating function.
Let
\begin{align}
f(x) &= a_{-1} + a_0 x + \frac{1}{2}a_1x^2 + \cdots = \sum_{k=-1}^\infty \frac{a_k}{(k+1)!}x^{k+1},\\
g(x) &= a_{-1}x + \frac{1}{2}a_0x^2 + \frac{1}{6}a_1x^3 + \cdots = \sum_{k=-1}^\infty \frac{a_k}{(k+2)!}x^{k+2}, \\
h(x) &= \frac{1}{2}a_{-1}x^2 + \frac{1}{6}a_0x^3 + \frac{1}{24}a_1x^4 + \cdots = \sum_{k=-1}^\infty \frac{a_k}{(k+3)!}x^{k+3}.
\end{align}
We have
\begin{align}
f(x) &= a_{-1} + a_0 x + \sum_{k=1}^\infty \frac{a_k}{(k+1)!}x^{k+1}\\
&= a_{-1} + a_0 x + \sum_{k=1}^\infty \frac{(n+k)a_{k-1} - ka_{k-2}}{(k+1)!}x^{k+1}\\
&= a_{-1} + a_0 x + \sum_{k=1}^\infty \frac{na_{k-1}}{(k+1)!}x^{k+1}
+ \sum_{k=1}^\infty \frac{ka_{k-1}}{(k+1)!}x^{k+1} - \sum_{k=1}^\infty \frac{ka_{k-2}}{(k+1)!}x^{k+1}\\
&= a_{-1} + a_0 x + \sum_{k=1}^\infty \frac{na_{k-1}}{(k+1)!}x^{k+1}
+ \sum_{k=1}^\infty \frac{(k+1)a_{k-1}}{(k+1)!}x^{k+1} - \sum_{k=1}^\infty \frac{a_{k-1}}{(k+1)!}x^{k+1}\\
&\qquad - \sum_{k=1}^\infty \frac{(k+1)a_{k-2}}{(k+1)!}x^{k+1} + \sum_{k=1}^\infty \frac{a_{k-2}}{(k+1)!}x^{k+1}\\
&= a_{-1} + a_0 x + n\sum_{k=1}^\infty \frac{a_{k-1}}{(k+1)!}x^{k+1}
+ x\sum_{k=1}^\infty \frac{a_{k-1}}{k!}x^{k} - \sum_{k=1}^\infty \frac{a_{k-1}}{(k+1)!}x^{k+1}\\
&\qquad - x\sum_{k=1}^\infty \frac{a_{k-2}}{k!}x^{k} + \sum_{k=1}^\infty \frac{a_{k-2}}{(k+1)!}x^{k+1}\\
&= a_{-1} + a_0 x + n(g(x) - a_{-1}x) + x(f(x) - a_{-1}) - (g(x) - a_{-1}x) - xg(x) + h(x).
\end{align}
By taking derivative on both sides, noting that $h'(x) = g(x), g'(x) = f(x)$, we get the ODE
$$f'(x) = \frac{n-x}{1-x}f(x) + \frac{a_0-na_{-1}}{1-x},\ f(0)=a_{-1}.\tag{1}$$
The general solution of $f'(x) = \frac{n-x}{1-x}f(x)$ is $f_c(x) = C_0\frac{\mathrm{e}^x}{(1-x)^{n-1}}$.
Using the method variation of constants, we get the particular solution
$$f_p(x) = \frac{(a_0-na_{-1})\mathrm{e}^x}{(1-x)^{n-1}}\int_0^x (1-t)^{n-2}\mathrm{e}^{-t} \mathrm{d} t.$$
Thus, the solution of (1) is
$$f(x) = f_c(x) + f_p(x) = C_0\frac{\mathrm{e}^x}{(1-x)^{n-1}}
+ \frac{(a_0-na_{-1})\mathrm{e}^x}{(1-x)^{n-1}}\int_0^x (1-t)^{n-2}\mathrm{e}^{-t} \mathrm{d} t.$$
Since $f(0) = a_{-1}$, we have $C_0 = a_{-1}$. Thus, the final solution is
$$f(x) = \frac{a_{-1}\mathrm{e}^x}{(1-x)^{n-1}}
+ \frac{(a_0-na_{-1})\mathrm{e}^x}{(1-x)^{n-1}}\int_0^x (1-t)^{n-2}\mathrm{e}^{-t} \mathrm{d} t.$$
Here are some examples.
Example 1: When $n=1$, $a_{-1} = 1, a_0 = 1$, we have $a_k = 1, \forall k$.
The ODE is: $f'(x)=f(x), f(0) = 1$. The solution is $f(x) = \mathrm{e}^x$.
Example 2: When $n=1$, $a_{-1} = 1, a_0 = 2$, we have $a_1 = 3, a_2 = 5, a_3 = 11, a_4 = 35, a_5 = 155, \cdots$.
The ODE is: $f'(x) = f(x) + \frac{1}{1-x}, f(0) = 1$.
The solution is
$f(x) = \mathrm{e}^{x-1}\mathrm{Ei}(1, x-1) - \mathrm{e}^{x-1}\mathrm{Ei}(1,-1) + \mathrm{e}^x$.
Example 3: When $n=2$, $a_{-1} = 1$, $a_0 = 3$, we have $a_1 = 8, a_2 = 26, a_3 = 106, a_4 = 532, a_5 = 3194, \cdots$.
The ODE is: $f'(x) = \frac{2-x}{1-x} f(x) + \frac{1}{1-x}$, $f(0) = 1$.
The solution is $f(x) = \frac{2\mathrm{e}^x - 1}{1-x}$.