Method of undetermined coefficients

Description of the method

Consider a linear non-homogeneous ordinary differential equation of the form

∑ i = 0 n c i y ( i ) + y ( n + 1 ) = g ( x ) ... where y ( i ) denotes the i-th derivate of y, and c i denotes a function of x

The method consists of finding the general homogeneous solution y c for the complementary linear homogeneous differential equation

∑ i = 0 n c i y ( i ) + y ( n + 1 ) = 0 ,

and a particular integral y p of the linear non-homogeneous ordinary differential equation based on g ( x ) . Then the general solution y to the linear non-homogeneous ordinary differential equation would be

y = y c + y p .

If g ( x ) consists of the sum of two functions h ( x ) + w ( x ) and we say that y p 1 is the solution based on h ( x ) and y p 2 the solution based on w ( x ) . Then, using a superposition principle, we can say that the particular integral y p is

y p = y p 1 + y p 2 .

Typical forms of the particular integral

In order to find the particular integral, we need to 'guess' its form, with some coefficients left as variables to be solved for. This takes the form of the first derivative of complementary function. Below is a table of some typical functions and the solution to guess for them.

If a term in the above particular integral for y appears in the homogeneous solution, it is necessary to multiply by a sufficiently large power of x in order to make the solution independent. If the function of x is a sum of terms in the above table, the particular integral can be guessed using a sum of the corresponding terms for y.

Examples

Example 1

Find a particular integral of the equation

y ″ + y = t cos ⁡ t .

The right side t cos t has the form

P n e α t cos ⁡ β t

with n = 1, α = 0, and β = 1.

Since α + iβ = i is a simple root of the characteristic equation

λ 2 + 1 = 0

we should try a particular integral of the form

y p = t [ F 1 ( t ) e α t cos ⁡ β t + G 1 ( t ) e α t sin ⁡ β t ] = t [ F 1 ( t ) cos ⁡ t + G 1 ( t ) sin ⁡ t ] = t [ ( A 0 t + A 1 ) cos ⁡ t + ( B 0 t + B 1 ) sin ⁡ t ] = ( A 0 t 2 + A 1 t ) cos ⁡ t + ( B 0 t 2 + B 1 t ) sin ⁡ t .

Substituting y_p into the differential equation, we have the identity

t cos ⁡ t = y p ″ + y p = [ ( A 0 t 2 + A 1 t ) cos ⁡ t + ( B 0 t 2 + B 1 t ) sin ⁡ t ] ″ + [ ( A 0 t 2 + A 1 t ) cos ⁡ t + ( B 0 t 2 + B 1 t ) sin ⁡ t ] = [ 2 A 0 cos ⁡ t + 2 ( 2 A 0 t + A 1 ) ( − sin ⁡ t ) + ( A 0 t 2 + A 1 t ) ( − cos ⁡ t ) ] + [ 2 B 0 sin ⁡ t + 2 ( 2 B 0 t + B 1 ) cos ⁡ t + ( B 0 t 2 + B 1 t ) ( − sin ⁡ t ) ] + [ ( A 0 t 2 + A 1 t ) cos ⁡ t + ( B 0 t 2 + B 1 t ) sin ⁡ t ] = [ 4 B 0 t + ( 2 A 0 + 2 B 1 ) ] cos ⁡ t + [ − 4 A 0 t + ( − 2 A 1 + 2 B 0 ) ] sin ⁡ t .

Comparing both sides, we have

4 B 0 = 1 2 A 0 + 2 B 1 = 0 − 4 A 0 = 0 − 2 A 1 + 2 B 0 = 0

which has the solution A 0 = 0 , A 1 = 1 / 4 , B 0 = 1 / 4 , B 1 = 0 . We then have a particular integral

y p = 1 4 t cos ⁡ t + 1 4 t 2 sin ⁡ t .

Example 2

Consider the following linear nonhomogeneous differential equation:

d y d x = y + e x .

This is like the first example above, except that the nonhomogeneous part ( e x ) is not linearly independent to the general solution of the homogeneous part ( c 1 e x ); as a result, we have to multiply our guess by a sufficiently large power of x to make it linearly independent.

Here our guess becomes:

y p = A x e x .

By substituting this function and its derivative into the differential equation, one can solve for A:

d d x ( A x e x ) = A x e x + e x A x e x + A e x = A x e x + e x A = 1.

So, the general solution to this differential equation is thus:

y = c 1 e x + x e x .

Example 3

Find the general solution of the equation:

d y d t = t 2 − y

f(t), t 2 , is a polynomial of degree 2, so we look for a solution using the same form,

y p = A t 2 + B t + C , where d y p d t = 2 A t + B

Plugging this particular integral with constants A, B, and C into the original equation yields,

2 A t + B = t 2 − ( A t 2 + B t + C ) , where t 2 − A t 2 = 0 and − B t = 2 A t and − C = B

Replacing resulting constants,

y p = t 2 − 2 t + 2

To solve for the general solution,

y = y p + y c

where y c is the homogeneous solution y c = c 1 e − t , therefore, the general solution is:

y = t 2 − 2 t + 2 + c 1 e − t