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

Bessel function and Bessel D.E.

  1. Mar 31, 2005 #1

    quasar987

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    I'm trying to show that the Bessel function of the first kind satisfies the Bessel differential equation for m greater of equal to 1.

    The Bessel function of the first kind of order m is defined by

    [tex]J_m(x) = \sum_{n=0}^{\infty} \frac{(-1)^n}{2^{m+2n}n!(n+m)!}x^{m+2n} = x^m \sum_{n=0}^{\infty} \frac{(-1)^n}{2^{m+2n}n!(n+m)!}(x^2)^n[/tex]

    and suposedly it satisfies the Bessel differential equation, which would then write

    [tex]x^2\frac{d^2 J}{dx^2} + x\frac{dJ}{dx} + (x^2-m^2)J = 0[/tex]

    I'm trying to help myself in doing that by following the steps laid down in the mathworld page on this subject. There is one detail however, it is that on the site, they proceed "backward". I.e. they start by saying that the Bessel function is defined by the solution to the Bessel D.E. and then proceed to find its form.

    My main trouble with the steps they take is that in equation (30), they state that a_1 = 0. But a_n, as they will find at the end, is

    [tex]a_n = \frac{(-1)^n}{2^{m+2n}n!(n+m)!}[/tex]

    and of course, a_1 is not 0 ! So what the greasy poop's going on here? :grumpy:
     
  2. jcsd
  3. Mar 31, 2005 #2

    dextercioby

    User Avatar
    Science Advisor
    Homework Helper

    I didn't see anything devious.You have to be careful with that change of summation variable from "n" to "l"...

    Daniel.

    P.S.If u still think there's something wrong,i'm sure u can find this matter in any ODE book treating Bessel's equation.
     
  4. Mar 31, 2005 #3

    quasar987

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    I didn't really get into any of the change of variable from n to l thing. My main concern is that they assert that a_1 is 0 but it's not.
     
  5. Apr 1, 2005 #4

    dextercioby

    User Avatar
    Science Advisor
    Homework Helper

    It is zero.Again,follow my advice and look into an ODE book which treats Bessel eq.(and implicitely the functions).

    Daniel.
     
  6. Apr 1, 2005 #5

    quasar987

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    I'm very far from having access to books right now. What's wrong with...

    [tex]a_1 = \frac{(-1)^1}{2^{m+2(1)}1!(1+m)!} = \frac{-1}{2^{m+2}(1+m)!} [/tex]

    ?
     
  7. Apr 1, 2005 #6

    dextercioby

    User Avatar
    Science Advisor
    Homework Helper

    There you go.It can't get any more clear than that:



    Daniel.
     

    Attached Files:

  8. Apr 1, 2005 #7

    quasar987

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    thanks dexter but this link really doesn't say anything mathworld doesn't already say and so I'm still oblivious to how we can simultaneously have a_1 = a_3 = ... = 0 and

    [tex]a_n = \frac{(-1)^n}{2^{m+2n}n!(n+m)!}[/tex]

    since it would imply

    [tex]a_1= \frac{-1}{2^{m+2}(1+m)!}\neq 0[/tex]


    (Unless of course, a_n is defined in this way only for n even, but then it would be false to write the Bessel function as

    [tex]J_m(x) = \sum_{n=0}^{\infty} \frac{(-1)^n}{2^{m+2n}n!(n+m)!}x^{m+2n}[/tex]

    with no mention of the crucial fact that n only takes even values.)
     
    Last edited: Apr 1, 2005
  9. Apr 1, 2005 #8

    dextercioby

    User Avatar
    Science Advisor
    Homework Helper

    Well,if u didn't understand how it was done there,then it's not my fault,nor the material's.I can't do anything more...

    Daniel.
     
  10. Apr 1, 2005 #9

    quasar987

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    You could say something like

    We don't, because ....
     
  11. Apr 1, 2005 #10
    [tex]J_m(x) = \sum_{n=0}^{\infty} \frac{(-1)^n}{2^{m+2n}n!(n+m)!}x^{m+2n}[/tex]

    [tex]= x^m\sum_{n=0}^{\infty} b_n x^{2n}[/tex]

    the terms with odd powers of [itex]x[/itex] aren't included in the series above.

    And I don't see where mathworld or dexter's document says

    [tex]a_n = \frac{(-1)^n}{2^{m+2n}n!(n+m)!}[/tex]
     
    Last edited: Apr 1, 2005
  12. Apr 1, 2005 #11

    shmoe

    User Avatar
    Science Advisor
    Homework Helper

    The series form you've written here does no have the same [tex]a_n[/tex] as you started with (in mathworlds (4)). Notice how the powers of this series are limited to m, m+2, m+4, ...., it's built in to this notation that 'every other' term drops out.

    Look carefully at (36)-(40), especially (36) which shows how they've seperated the even/odd terms.
     
  13. Apr 1, 2005 #12

    quasar987

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    Right. The a_n of (4) is not the same as the a_n I assumed because in (4), x is to the n, wherehas in the final form of the Bessel function, x is to the 2n.

    Ok, so this whole page from mathworld won't help me at all to show that the Bessel function of the first kind satisfies the Bessel differential equation for m greater of equal to 1, will it?

    Anyone sees how this can be done?

    I've tried playing with the following recurrence relations (about halfway on the page), but without success :grumpy:.
     
  14. Apr 1, 2005 #13

    dextercioby

    User Avatar
    Science Advisor
    Homework Helper

    Solving the equation and finding the initial solution (given by the problem) is an equivalent approach.

    Daniel.
     
  15. Apr 1, 2005 #14
    Why don't you just substitute it into the equation and show that it satisfies it? There's a bit of algebra, but it's straightforward and not that long anyhow.
     
  16. Apr 1, 2005 #15

    quasar987

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    It doesn't add up for me. What I did is set

    [tex]J_p(x) = x^p f(x)[/tex]

    Where

    [tex]f(x) = \sum_{n=0}^{\infty}a_n X^n[/tex]

    where

    [tex]a_n = \frac{(-1)^n}{2^{p+2n}n!(n+p)!}[/tex]

    and

    [tex]X = X(x) = x^2[/tex]

    That way, f(x) is a power serie with a radius of convergence [itex]R \neq 0[/itex], so according to a theorem seen in class,

    [tex]f ' (x) = \sum_{n=0}^{\infty} \frac{d}{dx}a_n x^{2n} = \sum_{n=0}^{\infty}2n a_n x^{2n-1}[/tex]

    and

    [tex]f '' (x) = \sum_{n=0}^{\infty}2n(2n-1) a_n x^{2n-2}[/tex]

    With this I construct [itex]x^2 J''(x)[/itex], [itex]x J'(x)[/itex] and [itex] (x^2-p^2)J(x)[/itex], and add them up. However, after double checking all my calculation, I am left with

    [tex] \sum_{n=0}^{\infty}(4np+4n^2+x^2) a_n x^{2n+p}[/tex]

    Is there anything wrong with what I did above?

    How did you do it Data?
     
  17. Apr 1, 2005 #16
    Well, if you're going to pull the [itex]x^p[/itex] out of the sum, then you'll need to use product rule later to evaluate [itex]J''[/itex]. Anyways, here's how I'd go about it:

    We have

    [tex]J_m(x) = \sum_{n=0}^{\infty} \frac{(-1)^n}{2^{m+2n}n!(n+m)!}x^{m+2n}, \; m \geq 0[/tex]

    which is a power series with radius of convergence [itex]\infty[/itex], so differentiating termwise is justified. Then

    [tex]J^\prime_m(x) = \sum_{n=0}^{\infty} \frac{(-1)^n(m+2n)}{2^{m+2n}n!(n+m)!}x^{m+2n-1}[/tex]

    [tex]J^{\prime\prime}_m(x) = \sum_{n=0}^{\infty} \frac{(-1)^n(m+2n)(m+2n-1)}{2^{m+2n}n!(n+m)!}x^{m+2n-2}[/tex]

    Our differential equation is

    [tex] x^2\frac{d^2\psi}{dx^2} + x\frac{d\psi}{dx} + (x^2-m^2)\psi = 0[/tex]

    so substituting [itex]J[/itex] for [itex]\psi[/itex] on the left side gives

    [tex] x^2\frac{d^2J}{dx^2} + x\frac{dJ}{dx} + (x^2-m^2)J[/tex]

    [tex] = x^2\sum_{n=0}^{\infty} \left[\frac{(-1)^n(m+2n)(m+2n-1)x^{m+2n-2}}{2^{m+2n}n!(n+m)!}\right] + x\sum_{n=0}^{\infty} \left[\frac{(-1)^n(m+2n)x^{m+2n-1}}{2^{m+2n}n!(n+m)!}\right] + (x^2 - m^2)\sum_{n=0}^{\infty} \left[\frac{(-1)^nx^{m+2n}}{2^{m+2n}n!(n+m)!}\right] [/tex]

    [tex] = \sum_{n=0}^{\infty} \left[\frac{(-1)^n(m+2n)(m+2n-1)x^{m+2n}}{2^{m+2n}n!(n+m)!}\right] + \sum_{n=0}^{\infty} \left[\frac{(-1)^n(m+2n)x^{m+2n}}{2^{m+2n}n!(n+m)!}\right] + \sum_{n=0}^{\infty} \left[\frac{(-1)^nx^{m+2n+2}}{2^{m+2n}n!(n+m)!}\right] - \sum_{n=0}^{\infty} \left[\frac{(-1)^nm^2x^{m+2n}}{2^{m+2n}n!(n+m)!}\right] [/tex]

    [tex] = \sum_{n=0}^\infty \left[\frac{(-1)^nx^{m+2n}}{2^{m+2n}n!(n+m)!}\left((m+2n)(m+2n-1)+(m+2n)-m^2\right)\right] + \sum_{n=0}^{\infty} \left[\frac{(-1)^nx^{m+2n+2}}{2^{m+2n}n!(n+m)!}\right][/tex]

    [tex] = \sum_{n=0}^\infty \left[\frac{(-1)^nx^{m+2n}}{2^{m+2n}n!(n+m)!}\left(4n^2 + 4nm)\right)\right] + \sum_{n=0}^{\infty} \left[\frac{(-1)^nx^{m+2n+2}}{2^{m+2n}n!(n+m)!}\right][/tex]

    by theorems, this is equal to zero if and only if the coefficient of each power of [itex]x[/itex] is zero. The smallest power of [itex]x[/itex] on the left side is [itex]m[/itex] when n=0 in the first series and its coefficient is (from the first, second, and last series) just

    [tex] \frac{(-1)^0(4(0)^2+4(0)m)}{2^m0!m!} = 0[/tex]

    as desired. We can thus ignore the first term in the first series and make a change of index, [itex] n \rightarrow n+1[/itex]. This gives

    [tex]\sum_{n=0}^\infty \left[\frac{(-1)^{n+1}x^{m+2n+2}}{2^{m+2n+2}(n+1)!(n+m+1)!}\left(4(n+1)^2 + 4(n+1)m)\right)\right] + \sum_{n=0}^{\infty} \left[\frac{(-1)^nx^{m+2n+2}}{2^{m+2n}n!(n+m)!}\right][/tex]

    [tex] =\sum_{n=0}^\infty \left(\frac{(-1)^{n+1}}{2^{m+2n+2}(n+1)!(n+m+1)!}\left(4(n+1)^2 + 4(n+1)m\right) + \frac{(-1)^n}{2^{m+2n}n!(n+m)!}\right)x^{m+2n+2}[/tex]

    [tex] =\sum_{n=0}^\infty \left(\frac{(-1)^{n+1}}{2^{m+2n+2}(n+1)!(n+m+1)!}\left(4n^2 + 4nm + 8n + 4m + 4\right) + \frac{(-1)^n}{2^{m+2n}n!(n+m)!}\right)x^{m+2n+2}[/tex]

    and again we look at coefficients. This time let's set

    [tex]k_n = \frac{(-1)^{n+1}}{2^{m+2n+2}(n+1)!(n+m+1)!}[/tex]

    so our sum reduces to just

    [tex] \sum_{n=0}^\infty k_n\left(\left(4n^2+4nm+8n+4m+4\right) + \left((-1)2^2(n+1)(n+m+1)\right)\right)x^{m+2n+2}[/tex]

    and thus the coefficient of [itex]x^{m+2n+2}[/itex] is

    [tex]k_n\left(\left( 4n^2 + 4nm + 8n + 4\right) + \left((-1)2^2(n+1)(n+m+1)\right)\right)[/tex]

    [tex] = k_n \left(4n^2 + 4nm + 8n +4m+ 4 - 4(n^2 + 2n + nm + m + 1)\right)[/tex]

    [tex] = k_n\left( (4n^2 - 4n^2) + (4nm-4nm) + (8n-8n) + (4m-4m) + (4-4) \right)[/tex]

    [tex] = 0[/tex]

    and thus every coefficient in the series is [itex]0[/itex], and we have found

    [tex] x^2\frac{d^2J}{dx^2} + x\frac{dJ}{dx} + (x^2-m^2)J = 0[/tex]

    as we wanted.
     
    Last edited: Apr 1, 2005
  18. Apr 1, 2005 #17

    quasar987

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    Phew! Thanks for the long and most complete reply ever :)
     
  19. Apr 1, 2005 #18
    Like I said, just a little arithmetic :wink:
     
  20. Apr 1, 2005 #19

    saltydog

    User Avatar
    Science Advisor
    Homework Helper

    Well, I worked it through also but I left the coefficients as [itex]a_n[/itex] and ended up with an expression:

    [tex]\sum_{n=1}^{\infty}[(m+2n)(m+2n-1)a_n+(m+2n)a_n+a_{n-1}-m^2a_n]x^{m+2n}[/tex]

    This reduces to:

    [tex]\sum_{n=1}^{\infty}[4n(m+n)a_n+a_{n-1}]x^{m+2n}[/tex]

    Now, substitute the expression for [itex] a_n[/itex] and [itex]a_{n-1}[/itex] and the results come out to zero.
     
  21. Apr 1, 2005 #20
    Yeah. Essentially the same thing I did, except without writing everything out each time! Of course I just copied and pasted anyways, so it doesn't matter too much~
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?



Similar Discussions: Bessel function and Bessel D.E.
Loading...