
040108, 12:36 AM #1
 Posts
 4
solving ODE using power series
I'm trying to solve an ordinary differential equation, with nonconstant coefficients, using a
power series solution. The question is as follows:

solve the initial value problem:
hint: since the initial condition is given at , it is best to write the
solution as a series centered at .

I have attempted the question, but I got stuck in one of the steps. Here's what I have so far:
assume the solution can be written as a power series centered about x = 1. Then the solution and
its derivatives are:
where the are constant coefficients.
Substitute y, y', and y'' into the ODE in order to determine the :
now I need to shift indices to make all the exponents of (x1) equal to n, and also to try to
make all the summations start at n=0:
I don't know how to proceed from here. I have 4 summations. Two of them start at n=0 and has
the factor. The other two summations start at n=2 or n=1. I can't get them
to start at n=0 because if I try to do that, then I get:
and
As you can see, these have the factors and
instead of . So I don't know how to combine all this, since the summations
do not start at the same number.

I'm told that in the end, the answer to the coefficients should be:
and
So I just don't know how to get from the point I'm stuck to this answer. I'm also quite confused
as to how to solve the initial value problem, after obtaining the expected answer. The initial
conditions are: , . Using the initial condition y(1)
= 1, I get:
This does not make sense...what am I doing wrong?

040108, 12:54 AM #2
You have to consider the first few terms separately. If one sum starts at 0 and the other starts at 1, it just means that the first coefficient in the latter sum is zero. So if you only consider the coefficient in front of the zero order term in the whole expression, that will consist of only the term coming from (the first term of) the first sum.

040108, 01:25 AM #3
On the second thought, you don't even need to write out the first terms separately, because the expressions in this case are so that if you just extend the summation limit down to n=0, exactly the right terms become zero. For example,
because the factor n(n1) is zero when n=0 or n=1.
Your attempt to make the sums start at 0 is not exactly how it should be done in this case, because that makes the power of x inconsistent with the other sums.
I would also suggest first make the variable substitution z=x1 to make computation simpler.

040108, 06:05 AM #4
 Posts
 4
Thanks for your response. So it seems that I can just rewrite the two summations that are giving me trouble as starting at zero, without shifting the indices.
I'm still not sure how to solve the initial value problem. Since the initial value is x=1, when I try to substitute x=1 into the expression for y, the whole sum becomes 0

040108, 06:31 AM #5
 Posts
 16,194
edit: oh, that's what we're doing. Never mind  
Last edited by iceaura; 040108 at 12:20 PM.

040108, 07:58 AM #6

040108, 08:18 AM #7

040108, 09:39 AM #8
 Posts
 2,257
The whole sum does not become zero when x=1. When you write power series in the form you are implicitly assuming for all x including . In short, when it comes to power series. Without this shorthand one must write the power series in the form .
To solve this particular problem, the first thing to do is to make the substitution . This will simplify things a lot; hence the hint. The derivatives with respect to z are easy to compute via the chain rule:
. Thus the same differential equation applies with the substitution, but now you have
This simplifies even more (I'll let you work that out). Next write y as a power series in z about z=0 (this is the same power series as the series in x about x=1) and compute the derivatives. Next collect terms with equal powers so you have something that looks like
This can be true for all z if and only if each b_{n} is identically zero (Jess: Justify this step.) Setting the terms to zero should give you a recursive relationship of the form
I'll let you figure out what that offset is. With this, you should be able to show that
and
Finally, solve for a_{0} and a_{1} using the given initial conditions.
Similar Threads

By Forceman in forum ChemistryLast Post: 041308, 12:49 AMReplies: 41

By coberst in forum General PhilosophyLast Post: 033108, 11:58 AMReplies: 2

By weed_eater_guy in forum Computer Science & CultureLast Post: 030608, 12:52 AMReplies: 5

By Buffalo Roam in forum PoliticsLast Post: 090207, 09:42 PMReplies: 21

By ghost7584 in forum Religion ArchivesLast Post: 062005, 09:59 AMReplies: 73
Bookmarks