Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Homework Help: 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
    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.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook