In talking about power series in a previous post, I mentioned one of their uses: as an aid in solving differential equations. This reminds me of a neat trick for solving some differential equations, which I will discuss in this post.
A standard method for solving linear differential equations with constant coefficients is to assume a trial solution of the form , run it through the differential equation, solve the resulting algebraic equation for r, and then take it from there.
For example, consider the differential equation
Let’s suppose that there is a solution to the differential equation of the form . Then
Inserting the expressions for y and into the differential equation, we get
which means that
Since the exponential factor is not equal to zero, it follows that
which is equivalent to
The conclusion is that a function of the form really is a solution to the differential equation, but only if r = –1 or r = 1. So we have found two solutions to the differential equation,
The next step in this standard solution is to invoke some theory: There is a theorem that states that because the differntial equation is second order, linear, and has constant coefficients, once we have found two linearly independent solutions, which we have done, then all possible solutions of the differential equation are of the form
where c1 and c2 are constants that are typically determined in each situation using initial conditions.
There is a nice linear algebra connection here, which dovetails nicely with our recent series of posts on the concept of a basis: The set of solutions of the differential equation is a two-dimensional vector space, and the two fundamental solutions (as they’re called in many textbooks), and , form a basis for the solution space.
An alternative method
Now there is something annoying about the standard solution method; at least it annoys some people. Even though the standard method amounts to the “guess and check” method, which is quite respectable, some people consider it an unmotivated “rabbit-out-of-the-hat” trick. They would prefer a more deductive procedure, where no assumptions about the solution have to be made.
OK, here goes. Adopt a more abstract approach, and write the differential equation in the form
In the previous line, D stands for the operation of differentiation. That is, Dy represents the operation of differentiating y, which results in . Similarly, D2 means to apply the D-operator twice, so that D2y results in . Thought of as operators, real numbers applied to functions result in multiplication of the function by the number. (One really ought to note that there is an identity operator present, as in , but it’s customary to omit it, and hope that readers somehow understand the implicit meaning.)
Now for the trick: Treat the operator expression as if it were an algebraic expression, and factor it as such. Then the differential equation takes the form
Now let z stand for the function (D + 1)y. Notice that we can thereby reduce the second order differential equation to a pair of first-order equations:
Being first-order equations, they will be easier to solve; what we’re banking on is that solving two first-order equations is easier than solving one second-order equation. The savings may not be much in this example, but it’s a general idea worth keeping in mind.
The first equation, , can be translated back into usual form as
This is a linear first order equation, and can be solved by multiplying each term by the integrating factor :
(where c is a constant)
Now take this expression, substitute it for z in the other first-order equation, and solve for y:
(multiplying each term by the integrating factor )
(multiplying each term by )
And this is exactly the same as the solution obtained by the standard method, if we identify the arbitrary constants c1 with b and c2 with c/2.
It’s worth mentioning that D is a linear operator, and from this it follows that the operator in the differential equation, D2 – 1I, is also a linear operator. Solving the differential equation, in our example, amounts to finding the null-space (or kernel, if you’d rather call it that) of the linear operator, which is itself a subspace of the vector space that D2 – 1I acts on. We never did specify what that space is, did we? But it lends support to the comment made earlier that the solution of this differential equation is a two-dimensional vector space, so all one need do is find a basis for the space; then all possible solutions are linear combinations of the two basic solutions.
The operator approach works well for linear ordinary differential equations with constant coefficients. Of course, the resulting operator (i.e, the polynomial in D) must be factorable. But the method still works if the factors are complex. For example, consider the differential equation
which can be written in operator form as
and in factored form as
Apply the operator method and see if you can solve this equation.
* * *
Here’s the result: The space of all possible solutions consists of functions of the form
where c1 and c2 are constants. Typically we would like to extract real solutions from the complex expression, and by using Euler’s formula, you can write the solution in the form
where A and B are constants.
* * *
Click here to read An Operator Method for Solving Second Order Differential Equations, Part 2.