 Home
 >
 Non Homogeneous Differential Equation – Solutions and Examples
Non Homogeneous Differential Equation – Solutions and Examples
Learning about nonhomogeneous differential equations is fundamental since there are instances when we’re given complex equations with functions on both sides of the equation. Laws of motion, for example, rely on nonhomogeneous differential equations, so it is important that we learn how to solve these types of equations.
Nonhomogenous differential equations contain functions on the righthand side of the equations. We can find their solutions by writing down the general solution of the associated homogeneous differential equation and the particular solution of the nonhomogeneous term.
Solving nonhomogeneous differential equations will still require our knowledge on solving second order homogeneous differential equations, so keep your notes handy on characteristic and second order homogeneous equations. This article covers the fundamentals needed to identify nonhomogeneous differential equations and two important methods that will help you find their solutions.
What Is a Non Homogeneous Differential Equation?
Nonhomogeneous differential equations are simply differential equations that do not satisfy the conditions for homogeneous equations. In the past, we’ve learned that homogeneous equations are equations that have zero on the righthand side of the equation.
This means that nonhomogenous differential equations are differential equations that have a function on the righthand side of their equation.
Here are some examples of homogeneous and nonhomogenous differential equations. Through these examples, we’ll learn how to identify differential equations based on their form.
Homogeneous Differential Equations  NonHomogenous Differential Equations 
\begin{aligned}y^{\prime \prime} + y^{\prime} – 2y = 0\end{aligned}  \begin{aligned}y^{\prime \prime} + y^{\prime} – 2y = 4x\end{aligned} 
\begin{aligned}y^{\prime \prime} 3y^{\prime} + y = 0\end{aligned}  \begin{aligned}y^{\prime \prime} 3y^{\prime} + y = 6 – 2x + 3x^2\end{aligned} 
\begin{aligned}y^{\prime \prime \prime} +4y^{\prime\prime} + 4y^{\prime} + y = 0\end{aligned}  \begin{aligned}y^{\prime \prime \prime} +4y^{\prime\prime} + 4y^{\prime} + y = 2e^x\end{aligned} 
From these three equations alone, we can clearly see that the righthand side of the equations will determine whether the differential equations are homogeneous or nonhomogeneous. We call the equations on the left column the associated homogeneous equations, since they share identical expressions with their nonhomogeneous counterparts.
Since we’ve been working on first order and second order nonhomogeneous linear differential equations, let us show you their general forms:
\begin{aligned}\textbf{First Order} &: y^{\prime} + P(x)y = f(x)\\\textbf{Second Order} &: y^{\prime} + P(x)y^{\prime} + Q(x)y = g(x)\\ \end{aligned}
In the next section, we’ll show you how to solve these types of equations by applying old techniques and even new methods!
How To Solve Non Homogeneous Differential Equations?
We can solve nonhomogeneous linear differential equations by finding the general solution of the associated homogeneous differential equation, $y_h$, and the particular solution of the nonhomogeneous equation, $y_p$.
Let’s work on a first order nonhomogeneous differential equation to start this section lightly.
\begin{aligned}y^{\prime} + P(x)y = f(x)\end{aligned}
In the past, we’ve learned how to deal with these types of equations. We’ve shown you how to use integrating factors to write the general equation for a first order nonhomogeneous differential equation. This process will eventually lead to a general solution shown below.
\begin{aligned}y &= h(x) + p(x) \end{aligned}
For this solution, $g(x)$ represents the general solution, while $h(x)$ represents the particular solution of the first order nonhomogeneous differential equation. Since we’ve set a special article for that, our discussion will focus on second order nonhomogeneous differential equations and those with a higher order as well.
SOLUTION FOR NON HOMOGENEOUS DIFFERENTIAL EQUATIONS Suppose that we have a second order nonhomogeneous linear differential equation shown below. \begin{aligned}y^{\prime \prime} +ay^{\prime} + by &= g(x)\end{aligned} For this equation, $a$ and $b$ are constants. When this happens, the nonhomogeneous differential equation with a general solution as shown below. \begin{aligned}\boldsymbol{y} &= \boldsymbol{y_h + y_p}\\\\y_h&: \text{General Solution}\\y_p&: \text{Particular Solution}\end{aligned} This can be extended with $n$thorder nonhomogeneous linear differential equation. The general solution of this linear differential equation is as shown below. \begin{aligned}c_1y_1 + c_2\end{aligned}

