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

and

#### 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,

and

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 *c*_{1} and *c*_{2} 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, *D*^{2} means to apply the *D*-operator twice, so that *D*^{2}*y* 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

or simply

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:

and

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 *c*_{1} with *b* and *c*_{2} 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, *D*^{2} – 1*I*, 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 *D*^{2} – 1*I* 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 *c*_{1} and *c*_{2} 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**.

Pingback: An Operator Method for Solving Second Order Differential Equations, Part 2 | QED Insight

beautifully explained, i now understand the D-operator method more it was taught in a lecture

Glad to hear the post was useful! All the best wishes in your studies!

Awesome!

thanks a lot sir…. helped me in deriving the equation of simple harmonic motion 🙂

My pleasure, glad it was helpful!

This is brilliant, it’s easy to understand and follow. Thank you so much.

Reblogged this on Math Online Tom Circle.

Reblogged this on Down the Quantum Hole.

Excellent explanation

Thank you soooo much!!