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: Indefinite integral (Hermite polynomials)

  1. Mar 30, 2010 #1
    1. The problem statement, all variables and given/known data
    I need to evaluate the following integral:

    [tex]\int_{-\infty}^{\infty}x^mx^ne^{-x^2}dx[/tex]

    I need the result to construct the first 5 Hermite polynomials.
    2. Relevant equations



    3. The attempt at a solution
    First I tried arbitrary values for "m" and "n". I was not able to evaluate the integral when (m+n) was even. I found that when (m+n) is odd, the integral is equal to 0. This makes sense since e^(-x^2) is an even function.

    Then I looked it up. Using wolfram integral calculator left me with something known as the incomplete gamma function (Not very sure what that is) in the solution. Maple gave me a very long answer that involved something labeled as Whittaker M.

    Finally I went to my professor to ask him how to evaluate it. He said that it is an integral that can be looked up and I should have a square root pi in my answer. He also said that the incomplete gamma function should not be there.

    I am not too sure of how to evaluate the integral myself. Maybe more importantly, I am not sure where to find the result of the integral if it can indeed be "looked up". Thanks for reading.
     
  2. jcsd
  3. Mar 30, 2010 #2

    vela

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper
    Education Advisor

    When the integrand is even, you can rewrite it with limits from 0 to infinity. Then use the substitution u=x^2, and it'll be in the form of a gamma function.
     
  4. Mar 30, 2010 #3
    Thanks vela!

    I am not sure how the gamma function works. I have just looked over articles pertaining to it online, but I am not confident enough in my understanding of it to be able to identify it in this case.

    Here is what I have following vela's advice:

    [tex]\int_{0}^{\infty}x^{2n}e^{-x^2}dx[/tex]

    Letting u=x^2 I get:

    [tex]\int_{0}^{\infty}\frac{u^{n-1}}{2}e^{-u}du[/tex]

    Integrating by parts yields:

    [tex]\int_{0}^{\infty}\frac{u^{n-1}}{2}e^{-u}du= -\frac{u^{n-1}e^{-u}}{2}+\int_{0}^{\infty}\frac{(n-1)u^{n-2}}{2}e^{-u}du[/tex]

    Integrating by parts yields more repeating integrals (as far as I can tell from here). Each time I integrate by parts I can pull out a factor. For the first one it was:

    [tex]\frac{n-1}{2}[/tex]

    The next factor (I guess I'm talking about the coefficient of the "uv" of the integration by parts...if you use uv-int(vdu)) I get is:

    [tex]\frac{(n-1)(n-2)}{2}[/tex]

    This keeps repeating. I see that according to wikipedia, the gamma function can be defined as:

    [tex]\Gamma(n)=(n-1)![/tex]

    As of this point, I do not see how this will get me where I want to go. I am looking for the result of this integral in order to establish rules to determine the orthogonality relationships between (x^m,x^n).

    For instance when finding the Legendre polynomials, I found that when (m+n) was odd, (x^m,x^n) = 0. When (m+n) was even, (x^m,x^n) = (2/(m+n+1)). I determined this by performing the integral in the first post. The difference being the weight function.

    Thanks again for taking the time to read this.
     
  5. Mar 30, 2010 #4

    vela

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper
    Education Advisor

    You made a mistake when you did the substitution. Because du=2x dx, you should get

    [tex]\frac{1}{2}\int_0^\infty u^{n-1/2}e^{-u} du = \frac{1}{2}\int_0^\infty u^{(n+1/2)-1} e^{-u} du = \frac{1}{2}\Gamma(n+\frac{1}{2})[/tex]

    You can look up what that evaluates to, or deduce it using the facts that [itex]\Gamma(n+1)=n\Gamma(n)[/itex] and [itex]\Gamma(1/2) = \sqrt{\pi}[/itex].

    I didn't realize you weren't familiar with the gamma function. I would have suggested a different way to go about this integral. For instance, you could use the trick where you say

    [tex]\int x^{2n} e^{-x^2} dx = \left[\int (-1)^n \frac{d^n}{d\alpha^n} e^{-\alpha x^2} dx\right]_{\alpha=0}[/tex]

    and then interchange the order of integration and differentiation. But then again, that's kind of tedious, so maybe learning about the gamma function isn't so bad. :)
     
  6. Mar 30, 2010 #5
    Thanks for catching my mistake! So after scouring the internet (I wasn't able to figure out the quantity by searching or using the 2 bits of information) and putting together all I found, I think:

    [tex]\Gamma(n+\frac{1}{2})=\frac{2n!}{n!2^{2n}}\sqrt\pi[/tex]

    I realize the above quantity needs to be multiplied by 1/2 to get the quantity I am looking for. Is this correct? If so, can it be simplified. I tried but so far have not been able to re-write it. In fact, I was only able to get the above equation by putting 3 or so different equations together. I don't fully understand why the (pi)^.5 can "come out" of the gamma function.

    Also, I didn't understand your other method of finding this relation. How could you do this without knowledge of the gamma function?
     
  7. Mar 30, 2010 #6

    vela

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper
    Education Advisor

    I think it's right. Wikipedia's page on the gamma function actually has the result you were looking for, although it's written slightly differently.

    When n=0, the integral reduces to a gaussian. That's why you get the square root of pi in the answer.
    It's just a common trick for doing these sorts of integrals. There's a mistake in what I wrote. You want to set alpha=1 at the end, not to zero.
     
  8. Mar 31, 2010 #7
    Now it seems that I either have the wrong solution, or I do not know what I am doing. When I use the relation :


    [tex]\int_0^\infty{x^m}x^n{e^{-x^2}dx=\frac{1}{2}((m+n)-1/2)![/tex] if (m+n) = even

    I'm assuming the gamma function below is equal to the right hand side above

    [tex]\frac{1}{2}\Gamma(n+\frac{1}{2})=\frac{1}{2}\frac{2n!}{n!2^{2n}}\sqrt\pi[/tex]

    I thought this value would allow me to generate the first few Hermite polynomials but I am not sure now how to do that.

    Using the Schmidt orthogonalization process, as outlined in my tutorial, I take:

    [tex]H_0(x)=\frac{x^0}{\sqrt{(x^0,x^0)}}[/tex]

    Using the above equations for n=0 gives me:

    [tex](x^0,x^0)=\frac{\sqrt{\pi}}{2}[/tex]

    Shouldn't the first polynomial = 1?

    To continue, the next step is to take

    [tex]H_1(x)=x^1+\alpha{H_0}[/tex]

    and solve for alpha to find H_1. Am I doing something incorrect here?
     
    Last edited: Mar 31, 2010
  9. Mar 31, 2010 #8

    vela

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper
    Education Advisor

    The difference is normalization. The Gram-Schmidt method needs the functions to be normalized such that (Hn,Hn)=1, but the Hermite polynomials are typically normalized to make them more convenient to use. The http://en.wikipedia.org/wiki/Hermite_polynomials" [Broken] discusses this a bit. What you will find using your method are polynomials that are proportional to the Hermite polynomials.
     
    Last edited by a moderator: May 4, 2017
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook