Find A Particular Solution Of The Differential Equation

Finding a Particular Solution: Your Key to Solving Real-World Differential Equations

Imagine you are an engineer designing a suspension system for a car. You know the general behavior of springs and dampers—they can oscillate in many ways depending on initial conditions. But your specific, pressing question is: How will the car’s body respond to a specific, ongoing bumpy road? You need the answer for that particular input, not all possible motions. This is the essence of finding a particular solution to a differential equation. While the general solution describes the complete family of all possible motions (the "homogeneous" or natural response), a particular solution is a single, specific function that satisfies the entire non-homogeneous equation, accounting for an external driving force or input. It represents the steady-state response or forced response of the system to that specific stimulus. Mastering how to find this particular solution is what transforms abstract mathematics into a practical tool for modeling physics, engineering, economics, and biology.

Detailed Explanation: The Two-Part Symphony of Solutions

To understand a particular solution, we must first frame the problem it solves. A linear non-homogeneous differential equation has the standard form: L[y] = g(x) Here, L is a linear differential operator (like d²/dx² + 3d/dx + 2), y is the unknown function we seek, and g(x) is a non-zero function representing an external input, force, or source term. The equation is "non-homogeneous" precisely because of this g(x).

The Superposition Principle for linear operators is our guiding light. It tells us that the general solution y_general to the non-homogeneous equation is the sum of two distinct parts: y_general = y_c + y_p

1. y_c (Complementary Solution): This is the solution to the associated homogeneous equation L[y] = 0. It represents the system's natural, unforced behavior—the "free" oscillations or decay that depend entirely on initial conditions (like the initial displacement and velocity of our car's suspension). It contains arbitrary constants (C₁, C₂, etc.).
2. y_p (Particular Solution): This is any single function that satisfies the full, original non-homogeneous equation L[y] = g(x). It is not arbitrary; it is specifically shaped by the form of the input g(x). It represents the system's long-term, sustained response to the persistent external force g(x). It contains no arbitrary constants.

Therefore, finding y_p is the critical step in solving the non-homogeneous equation. Once y_p is found, the full general solution is completed by adding the complementary solution y_c. The arbitrary constants in y_c are then determined by applying the initial conditions or boundary conditions of the specific problem, yielding the unique solution for that scenario.

Step-by-Step Breakdown: Two Primary Methods

There are two primary, systematic methods for finding a particular solution for linear equations with constant coefficients, which are the most common in applications.

Method 1: Method of Undetermined Coefficients

This method is elegant and efficient but works only when g(x) is a combination of polynomials, exponentials (e^(αx)), sines (sin(βx)), and cosines (cos(βx)). The strategy is to guess the form of y_p based on the form of g(x), using "undetermined coefficients" as placeholders, then plug the guess into the equation to solve for those coefficients.

Step-by-Step Process:

1. Examine g(x). Identify its components (e.g., a polynomial, e^(2x), sin(3x)).
2. Make an initial guess for y_p. For each component in g(x), use a corresponding template:
   • If g(x) is a polynomial of degree n, guess a polynomial of degree n: A_n x^n + ... + A_1 x + A_0.
   • If g(x) is P(x)e^(αx) (where P is a polynomial), guess Q(x)e^(αx), where Q is a polynomial of the same degree as P.
   • If g(x) is P(x)e^(αx)cos(βx) or P(x)e^(αx)sin(βx), guess e^(αx)[R(x)cos(βx) + S(x)sin(βx)], where R and S are polynomials of the same degree as P.
3. Check for Resonance (The Crucial Adjustment). If any term in your initial guess for y_p is already a solution to the homogeneous equation (i.e., appears in y_c), you must multiply your entire guess by x^k, where k is the smallest positive integer that eliminates this duplication. This accounts for the "resonance" or repeated root problem.
4. Substitute and Solve. Take the derivatives of your adjusted guess, substitute y_p and its derivatives into the original differential equation L[y] = g(x).
5. Equate Coefficients. Collect like terms on the left side and set the coefficients equal to those on the right side (g(x)). This yields a system of linear algebraic equations for your undetermined coefficients (A, B, etc.).
6. Solve for Coefficients. Solve the system to find the numerical values of the coefficients. Substitute them back into your guess to obtain the explicit y_p.

