How To Find Particular Solution Of A Differential Equation

How to Find the Particular Solution of a Differential Equation: A Complete Guide

Solving a differential equation is often described as finding its general solution, a family of curves containing all possible solutions. However, for most real-world applications—from predicting the motion of a pendulum to modeling population growth—we need one specific solution that fits our initial conditions or physical constraints. This specific solution is the particular solution. Understanding how to isolate it from the general solution is a fundamental skill in applied mathematics, physics, and engineering. This guide will demystify the process, breaking down the theory, methods, and common pitfalls into a clear, actionable framework. By the end, you will not only know how to find a particular solution but also why each step matters and how to apply it to tangible problems.

Detailed Explanation: What is a Particular Solution?

To grasp the particular solution, we must first distinguish between two core components of solving non-homogeneous linear differential equations. A standard second-order linear equation has the form: a(x)y'' + b(x)y' + c(x)y = g(x) Here, g(x) is the non-homogeneous term or forcing function, representing an external influence. The equation is split into two parts:

1. The homogeneous equation: a(x)y'' + b(x)y' + c(x)y = 0. Its solution, y_h(x) or y_c(x) (the complementary solution), contains arbitrary constants (e.g., C1, C2) and describes the system's natural, unforced behavior.
2. The particular solution, y_p(x). This is any single function that satisfies the entire non-homogeneous equation. Crucially, it contains no arbitrary constants. It models the specific, steady-state response of the system directly caused by the forcing function g(x).

The general solution to the non-homogeneous equation is the sum of these two parts: y(x) = y_h(x) + y_p(x) The complementary solution y_h accounts for the system's intrinsic dynamics (like a pendulum's natural frequency), while y_p accounts for the driven response (like a periodic push). Finding y_p is the critical step that connects the abstract equation to the physical forcing. Without it, we have a description of all possible motions but not the one that actually occurs under a given push.

Step-by-Step Breakdown: Primary Methods for Finding y_p

The method chosen depends heavily on the form of the forcing function g(x). For linear differential equations with constant coefficients, two primary techniques dominate.

Method 1: Method of Undetermined Coefficients

This elegant guessing method works when g(x) is a simple, differentiable function: a polynomial, an exponential (e^(αx)), a sine/cosine (sin(βx), cos(βx)), or a product of these. The core idea is to assume a trial form for y_p that mirrors g(x), but with undetermined coefficients.

Step-by-Step Process:

1. Examine g(x). Identify its "type" (polynomial, exponential, sinusoidal, or combination).
2. Write a trial y_p. For g(x) = 5x^3, guess y_p = Ax^3 + Bx^2 + Cx + D. For g(x) = 7e^(2x), guess y_p = Ae^(2x). For g(x) = 3sin(4x), guess y_p = Asin(4x) + Bcos(4x) (you often need both sin and cos even if only sin appears in g(x)).
3. Check for Resonance (The Crucial Step). If any term in your trial y_p already appears in the complementary solution y_h, you must multiply your entire trial y_p by x^s, where s is the smallest positive integer that makes the new trial form linearly independent from y_h. This is the most common source of error.
4. Substitute and Solve. Plug the trial y_p (with its derivatives) into the original non-homogeneous equation. All terms will collapse to an expression involving only your undetermined coefficients (A, B, C, etc.) and x. Collect like terms.
5. Equate Coefficients. Set the coefficients of each power of x (or e^(αx), sin(βx), etc.) on the left equal to those on the right (which come from g(x)).