1. Not finding help here? Sign up for a free 30min tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Series solutions of differential equations

  1. May 2, 2008 #1
    i haven't been able to solve this question in my textbook:

    --> Find a power series solution for y''(x) - y(x)=0

    First i attempted to find singular and ordinary points for the differential equation and i found that x=0 is an ordinary point. then i set y=[tex]\sum[/tex]A(n)x^n for n=0 to infinity
    subsituted into the differential equation after finding the second derivative.
    then after some manipulations i found that
    y=A(0)[tex]\sum[/tex]x^(2n)/(2n)! + A(1)[tex]\sum[/tex]x^(2n+1)/(2n+1)!

    the solution says it is equal to e^x+e^(-x)

    I tried using the series of e^x=[tex]\sum[/tex]x^n/n! but i couldnt get the solution

    help people

  2. jcsd
  3. May 2, 2008 #2


    User Avatar

    Are there any initial conditions for this problem (otherwise there is no way you can get that solution from the series).

    I get the following: go through the first steps (derivatives, etc...)
    get the recurrence relation of A(n+2)=A(n)/[(n+2)(n+1)], solving the recurrence relation I get A(2k)=A(0)/k! and A(2k+1)=A(1)/k!

    [itex]y(x)=\sum{A(2k)x^{2k}+A(2k+1)x^{2k+1}} = \sum{\frac{A(0)x^{2k}}{k!}+\frac{A(1)x^{2k+1}}{k!}}[/itex]

    Which is certainly going to be 2 functions of the e^(x) type. However, without initial conditions you can't get A(0) and A(1), although it's possible to see (given that you know the final solution) that A(0) probably = 1 and A(1) probably = -1 (since that -1 can be changed into (-x)^(2k+1) preserving the -1 since the power is always odd).

    Note that the summation is still from 0 to infinity.
    Last edited: May 2, 2008
  4. May 2, 2008 #3
    exk there are no initial conditions the question asks for a general solution, i'm sorry for not specifying that before.
    the question is basically trying to show that using characteristic equation and a power series both give the same answer
  5. May 2, 2008 #4
    but how do you go from [itex]y(x)=\sum{A(2k)x^{2k}+A(2k+1)x^{2k+1}} = \sum{\frac{A(0)x^{2k}}{k!}+\frac{A(1)x^{2k+1}}{k!} }[/itex] to an e^x and e^(-x) disregarding the coefficients of the two functions?
  6. May 2, 2008 #5


    User Avatar
    Science Advisor
    Homework Helper

    If the problem asks for a general solution and doesn't give initial conditions then e^(x)+e^(-x) isn't the general solution. A*e^x+B*e^(-x) is as you know from the characteristic equation. So the series solution should also have two undetermined constants in it.
  7. May 2, 2008 #6


    User Avatar

    Erm, I am sorry I was in too much of a hurry before apparently.
    You solved the recurrence correctly getting:
    y(x)=\sum{A(2k)x^{2k}+A(2k+1)x^{2k+1}} = \sum{\frac{A(0)x^{2k}}{(2k)!}+\frac{A(1)x^{2k+1}}{(2k+1)!} }

    With a little manipulation I get:


    Now we can expand both summations to get:

    This can be re-written as:

    Now suppose you know before hand that the solution should be e^x+e^-x (lets assume the undetermined coefficients are 1s for now).

    ok now expand that solution as a series (the only difference between the two being that e^-x produces a negative term for odd powers of x). Add them together to get:

    rewrite that as [itex]2+2\sum_{1}^{\infty}{\frac{x^{2k}}{(2k)!}}[/itex] which looks very similar to our y(x) above if A(0)=2 and A(1)=0.

    However your general solution is A*e^x+B*e^-x as Dick pointed out. So going back to addition of the two expansions you get:
    Now recall the original solution progress:
    Re-write as:

    Finally we can see that if A(0)=A+B and A(1)=A-B you get:

    This doesn't completely agree since 2A=A+B iff A=B, but lets say A could be = to B. Then you are back with A(0)=2*A and A(1)=0, which corresponds to what I showed above.

    I hope I didn't make a mistake again, all these summations are confusing. However, I certainly have no idea what would prompt me to do the above steps without actually knowing the end solution. Personally I would have probably left the solution in terms of summations.
  8. May 3, 2008 #7
    the sums all look correct and are quite similar to what i have on paper but i'm still stuck from going from the sums to e^x
    i have no clue on how to manipulate the sums to get it to e^x by using the e^x series (e^x= [tex]\sum[/tex]x^n/n!)....or by any other method
  9. May 3, 2008 #8


    User Avatar

    Take a second look at what I did above. The conversion to e^x can be seen fairly clearly in tex lines 3-5. The actual solution conversion is seen in this:
    y(x)=(A+B)+(A+B)\sum_{1}^{\infty}{\frac{x^{2k}}{(2 k)!}}+(A-B)\sum_{1}^{\infty}{\frac{x^{2k+1}}{(2k+1)!}}
  10. May 3, 2008 #9


    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    Even if you don't see the conversion to exponentials... you should recognize the power series
    [tex]\sum_{x =0}^{+\infty} \frac{x^{2n}}{(2n)!}[/tex]
    as the Taylor series for familiar function.
  11. May 3, 2008 #10


    User Avatar
    Staff Emeritus
    Science Advisor

    I'm not sure that many students are really familiar with that particular function.
  12. May 3, 2008 #11


    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    They should at least recognize the other function it's similar to -- and then they could derive the needed change of variable.
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?

Similar Discussions: Series solutions of differential equations