Power series solutions to diffeq

In summary: And yes, Daniel, you're right. I was editing my post to make the necessary corrections before I saw your post.It appears that we have the same recurrence relation, but different series. Anyway, it's still the same problem, and we are both right about it.
  • #1
Townsend
232
0
I am suppose to find the power series solutions to some diffeqs.

[tex]y'=xy[/tex]

The method is to assume that
[tex]y=\sum_{n=0}^{\infty}{c_nx^n}[/tex]
is a solution to the diffeq.Then since we can differenitate term by term we have
[tex]
\frac{d}{dx} \left[ \sum_{n=0}^{\infty}{c_nx^n} \right]=\sum_{n=1}^{\infty}{nc_nx^{n-1}}
[/tex]

Then we sub these back into our differential equation to get
[tex]
\sum_{n=1}^{\infty}{c_{n}nx^{n-1}}-x\sum_{n=0}^{\infty}{c_nx^n}=
[/tex]

[tex]
\sum_{n=1}^{\infty}{c_{n}nx^{n-1}}-\sum_{n=0}^{\infty}{c_nx^{n+1}}=0
[/tex]
So now I need to try to get the powers of [tex]x^n[/tex] and the starting index the same for both sums. But I cannot seems to figure how to push things into that form.

If I try adding one to the index of the left most series so I have the same starting index I get to
[tex]
\sum_{n=0}^{\infty}c_{n+1}(n+1)x^{n-1} -\sum_{n=0}^{\infty}{c_nx^{n+1}}=0
[/tex]

But without having the same powers for x I cannot factor all the x's. How will I get my recurrance relation then?

Thanks for the help

Jeremy
 
Last edited:
Physics news on Phys.org
  • #2
You'll have to write a few terms of one of the series seperately out front and then start the sum. Actually, the usual method has a sum over xn+s, where n goes from 0 to infinity and s is the starting power. s would be determined by these terms out front. If this problem was given to you and s wasn't mentioned, it should probably work out by itself.
 
  • #3
It's very simple.

[tex] y= C \sum_{n=0}^{+\infty} c_{n}x^{n} [/tex](1)

[tex] \frac{dy}{dx}=C\sum_{n=0}^{+\infty} nc_{n}x^{n-1} [/tex] (2)

[tex] xy=C\sum_{n=0}^{+\infty} c_{n}x^{n+1} [/tex] (3)

Equate (2) & (3) and write the term

[tex] c_{1}+2c_{2}x+3c_{3}x^{2}+...=c_{0}x+c_{1}x^{2}+c_{2}x^{3}+... [/tex] (4)

In order for the 2 sides of (4) to be equal,one gets

[tex] c_{1}=0 [/tex] (5)

[tex] c_{n}=(n+2)c_{n+2} [/tex] (6)

It's preferable to choose [tex] c_{0}=1 [/tex] (why?,because i put a C (arbitrary constant) in front of the series,so [itex] y(0)=Cc_{0}[/itex],i.e.a product of arbitrary constants).

So the series (solution) becomes