Example: Find a particular solution for y'' + y = 3 sin(2x).

• g(x) = 3 sin(2x). Initial guess: y_p = A cos(2x) + B sin(2x).
• Check y_c: The homogeneous equation y'' + y = 0 has solutions cos(x) and sin(x). Our guess cos(2x), sin(2x) is

Our guess (y_p = A\cos(2x)+B\sin(2x)) does not contain any term that solves the homogeneous equation—(\cos(x)) and (\sin(x)) are distinct from (\cos(2x)) and (\sin(2x)). Therefore no multiplication by (x) is required.

Compute the derivatives:

[ \begin{aligned} y_p' &= -2A\sin(2x)+2B\cos(2x),\[2pt] y_p'' &= -4A\cos(2x)-4B\sin(2x). \end{aligned} ]

Substitute into the left‑hand side of the differential equation:

[ y_p''+y_p = \bigl[-4A\cos(2x)-4B\sin(2x)\bigr] + \bigl[A\cos(2x)+B\sin(2x)\bigr] = (-3A)\cos(2x)+(-3B)\sin(2x). ]

Set this equal to the right‑hand side (3\sin(2x)):

[ (-3A)\cos(2x)+(-3B)\sin(2x)=0\cdot\cos(2x)+3\sin(2x). ]

Equating coefficients gives the linear system

[ \begin{cases} -3A = 0,\ -3B = 3, \end{cases} \qquad\Longrightarrow\qquad A = 0,; B = -1. ]

Thus a particular solution is

[ y_p = -\sin(2x). ]

Method 2: Variation of Parameters

When the forcing term (g(x)) does not belong to the elementary families listed above—or when the coefficients of the differential equation are not constant—undetermined coefficients either fails or becomes cumbersome. In such cases variation of parameters provides a systematic alternative that works for any continuous (g(x)).

Core Idea

Assume a particular solution of the form

[ y_p(x)=u_1(x)y_1(x)+u_2(x)y_2(x)+\dots+u_n(x)y_n(x), ]

where (y_1,\dots,y_n) are a fundamental set of solutions to the homogeneous equation and (u_i(x)) are unknown functions to be determined.

Derivation of the System for (u_i)

To avoid an overdetermined system, impose the constraints

[ \begin{aligned} u_1' y_1 + u_2' y_2 + \dots + u_n' y_n &= 0,\ u_1' y_1' + u_2' y_2' + \dots + u_n' y_n' &= 0,\ &\ \vdots\ u_1' y_1^{(n-1)} + u_2' y_2^{(n-1)} + \dots + u_n' y_n^{(n-1)} &= 0. \end{aligned} ]

These (n-1) equations, together with the original differential equation, yield a linear system for the derivatives (u_i'). Solving this system (often via Cramer's rule) gives each (u_i'), which can then be integrated to obtain (u_i(x)). Finally, substitute the (u_i) back into the expression for (y_p).

Example (Second‑order case)

Consider the same equation (y''+y=3\sin(2x)) but now treat it with variation of parameters.

1. Homogeneous solutions: (y_1=\cos x,; y_2=\sin x).
2. Form the Wronskian: [ W(y_1,y_2)=\begin{vmatrix} \cos x & \sin x\ -\sin x & \cos x \end{vmatrix}= \cos^2 x+\sin^2 x = 1. ]
3. Compute (u_1') and (u_2') using the standard formulas for a second‑order equation: [ u_1' = -,\frac{y_2,g(x)}{W}= -\sin x \cdot 3\sin(2x)= -3\sin x\sin(2x), ] [ u_2' = ;;\frac{y_1,g(x)}{W}= \cos x \cdot 3\sin(2x)= 3\cos x\sin(2x). ]
4. Integrate: [ u_1(x)= -\int 3\sin x\sin(2x),dx,\qquad u_2(x)= \int 3\cos x\sin(2x),dx. ] Using product‑to‑sum identities, [ \sin x\sin(2x)=\tfrac12[\cos(x)-\cos(3x)],\qquad \cos x\sin(2x)=\tfrac12[\sin(3x)+\sin(x)]. ] Hence [ u_1(x)= -\tfrac32!\int![\cos x-\cos 3x]dx =