Paul's Online Math Notes
     
 
Online Notes / Differential Equations / Series Solutions to DE`s / Series Solutions
Differential Equations

You can navigate through this E-Book using the menu to the left. For E-Books that have a Chapter/Section organization each option in the menu to the left indicates a chapter and will open a menu showing the sections in that chapter. Alternatively, you can navigate to the next/previous section or chapter by clicking the links in the boxes at the very top and bottom of the material.

Also, depending upon the E-Book, it will be possible to download the complete E-Book, the chapter containing the current section and/or the current section. You can do this be clicking on the E-Book, Chapter, and/or the Section link provided below.

For those pages with mathematics on them you can, in most cases, enlarge the mathematics portion by clicking on the equation. Click the enlarged version to hide it.
   
Review : Taylor Series E-Book   Chapter   Section Euler Equations

 Series Solutions to Differential Equations

Before we get into finding series solutions to differential equations we need to determine when we can find series solutions to differential equations.  So, let’s start with the differential equation,

(1)

 

This time we really do mean nonconstant coefficients.  To this point we’ve only dealt with constant coefficients.  However, with series solutions we can now have nonconstant coefficient differential equations.  Also, in order to make the problems a little nicer we will be dealing only with polynomial coefficients.

 

Now, we say that x=x0 is an ordinary point if provided both

 

 

are analytic at x=x0.  That is to say that these two quantities have Taylor series around x=x0.  We are going to be only dealing with coefficients that are polynomials so this will be equivalent to saying that

 

 

 

for most of the problems.

 

If a point is not an ordinary point we call it a singular point. 

 

The basic idea to finding a series solution to a differential equation is to assume that we can write the solution as a power series in the form,

(2)

and then try to determine what the an’s need to be.  We will only be able to do this if the point x=x0, is an ordinary point.  We will usually say that (2) is a series solution around x=x0.

 

Let’s start with a very basic example of this.  In fact it will be so basic that we will have constant coefficients.  This will allow us to check that we get the correct solution.

 

Example 1  Determine a series solution for the following differential equation about .

                                                                 

Solution

Notice that in this case p(x)=1 and so every point is an ordinary point.   We will be looking for a solution in the form,

                                                             

 

We will need to plug this into our differential equation so we’ll need to find a couple of derivatives.

                             

 

Recall from the power series review section on power series that we can start these at n=0 if we need to, however it’s almost always best to start them where we have here.  If it turns out that it would have been easier to start them at n=0 we can easily fix that up when the time comes around.

 

So, plug these into our differential equation.  Doing this gives,

                                                 

 

The next step is to combine everything into a single series.  To do this requires that we get both series starting at the same point and that the exponent on the x be the same in both series. 

 

We will always start this by getting the exponent on the x to be the same.  It is usually best to get the exponent to be an n.  The second series already has the proper exponent and the first series will need to be shifted down by 2 in order to get the exponent up to an n.  If you don’t recall how to do this take a quick look at the first review section where we did several of these types of problems.

 

Shifting the first power series gives us,

                                            

 

Notice that in the process of the shift we also got both series starting at the same place.  This won’t always happen, but when it does we’ll take it.  We can now add up the two series.  This gives,

                                              

 

Now recalling the fact from the power series review section we know that if we have a power series that is zero for all x (as this is) then all the coefficients must have been zero to start with.  This gives us the following,

                                        

 

This is called the recurrence relation and notice that we included the values of n for which it must be true.  We will always want to include the values of n for which the recurrence relation is true since they won’t always start at n = 0 as it did in this case.

 

Now let’s recall what we were after in the first place.  We wanted to find a series solution to the differential equation.  In order to do this we needed to determine the values of the an’s.  We are almost to the point where we can do that.  The recurrence relation has two different an’s in it so we can’t just solve this for an and get a formula that will work for all n.  We can however, use this to determine what all but two of the an’s are.

 

To do this we first solve the recurrence relation for the an that has the largest subscript.  Doing this gives,

                                        

 

Now, at this point we just need to start plugging in some value of n and see what happens,

 

n=0

 

n=1

 

n=2

 

n=3

 

n=4

 

n=5

 

 

 

 

 

 

 

 

 

 

Notice that at each step we always plugged back in the previous answer so that when the subscript was even we could always write the an in terms of a0 and when the coefficient was odd we could always write the an in terms of a1.  Also notice that, in this case, we were able to find a general formula for an’s with even coefficients and an’s with odd coefficients.  This won’t always be possible to do.

 

There’s one more thing to notice here.  The formulas that we developed were only for k=1,2,… however, in this case again, the will also work for k=0.  Again, this is something that won’t always work, but does here.

 

Do not get excited about the fact that we don’t know what a0 and a1 are.  As you will see, we actually need these to be in the problem to get the correct solution. 

 

Now that we’ve got formulas for the an’s let’s get a solution.  The first thing that we’ll do is write out the solution with a couple of the an’s plugged in.

                  

 

The next step is to collect all the terms with the same coefficient in them and then factor out that coefficient.

              

 

In the last step we also used the fact that we knew what the general formula was to write both portions as a power series.  This is also our solution.  We are done.

 

Before working another problem let’s take a look at the solution to the previous example.  First, we started out by saying that we wanted a series solution of the form,

 

 

and we didn’t get that.  We got a solution that contained two different power series.  Also, each of the solutions had an unknown constant in them.  This is not a problem.  In fact, it’s what we want to have happen.  From our work with second order constant coefficient differential equations we know that the solution to the differential equation in the last example is,

 

 

Solutions to second order differential equations consist of two separate functions each with an unknown constant in front of them that are found by applying any initial conditions.  So, the form of our solution in the last example is exactly what we want to get.  Also recall that the following Taylor series,

 

 

 

Recalling these we very quickly see that what we got from the series solution method was exactly the solution we got from first principles, with the exception that the functions were the Taylor series for the actual functions instead of the actual functions themselves.

 

Now let’s work an example with nonconstant coefficients since that is where series solutions are most useful.

 

Example 2  Find a series solution around  for the following differential equation.

                                                                 

Solution

As with the first example p(x)=1 and so again for this differential equation every point is an ordinary point.  Now we’ll start this one out just as we did the first example.  Let’s write down the form of the solution and get its derivatives.

                 

 

Plugging into the differential equation gives,

                                                

 

Unlike the first example we first need to get all the coefficients moved into the series.

                                                

 

Now we will need to shift the first series down by 2 and the second series up by 1 to get both of the series in terms of xn.

                                           

 

Next we need to get the two series starting at the same value of n.  The only way to do that for this problem is to strip out the n=0 term.

                                 

 

We now need to set all the coefficients equal to zero.  We will need to be careful with this however.  The n=0 coefficient is in front of the series and the n=1,2,3… are all in the series.  So, setting coefficient equal to zero gives,

                                       

 

Solving the first as well as the recurrence relation gives,

                                             

 

Now we need to start plugging in values of n.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

There are a couple of things to note about these coefficients.  First, every third coefficient is zero.  Next, the formulas here are somewhat unpleasant and not all that easy to see the first time around.  Finally, these formulas will not work for k=0 unlike the first example.

 

Now, get the solution,

                  

 

Again, collect up the terms that contain the same coefficient, factor the coefficient out and write the results as a new series,

 

 

We couldn’t start our series at k=0 this time since the general term doesn’t hold for k=0.