7.3. Steps in the Method of Undetermined Coefficients

Solve the homogeneous equation: ay′′ + by′ + cy = 0. Find the complementary function y_c by solving the auxiliary equation: ar² + br + c = 0. The nature of roots (real and distinct, real and equal, or complex) determines the form of y_c.

In this first step, we begin by solving the associated homogeneous equation, which is an equation without the non-homogeneous term. The goal is to find the complementary function, y_c. To do this, we write down the auxiliary equation (which is derived from the homogeneous equation) and solve for the roots. These roots can either be real and different, real and the same, or complex. Each type of root gives a different form for y_c. If the roots are real and distinct, for example, the complementary function takes the form of a linear combination of exponentials. If the roots are complex, it will include sine and cosine functions.

You can think of finding the complementary function as discovering the basic modes of vibration of a musical instrument. Different instruments produce different tones based on their construct (analogous to the nature of the roots) — for instance, a guitar and a violin can both vibrate but produce distinct sounds based on their design.

Assume a trial solution for y_p based on the form of f(x), using undetermined coefficients. The forms of trial solutions are:

• P(x) (polynomial of degree n) => y_p = Ax^n + Bx^{n-1} + ... + C
• e^(kx) => y_p = Ae^(kx)
• sin(ax), cos(ax) => y_p = Acos(ax) + Bsin(ax)
• x^m e^(kx) => y_p = (Ax^m + ... + C)e^(kx)
• e^(kx)sin(ax), e^(kx)cos(ax) => y_p = e^(kx)(Acos(ax) + Bsin(ax)).

In this step, we propose a form for the particular integral solution, y_p, based on the type of the forcing function f(x). We create a trial function that mimics the same structure as f(x). For example, if f(x) is a polynomial of degree n, our trial solution will be a polynomial of the same degree with undetermined coefficients (A, B, C, etc.). This trial function will allow us to approximate the response of the system to the external forces represented by f(x). The trial form changes based on whether f(x) is polynomial, exponential, or sinusoidal.

Imagine you're trying to predict the future sales of a product based on past performance. If sales data shows a quadratic trend, you'd build a model (a polynomial form) that fits that trend to forecast future sales. In this method, we craft our 'guess' in a similar way, based on the visible trend indicated by f(x).

If any term in the trial solution already appears in y_c, multiply the entire trial solution by x (or x², if needed) to eliminate duplication. This adjustment is known as the 'annihilator approach' or 'repetition rule.'

In this step, we need to ensure that our trial solution for y_p does not overlap with the complementary function y_c. If it does, we modify the trial solution by multiplying it by the smallest power of x necessary to eliminate that overlap. This technique is referred to as the 'annihilator approach,' and it helps ensure that the different parts of our solution are linearly independent. This adjustment is critical because if we did not make this change, we would not properly account for the effects of the non-homogeneous part on the system.

Think of stacking books on a shelf. If one book is too similar in size to another already on the shelf, it wouldn't fit well without creating clutter. So you might stack it on top of the other, using a box underneath to create space. In the context of the trial solution, we create space by multiplying our guessed function to avoid 'clutter' or duplicates with y_c.

Substitute the guessed y_p into the original non-homogeneous differential equation. Compare both sides and equate the coefficients of like terms to determine the unknown constants.

During this crucial step, we substitute our modified trial solution y_p back into the original non-homogeneous differential equation. Once substituted, we compare the left-hand side and the right-hand side of the equation. By aligning like terms, we can solve for the unknown coefficients (such as A, B, C) in our trial solution. This process of equating coefficients helps establish the exact influence of the non-homogeneous term on the solution.

Consider a recipe where you need to mix ingredients to achieve the right taste. After guessing the ratio of sugar to flour, you taste the mixture (substitution). If it's too sweet, you adjust the sugar content until it matches your ideal flavor (comparing and adjusting coefficients). In mathematics, this adjustment helps us refine our guess to fit the equation perfectly.

Finally, the general solution is obtained by combining the complementary function and the particular solution: y(x) = y_c(x) + y_p(x).

In the final step, we combine the solutions found in the previous steps. The general solution to the non-homogeneous differential equation is created by adding the complementary function y_c (found from the homogeneous part) to the particular solution y_p (derived through our guesses). This gives us a complete solution that incorporates both the natural behavior of the system and the effect of the external forces.

This step is like putting together a complete puzzle after finding individual pieces. Each piece represents a different aspect: the complementary function captures the inherent characteristics of the system, while the particular solution shows how it reacts to external influences. Together, they create the full picture of the situation you are analyzing.

--