[tex] \left\{\begin{array}{c}y(x)=C \sum_{0}^{+\infty} c_{n}x^{n}\\c_{n}=(n+2)c_{n+2}\\c_{0}=1\\c_{1}=0 \end{array}\right [/tex] (5)



Daniel.
 
  • #4
(5) is the Taylor expansion around 0 (okay,some of you may know it a MacLaurin expansion) of the solution of the ODE

[tex] \frac{dy}{dx}=xy [/tex] (6)

,i.e.it is the Taylor expansion around 0 of the function

[tex] y(x)=Ce^{\frac{x^{2}}{2}} [/tex] (7)

Daniel.

P.S.Note the C in (7) is the same with the C in (1-3) and (5).
 
  • #5
Jeremy:

We seem to be using the same methods in our class as you are. I think a general rule of thumb is this: you cannot combine two power series (by putting the two general terms of each under the same summation sign) unless the degree of the two series is the same, which requires both that x is raised to the same power (in terms of n), and that the starting index is the same in both cases. However, if the first condition is satisfied, but the starting indices are not the same, then it is not a big deal because you can just remove terms from (take them outside) series with lower starting indices until they match that of the series with the highest one, just as StatusX said. However, if the first condition is not satisfied (one series is in terms of x^n, and the other is in terms of, say, x^n+1), then you're going to have much more trouble. So your efforts should focus first on having both series with x raised to the same power (worry about starting indices later). That having been said, here's how we usually dealt with the y' terms:

[tex]\frac{d}{dx} \left[ \sum_{n=0}^{\infty}{c_nx^n} \right]=\sum_{n=1}^{\infty}{nc_nx^{n-1}} = \sum_{n=0}^{\infty}{(n+1)c_{n+1}x^{n}} [/tex]

Note that last shift.

Then the DE becomes:

[tex] \sum_{n=0}^{\infty}{(n+1)c_{n+1}x^{n}} - \sum_{n=0}^{\infty}{c_{n}x^{n+1}} = 0 [/tex]

drat! Now we have an x^n+1 because of that pesky x multiplier:

[tex] \sum_{n=0}^{\infty}{(n+1)c_{n+1}x^{n}} - \sum_{n=1}^{\infty}{c_{n-1}x^{n}} = 0 [/tex]

It's all good...except starting indices are not the same. Throw out a term from the left series:

[tex] c_1 + \sum_{n=1}^{\infty}{(n+1)c_{n+1}x^{n}} - \sum_{n=1}^{\infty}{c_{n-1}x^{n}} = 0 [/tex]

[tex] c_1 + \sum_{n=1}^{\infty}{[(n+1)c_{n+1}x^{n} - c_{n-1}x^{n}]} = 0 [/tex]

Which, if I remember right, tells you that c1 = 0 (don't quote me on that), and that the rest of it is zero as well:

[tex] \sum_{n=1}^{\infty}{[(n+1)c_{n+1} - c_{n-1}]x^{n}} = 0 [/tex]

[tex] (n+1)c_{n+1} - c_{n-1} = 0 [/tex]

[tex] c_{n+1} = \frac{c_{n-1}}{n+1} [/tex]

[tex] n = 1, 2, 3... [/tex]

This is the same as:

[tex] c_{m+2} = \frac{c_{m}}{m+2} [/tex]

[tex] m = 0, 1, 2... [/tex]

I hope that is the correct recurrence relation. Although I can sort of see that this might result in a Taylor series of the function Daniel came up with, I have no idea how he just stated the solution by inspection. Anyways, I thought the point was to actually come up with the power series solution. Can you take it from here?
 
Last edited:
  • #6
Cepheid,if I'm not mistaking,your solution is

[tex] y(x)^{Cepheid}=c_{1}x+c_{3}x^{3}+... [/tex]

which is not okay,because it is not the Taylor expansion of [itex] e^{x^{2}/2} [/itex],moreover,it doesn't verify identically the ODE,not to mention that y(0)=0,which is incorrect,as it wasn't a Cauchy problem to begin with...

Daniel.
 
Last edited:
  • #7
One mistake I made that I can see is that when you "throw out" the n=0 term from the first series, the term you throw out is actually c1, not c0. (just plug in n=0 to check). So only the even terms of my series survive, including c0.

I will edit that post to make the correction.
 
  • #8
cepheid said:
Jeremy:

We seem to be using the same methods in our class as you are. I think a general rule of thumb is this: you cannot combine two power series (by putting the two general terms of each under the same summation sign) unless the degree of the two series is the same, which requires both that x is raised to the same power (in terms of n), and that the starting index is the same in both cases. However, if the first condition is satisfied, but the starting indices are not the same, then it is not a big deal because you can just remove terms from (take them outside) series with lower starting indices until they match that of the series with the highest one, just as StatusX said. However, if the first condition is not satisfied (one series is in terms of x^n, and the other is in terms of, say, x^n+1), then you're going to have much more trouble. So your efforts should focus first on having both series with x raised to the same power (worry about starting indices later). That having been said, here's how we usually dealt with the y' terms:

[tex]\frac{d}{dx} \left[ \sum_{n=0}^{\infty}{c_nx^n} \right]=\sum_{n=1}^{\infty}{nc_nx^{n-1}} = \sum_{n=0}^{\infty}{(n+1)c_{n+1}x^{n}} [/tex]

Note that last shift.

Then the DE becomes:

[tex] \sum_{n=0}^{\infty}{(n+1)c_{n+1}x^{n}} - \sum_{n=0}^{\infty}{c_{n}x^{n+1}} = 0 [/tex]

drat! Now we have an x^n+1 because of that pesky x multiplier:

[tex] \sum_{n=0}^{\infty}{(n+1)c_{n+1}x^{n}} - \sum_{n=1}^{\infty}{c_{n-1}x^{n}} = 0 [/tex]

It's all good...except starting indices are not the same. Throw out a term from the left series:

[tex] c_0 + \sum_{n=1}^{\infty}{(n+1)c_{n+1}x^{n}} - \sum_{n=1}^{\infty}{c_{n-1}x^{n}} = 0 [/tex]

[tex] c_0 + \sum_{n=1}^{\infty}{[(n+1)c_{n+1}x^{n} - c_{n-1}x^{n}]} = 0 [/tex]

Which, if I remember right, tells you that c0 = 0 (don't quote me on that), and that the rest of it is zero as well:

[tex] \sum_{n=1}^{\infty}{[(n+1)c_{n+1} - c_{n-1}]x^{n}} = 0 [/tex]

[tex] (n+1)c_{n+1} - c_{n-1} = 0 [/tex]

[tex] c_{n+1} = \frac{c_{n-1}}{n+1} [/tex]

[tex] n = 1, 2, 3... [/tex]

This is the same as:

[tex] c_{m+2} = \frac{c_{m}}{m+2} [/tex]

[tex] m = 0, 1, 2... [/tex]

I hope that is the correct recurrence relation. Although I can sort of see that this might result in a Taylor series of the function Daniel came up with, I have no idea how he just stated the solution by inspection. Anyways, I thought the point was to actually come up with the power series solution. Can you take it from here?


Yes, that is what I was looking for. It would be easy to solve since it is seperable but I really wanted to work out the recurrsion relation the way the book suggest to do it and I could not understand how to do this until now.

As far as seeing the solution Daniel came up with all you would have to know is that the power series solution for the exponential function is
[tex]
e^x=\sum_{n=0}^{\infty} \frac{x^n}{n!}
[/tex]

Since the solution is
[tex]
y=C_0e^{\frac{x^2}{2}}
[/tex]

All you would have to do is push this into the series given above and you're there.

[tex]
C_0 \sum_{n=0}^{\infty} \frac{x^{2n}}{(2n)!}
[/tex]

And after solving the recurrance relation you gave me I found
[tex]
c_{2n}= \frac{c_0}{(2n)!}
[/tex]
So our solution checks out perfectly.

But seriously thank you for your sol. it was just what I was hoping for.

Jeremy
 
  • #9
Nope.

[tex] \frac{x^{2n}}{(2n)!}\neq \frac{ x^{2n}}{2^{n}n!} [/tex]

The RHS correctly describes the general term of the Taylor expansion of [itex] e^{x^{2}/2} [/itex] around 0.


Daniel.
 
  • #10
dextercioby said:
Nope.

[tex] \frac{x^{2n}}{(2n)!}\neq \frac{ x^{2n}}{2^{n}n!} [/tex]

The RHS correctly describes the general term of the Taylor expansion of [itex] e^{x^{2}/2} [/itex] around 0.


Daniel.

You're right...

I think I am getting too tired to work on this any more tonight. Thanks for the help. I will look at this in the morning.

Jeremy
 
  • #11
Jeremy

From the recurrence relation I derived:

[tex] c_2 = \frac{c_0}{2} [/tex]

[tex] c_4 = \frac{c_2}{4} = \frac{c_0}{2 \cdot 4} [/tex]

[tex] c_6 = \frac{c_4}{6} = \frac{c_0}{2 \cdot 4 \cdot 6} [/tex]

[tex] c_8 = \frac{c_6}{8} = \frac{c_0}{2 \cdot 4 \cdot 6 \cdot 8} [/tex]

As you can see, the denominator is not (2n)!, for that would be 8*7*6*...*2*1 in that last example. Instead we have:

[tex] c_{2n} = \frac{c_0}{2n \cdot (2n-2) \cdot ... \cdot 2} [/tex]

But that denominator can be expressed as a factorial. Hint: what can you factor out from each of the factors in that product? So how many times does it appear? To answer that, you just need to know: how many factors are in that product?

Hint that will give it away: In the end you will of course end up with the result dexter gave.
 

FAQ: Power series solutions to diffeq

What is a power series solution to a differential equation?

A power series solution to a differential equation is a representation of a function as an infinite sum of terms, each of which is a power of the independent variable. It is a useful method for solving differential equations that cannot be solved using traditional algebraic methods.

How do you determine the convergence of a power series solution?

The convergence of a power series solution can be determined by using the ratio test or the root test. These tests evaluate the limit of the absolute value of the ratio or root of consecutive terms in the series. If the limit is less than 1, the series converges.

What are the advantages of using a power series solution?

One advantage of using a power series solution is that it can be used to approximate solutions to differential equations with high accuracy. Additionally, it can be used to solve nonlinear differential equations that cannot be solved using traditional methods.

Are there any limitations to using a power series solution?

Yes, there are limitations to using a power series solution. It can only be used for certain types of differential equations, such as those with constant coefficients. It is also not suitable for solving differential equations with singularities or discontinuities.

How do you find the coefficients in a power series solution?

The coefficients in a power series solution can be found by using the Taylor series expansion of the function. This involves taking derivatives of the function and evaluating them at a specific point, typically the point where the series is centered. The coefficients can also be found by using recurrence relations or by equating coefficients in the original differential equation.

Similar threads

Back
Top