Generally speaking, you can solve any non-homogeneous linear recurrence relations with constant coefficients using several methods depending on the recurrence formula.
First-order non-homogeneous linear recurrences: guessing, iteration or telescoping
A first-order non-homogeneous linear recurrence has the form $u_n=au_{n-1}+p\left(n\right)\ \text{with}\ a\in\mathbb{R}$ and you can often guess the closed formula or use iteration to find a recurring pattern. Also, first-order recurrence relations always telescope to a simple sum, starting with $u_n-u_{n-1}$ you can build the corresponding telescoping sum and determine the closed formula. Then, induction is naturally used to prove the found formula.
There is also a general method to determine the closed formula of non-homogeneous linear recurrences of the form $u_n=au_{n-1}+b\ \text{with}\ a,b\in\mathbb{R}$ using the fixed point theorem, let the fixed point $\alpha$ satisfies $\alpha=a\alpha+b$ and you have
$$u_n-\alpha=a\left(u_{n-1}-\alpha\right)$$
This theorem can then be used to construct a translated sequence $v_n=u_n-\alpha$ to determine the closed formula for $u_n$ depending on the initial value. However, this elegant method is not mentioned in most calculus books and there is not many online resources on the topic (except this blog article). I was introduced to this method during my high school in France.
Note: this approach is very similar to the one used for solving second-order and higher non-homogeneous linear recurrences. Namely solving the characteristic polynomial equation for the homogeneous recurrence relation and then finding the particular solution, which turns out to be $\alpha$ in this situation.
Second-order and higher non-homogeneous linear recurrences: the characteristic polynomial equation
For second-order and higher order recurrence relations, trying to guess the formula or use iteration will usually result in a lot of frustration. Also, these recurrence relations will usually not telescope to a simple sum. The idea here is to solve the characteristic polynomial equation associated with the homogeneous recurrence relation.
Note: the characteristic polynomial equation is directly derived from generating series, the proof is admitted here.
A second-order non-homogeneous linear recurrence has the form $u_n=au_{n-2}+bu_{n-2}+p\left(n\right)\ \text{with}\ a,b\in\mathbb{R}\ \text{and}\ p\left(n\right)\neq0$ and from there you can distinguish its associated homogeneous relation $u_n=au_{n-2}+bu_{n-2}$ and the non-homogeneous part $p\left(n\right)$. The solution of a non-homogeneous linear recurrence relation has thus two parts.
- The solution $\left\{u_n^H\right\}$ of the associated homogeneous recurrence relation $u_n=au_{n-2}+bu_{n-2}$
- The solution $\left\{u_n^P\right\}$ of the non-homogeneous part $p\left(n\right)$ called the particular solution
We eventually have the final solution $\left\{u_n^H+u_n^P\right\}$ as a combination of the two previous solutions.
The homogeneous solution can be found using the roots of the associated characteristic polynomial equation. The particular solution is then found based on the homogeneous solution and chosen to satisfy the initial values $u_0$ and $u_1$.
Another method to solve a non-homogeneous recurrence is the method of symbolic differentiation. If $p\left(n\right)$ is a polynomial with degree $r$, then this non-homogeneous recurrence can be reduced to a homogeneous recurrence by applying the method of symbolic differentiation $r$ times. You can then find the roots of the characteristic polynomial in order to solve the relation.
I will edit this post and add more content when I have more free time.