In previous articles, we have also learned how to find the general solutions for homogeneous differential equations. We’ll apply the same methods and techniques when finding the general solution for our associated homogeneous equation. This means that we’ll be focusing on techniques to find the particular solution for these nonhomogeneous equations.
How To Find the Particular Solution of a Non Homogeneous Differential Equation
The two most common methods when finding the particular solution of a nonhomogeneous differential equation are: 1) the method of undetermined coefficients and 2) the method of variation of parameters.
Method of Undetermined Coefficients
Let’s first break down the steps for the method of undetermined coefficients and know when it’s best to use this technique. The undetermined coefficient method works best when the righthand side of our nonhomogeneous differential equation is a function that can be written as sum or product of the following functions: $x^n$, $e^{ax}$,$\sin \beta x$, or $\cos \beta x$.
\begin{aligned}y^{\prime \prime} +ay^{\prime} + by &= g(x)\end{aligned}
Once we identify the form of the $g(x)$, use a strategic guess for the particular solution, $y_p$. Let us show you some common examples for a smart guess given an expression for $g(x)$:
Example of $\boldsymbol{g(x)}$  Particular Solution ($\boldsymbol{y_p}$) 
\begin{aligned}g(x) = 2x^2\end{aligned}  \begin{aligned}y_p = ax^2 + bx + c\end{aligned} 
\begin{aligned}g(x) = 3xe^x\end{aligned}  \begin{aligned}y_p = axe^x + bex^x c\end{aligned} 
\begin{aligned}g(x) = x + 4\cos 3x\end{aligned}  \begin{aligned}y_p = (ax + b) + c \sin 3x + d + \cos 3x\end{aligned} 
These three equations should give you an idea of how we come up with particular solutions we can use to complete the general solution for our second order homogeneous differential equation. Find the first and second derivatives of $y_p$ and use them inside the equation then solve for the values of the constants.
To understand how this method works, let’s try to solve the nonhomogeneous differential equations, $y^{\prime \prime} + 6y^{\prime} + 5y = 4x$. First, let’s find $y_h$, the general solution of the equation’s associated homogeneous equation.
\begin{aligned}y^{\prime \prime} + 6y^{\prime} + 5y = 0\end{aligned}
This means that the characteristic equation (or auxiliary equation) is a quadratic equation, $r^2 + 6r + 5= 0$. The roots of this equation are $r_1 = 1$ and $r_2 = 5$, so the general form of $y_h$ is equal to the equation shown below.
\begin{aligned}y_h = C_1 e^{x} + C_2 e^{5x} \end{aligned}
Since our equation’s righthand side is equal to $g(x) = 4x$, the general form of our particular solution would be equivalent to $y_p = ax + b$. This means that $y_p^{\prime } = a$ and $y_p^{\prime } = 0$. We’ve learned in the past that the particular solution should be a solution to the differential equation, so substitute these expressions and solve for $a$ and $b$.
\begin{aligned}y^{\prime \prime}_p + 6y^{\prime}_p + 5y_p &= 4x\\0 + 6a + 5(ax +b) &= 4x\\5ax + (6a +b) &= 4x\\\\5a= 4, 6a &+b =0\\a = \dfrac{4}{5} ,b =& \dfrac{24}{5}\end{aligned}
This means that $y_p = \dfrac{4}{5}x – \dfrac{24}{5}$. Now that we have both components, $y_h$ and $y_p$, let’s complete the general solution of our second order nonhomogeneous differential equation is:
\begin{aligned}y_h&= C_1 e^{x} + C_2 e^{5x} \\y_p&=\dfrac{4}{5}x – \dfrac{24}{5}\\\\ \boldsymbol{y(x)} &\boldsymbol{=} \boldsymbol{C_1 e^{x} + C_2 e^{5x}+\dfrac{4}{5}x – \dfrac{24}{5}} \end{aligned}
We’ve now shown you how the first method works, so let’s move on to our second technique: the method of variation of parameters.
Method of Variation Parameters
For more complex expressions for $g(x)$, it may be difficult to come up with a particular solution that we can easily use. When that happens, it is much better for us to apply the method of variation parameters, where we assume that the general and particular solution components have the same forms. Below are the guidelines to remember when using this method:
 Find the expression for $y_h = C_1y_1 + C_2y_2$, using previous techniques.
 Write down the particular solution, $y_p = ay_1 + by_2$, based on the form of $y_h$.
 Set up the system of equations to solve for $a^{\prime}$ and $b^{\prime}$.
\begin{aligned}a^{\prime}y_1 + b^{\prime}y_2 &= 0\\a^{\prime}y_1^{\prime} + b^{\prime}y_2^{\prime} &= g(x) \end{aligned}
 Integrate the resulting expressions for $a^{\prime}$ and $b^{\prime}$ to find $a$ and $b$.
 Write down the general solution for the nonhomogeneous differential equation: $y = y_h + y_p$.
