Particular Solution Of Nonhomogeneous Differential Equation

Finding a particular solutionof a nonhomogeneous differential equation is a crucial step in solving many physical and engineering problems where external forces, inputs, or source terms are present. Unlike the homogeneous case, where the solution depends only on the system’s inherent properties, the particular solution captures the response driven by the non‑homogeneous term. Mastering this concept allows students and professionals to construct the full solution by combining the complementary (homogeneous) part with a specific answer that satisfies the forcing function.

Introduction

A linear differential equation can be written in the form

[ L[y] = y^{(n)} + a_{n-1}(x)y^{(n-1)} + \dots + a_0(x)y = g(x), ]

where (L) is a linear differential operator and (g(x)\neq 0) is the nonhomogeneous term. The general solution is

[ y(x)=y_h(x)+y_p(x), ]

with (y_h) solving the associated homogeneous equation (L[y]=0) and (y_p) a particular solution of the full equation. The particular solution is not unique; any two particular solutions differ by a solution of the homogeneous equation. So naturally, the goal is to find one convenient (y_p) that satisfies (L[y_p]=g(x)).

Methods for Finding a Particular Solution

Two widely used techniques are the Method of Undetermined Coefficients and Variation of Parameters. Each has its own domain of applicability and procedural steps.

Method of Undetermined Coefficients

This method works when (g(x)) is a linear combination of functions whose derivatives are of the same type—namely, exponentials, polynomials, sines, and cosines. The procedure is:

1. Solve the homogeneous equation to obtain (y_h).
2. Guess the form of (y_p) based on (g(x)). As an example, if (g(x)=e^{2x}), try (y_p=Ae^{2x}); if (g(x)=x^2), try (y_p=Ax^2+Bx+C).
3. Adjust for duplication: if any term of the guess appears in (y_h), multiply the guess by (x) (or a higher power of (x)) until no overlap remains.
4. Substitute the guess into the original differential equation and equate coefficients of like terms to solve for the unknown constants.
5. Write the particular solution with the determined constants.

Because the guess is anchored in the same function family as (g(x)), the algebra remains straightforward, making this method ideal for constant‑coefficient equations with simple forcing functions.

Variation of Parameters

When (g(x)) contains functions outside the undetermined‑coefficients repertoire (e.g., (\ln x), (\tan x), or variable coefficients), variation of parameters provides a general formula.

1. Find a fundamental set ({y_1(x),y_2(x),\dots,y_n(x)}) of solutions to the homogeneous equation.
2. Construct the Wronskian (W(x)=\det\big[y_i^{(j-1)}\big]).
3. Set up the particular solution as

[ y_p(x)=\sum_{i=1}^{n} y_i(x)\int \frac{g(x),W_i(x)}{W(x)},dx, ]

where (W_i) is the Wronskian obtained by replacing the (i)-th column of (W) with the vector ([0,\dots,0,g(x)]^T).
Also, 4. Evaluate the integrals (often requiring integration techniques such as substitution or integration by parts).
5. Add (y_p) to the homogeneous solution to obtain the general answer.

Although more computationally intensive, variation of parameters is universally applicable to linear ODEs with continuous coefficients.

For more on this topic, read our article on year 11 biology syllabus 2024 or check out who was the first to propose the existence of atoms.

Scientific Explanation

Why the Particular Solution Matters

The homogeneous solution (y_h) describes the system’s natural behavior—its response to initial conditions in the absence of external influence. The particular solution (y_p), however, encodes the forced response directly attributable to the source term (g(x)). In physical terms, if the differential equation models a mass‑spring‑damper system subjected to an external force (F(t)), then (y_h) captures free vibrations, while (y_p) represents the steady‑state motion caused by (F(t)). Ignoring (y_p) would omit the effect of the forcing altogether, leading to incomplete or incorrect predictions.

Superposition Principle

Because the differential operator (L) is linear, the principle of superposition holds: if (y_{p1}) solves (L[y]=g_1(x)) and (y_{p2}) solves (L[y]=g_2(x)), then (c_1y_{p1}+c_2y_{p2}) solves (L[y]=c_1g_1(x)+c_2g_2(x)). This property justifies breaking down a complicated (g(x)) into simpler components, finding a particular solution for each, and summing the results. It also explains why the method of undetermined coefficients can treat sums of exponentials, polynomials, and trigonometric functions term by term.

Step‑by‑Step Example

Consider the second‑order constant‑coefficient equation

[ y'' - 3y' + 2y = 5e^{2x}. ]

1. Homogeneous solution: Solve (r^2-3r+2=0) → (r=1,2). Hence (y_h=C_1e^{x}+C_2e^{2x}).
2. Guess for (y_p): Since (g(x)=5e^{2x}) is an exponential, we would normally try (y_p=Ae^{2x}). On the flip side, (e^{2x}) already appears in (y_h), so we multiply by (x): (y_p=Axe^{2x}).
3. Compute derivatives:

[ y_p' = A e^{2x} + 2Ax e^{2x},\qquad y_p'' = 4A e^{2x} + 4Ax e^{2x}.

4. Substitute into the equation:

[ (4A e^{2x} + 4Ax e^{2x}) - 3(A e^{2x} + 2Ax e^{2x}) + 2(Ax e^{2x}) = 5e^{2x}. ]

Simplify:

[ 4A e^{2x} + 4Ax e^{2x} - 3A e^{2x} - 6Ax e^{2x} + 2Ax e^{2x} = 5e^{2x}. ]

Group terms:

[ (4A - 3A) e^{2x} + (4A - 6A + 2A)x e^{2x} = 5e^{2x}. ]

The (x e^{2x}) terms cancel, leaving:

[ A e^{2x} = 5e^{2x} \quad\Rightarrow\quad A = 5. ]

Thus (y_p = 5x e^{2x}).

5. General solution:

[ y(x) = C_1 e^{x} + C_2 e^{2x} + 5x e^{2x}. ]

This example illustrates how the particular solution captures the system's forced response, while the homogeneous part reflects natural dynamics.

Conclusion

Finding a particular solution to a linear differential equation is a fundamental step in solving initial value and boundary value problems. In practice, whether through the straightforward method of undetermined coefficients or the universally applicable variation of parameters, the goal is to identify the forced response (y_p) that complements the homogeneous solution (y_h). Understanding the physical meaning of (y_p) as the system's reaction to external inputs, along with the mathematical structure provided by linearity and superposition, equips one to tackle a wide array of problems in physics, engineering, and applied mathematics. Mastery of these techniques ensures accurate modeling of dynamic systems, from mechanical vibrations to electrical circuits, and lays the groundwork for more advanced analytical and numerical methods.