These steps are straightforward but can be complex depending on the resulting expressions. Just apply the appropriate techniques learned in the past to find the solutions using variations of parameters. Our first example below shows how we can use this method to solve more interesting nonhomogeneous differential equations.
Example 1
Find the general solution of the nonhomogeneous differential equation, $y^{\prime\prime} +y = \tan x$.
Solution
Since we don’t have a guiding rule for $g(x) = \tan x$, we can’t use the method of undetermined coefficients. We can instead use the second method beginning with finding the general solution for the associated homogeneous equations. This means that the characteristic equation is equal to $r^2 + 1= 0 \rightarrow r = \pm i$, so the homogeneous solution is equal to
\begin{aligned}y_h = C_1\cos x + C_2 \sin x\end{aligned}
Using the form of $y_h$, let’s use the particular solution, $y_p = a\cos x +b\sin x$. Write down the system of linear equations:
\begin{aligned}a^{\prime}\cos x+ b^{\prime}\sin x &= 0\\a^{\prime}\sin x+ b^{\prime}\cos x &= \tan x \end{aligned}
Apply the elimination method by multiplying appropriate factors then find $a$.
\begin{aligned}{\color{blue}\sin x}(a^{\prime}\cos x+ b^{\prime}\sin x) &= {\color{blue}\sin x}(0)\\ {\color{blue}\cos x}(a^{\prime}\sin x+ b^{\prime}\cos x ) &= {\color{blue}\cos x}(\tan x)\\\\a^{\prime}\sin x\cos x + b^{\prime}\sin^2x&= 0\\a^{\prime}\sin x\cos x + b^{\prime}\cos^2x&= \sin x \end{aligned}
Simplify the equation further and you should come up with $b^{\prime} = \cos x – \sec x $. This eventually leads to $a^{\prime} = \dfrac{\sin^2 x}{\cos x} = \cos x – \sec x$. Integrate the expressions to find the constants for $y_p$.
\begin{aligned}\boldsymbol{a^{\prime} = \cos x – \sec x}\end{aligned}  \begin{aligned}\boldsymbol{\sin x}\end{aligned} 
\begin{aligned} a^{\prime} &= \int (\cos x – \sec x ) \phantom{x}dx\\&= \sin x – \ln \sec x +\tan x\end{aligned}  \begin{aligned} b^{\prime} &= \int \sin x \phantom{x}dx\\&= \cos x\end{aligned} 
\begin{aligned} y_p &= a\cos x +b\sin x\\&= \cos x(\sin x – \ln \sec x +\tan x) + \sin x(\cos x)\\&= \cos x \ln \sec x +\tan x\end{aligned} 
Combining the expressions for $y_h$ and $y_p$, we have $y = y_h + y_p = C_1\cos x + C_2 \sin x \cos x \ln \sec x +\tan x$.
Example 2
Find the general solution of the nonhomogeneous differential equation, $y^{\prime\prime\prime} + 6y^{\prime\prime} + 12y^{\prime} +8y = 4x$.
Solution
Our righthand side this time is $g(x) = 4x$, so we can use the first method: undetermined coefficients. We begin by finding the general solution for the associated homogeneous equations, $ y^{\prime\prime\prime} + 6y^{\prime\prime} + 12y^{\prime} +8y = 0$. This type of thirdorder differential equation will have a general solution of $y_h=C_1e^{r_1x}+ C_2xe^{r_2x}+ C_3x^2e^{r_3x}$, where $\{r_1, r_2, r_3\}$ are roots of the characteristic equation.
\begin{aligned} y_h= C_1e^{2x}+ C_2xe^{2x}+ C_3x^2e^{2x}\end{aligned}
Since the righthand side of the equation is equal to $4x$, the particular solution, $y_p$ will have a general form of $ax + b$. This means that $y_p^{\prime} =a$ and $y_p^{\prime \prime} =0$.
\begin{aligned} y^{\prime\prime\prime} + 6y^{\prime\prime} + 12y^{\prime} + 8y &= 4x\\0 + 6(0) + 12a + 8(ax +b) &= 4x\\8ax + (12a +8b) &= 4x\\\\8a= 4 &\rightarrow a = \dfrac{1}{2}\\6+8b =0 &\rightarrow b=\dfrac{3}{4}\end{aligned}
This means that we have the particular solution, $y_p = \dfrac{1}{2}x – \dfrac{3}{4}$. Combine the two components of our equation to find the general equation of the third order nonhomogeneous differential equation.
\begin{aligned}y&= y_h + y_p\\&= C_1e^{2x}+ C_2xe^{2x}+ C_3x^2e^{2x} +\dfrac{1}{2}x – \dfrac{3}{4}\end{aligned}
Practice Questions
1. Find the general solution of the nonhomogeneous differential equation, $y^{\prime\prime} +y^{\prime}6 = x^3$.
2. Find the general solution of the nonhomogeneous differential equation, $y^{\prime\prime} 2y^{\prime} = \dfrac{e^x}{4x}$.
3. Find the general solution of the nonhomogeneous differential equation, $ y^{\prime\prime\prime}+3y^{\prime\prime}+3y^{\prime} +y = 2x $.
Answer Key
1. $ y=C_1+C_2e^{x}+ \dfrac{x^4}{4} x^3+3x^2$
2.$ y= C_1e^x+C_2xe^x+ \dfrac{e^xx+e^xx\ln \left(x\right)}{4}$
3.$ y=C_1e^{x}+C_2xe^{x}+C_3x^2e^{x}+2x6$