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!

Tricky integral

  1. Sep 27, 2005 #1
    I just can't solve this:
    [tex]I(\gamma)=\int_{-\infty}^{\infty} e^{-\sqrt{x^2+\gamma^2}}dx[/tex]
    I have tried a trigonometric substitution
    [tex]x=\gamma \tan(\theta)[/tex]
    but I am not happy with the result :grumpy:.
    Anyone got a hint?
     
  2. jcsd
  3. Sep 27, 2005 #2

    hotvette

    User Avatar
    Homework Helper

    Try:

    [tex]y=\sqrt{x^2+\gamma^2}[/tex]

    hotvette
     
    Last edited: Sep 27, 2005
  4. Sep 27, 2005 #3
    Then
    [tex]dy=\frac{x}{\sqrt{x^2+\gamma^2}}dx[/tex]
    and I don't quite see it happening...
    Can you be more precise?
     
  5. Sep 27, 2005 #4

    hotvette

    User Avatar
    Homework Helper

    Just go a bit further. You can get dx in terms of dy, y, and [tex]\gamma[/tex] by using the substitution and some algebra. By the way, the suggestion was to simplify the expontential. Hopefully, the resulting integrand will be more reasonable to deal with, such as integration by parts. I haven't solved it, but the substitution looks like it takes you in a promising direction.

    hotvette
     
  6. Sep 27, 2005 #5
    Looks dirty. :D
    Numerical integration time! :rofl:
     
  7. Sep 27, 2005 #6

    Tom Mattson

    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    Is this really a high school calculus problem? :eek:

    How about taking the natural log of both sides, integrating, and then exponentiating the result? Is it permissible to interchange the operations of ln and integration here?
     
  8. Sep 27, 2005 #7
    Confusing the function is govern by gamma, but you are doing an intergration of x. Do you mean that the function is govern both by x and gamma?
     
  9. Sep 27, 2005 #8

    Tom Mattson

    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    No, the x is integrated out. To see this, consider the case [itex]\gamma=1[/itex]:

    [tex]I(1)=\int_{-\infty}^{\infty} e^{-\sqrt{x^2+1}}dx[/tex]

    When that integral is evaluated, you will get a number, not a function of x.
     
  10. Sep 27, 2005 #9

    Tom Mattson

    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    Apparently, it isn't. I tried it for the indefinite integral, and the antiderivative I obtained is wrong.

    Actually it should have been obvious to me that this won't work, otherwise someone would have used it to find an antiderivative of [itex]exp(-x^2)[/itex]. Shoulda knowed better.
     
  11. Sep 27, 2005 #10
    Lots of answers, thanks guys! :smile: but still no result :cry:

    hotvette: I am not sure about this but I'll try. The pain is you have to split the integral because of the term [tex]\sqrt{y^2-\gamma^2}[/tex]

    Pseudo Statistic: I do not want a numerical result this time... sorry.

    Tom Mattson: No this is not a high school problem, any recommendation on where to post it? Good try with the log but integrals unfortunately do not have this property, I can not even think of any non trivial integral (that is the integral of 1 or that of 0) that would have this property.
     
  12. Sep 27, 2005 #11

    Tom Mattson

    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    If it's a homework assignment, then it is most likely at the College Level. (In the USA, College=University).

    I had another attempt that didn't work, but maybe you can play with it to get something. The trick is called "parameterizing the integral".

    I started by generalizing your function to include a second variable.

    [tex]F(\gamma,a)=\int_{-\infty}^{\infty} e^{-\sqrt{ax^2+\gamma^2}}dx[/tex]

    Then I differentiated with respect to [itex]a[/itex]:

    [tex]\frac{\partial F(\gamma,a)}{\partial a}=\int_{-\infty}^{\infty}\frac{\partial}{\partial a} e^{-\sqrt{ax^2+\gamma^2}}dx[/tex]

    [tex]\frac{\partial F(\gamma,a)}{\partial a}=\int_{-\infty}^{\infty}\frac{-a}{\sqrt{ax^2+\gamma^2}} e^{-\sqrt{ax^2+\gamma^2}}dx[/tex]

    Then I tried to do a u-substitution:

    [tex]u=\sqrt{ax^2+\gamma^2}[/tex]
    [tex]du=\frac{ax}{\sqrt{ax^2+\gamma^2}}[/tex]

    The idea is to integrate the new integral with respect to [itex]x[/itex], then integrate that with respect to [itex]a[/itex], and then recognize that [itex]I(\gamma)=F(\gamma,1)[/itex].

    Of course, I couldn't get it to work because of that [itex]x[/itex] in my [itex]du[/itex]. You can solve for [itex]x[/itex] in terms of [itex]u[/itex], but it gets ugly again. But maybe you can find a parameterization that works.

    And then again, there's always numerical integration. :frown:
     
    Last edited: Sep 27, 2005
  13. Sep 27, 2005 #12
    hotvette: Now I see what you meant :smile: I am taking that direction.

    Tom Mattson: Funny you would think of that, I was also looking into finding a differential equation for which "my" integral would be a solution.

    I'll let you know guys, I've got to break this one!
     
  14. Sep 27, 2005 #13

    hotvette

    User Avatar
    Homework Helper

    Using the substitution [tex]y = \sqrt{u^2+\gamma^2}[/tex] I as able to come up with:

    [tex]\int\frac{ye^{-y}}{\sqrt{y^2-\gamma^2}}dy[/tex]

    which seems a bit more reasonable. I then tried integration by parts, which resulted in:

    [tex]e^{-y}\sqrt{y^2 - \gamma^2} + \int e^{-y}\sqrt{y^2-\gamma^2}dy[/tex]

    but got stuck there.

    hotvette
     
  15. Sep 28, 2005 #14
    That's where I stand too... noting that given the bounds the first term is zero.
    I'm trying to crack it with Maple now (o:)shame on me!) but nothing good so far.
     
  16. Sep 28, 2005 #15
    So far, I've figured out that since the integrand is an even function, instead of integrating from negative infinity to positive infinity you can integrate from 0 to positive infinity and multiply the result by two. Also, when gamma = 0, it's pretty easy to show that I(0) = 2*1 = 2. I haven't figured out the rest, though.
     
  17. Sep 28, 2005 #16
    Thanks JoAuSc, it's a very wise thing to search for simplifications and particular cases. The "rest" is the tough part, though.

    Anyway, I have the answer! :smile:
    I just couldn't find it in text books but this integral is well known as the "modified Bessel function of second kind". As expected it does not have a primitive :cry: but that's all right because it is the answer to my original question.

    Thanks all for helping out and sorry I had you sweat on such a tough problem.
     
  18. Sep 28, 2005 #17

    hotvette

    User Avatar
    Homework Helper

    Interesting. I'm curious as to how you found out. Can you post a web link?

    hotvette
     
  19. Sep 29, 2005 #18
    A friend a mine is a researcher in Maths and eventhough he is not a specialist in integration, a friend of his just recognised the function right away! Like: "Hey dude, this looks like [tex]K_1(\gamma)[/tex], one of the modified Bessel function of second kind :cool: "... nevermind.

    Anyway you can find some information here:
    http://en.wikipedia.org/wiki/Bessel_function
    and some more there:
    http://mathworld.wolfram.com/ModifiedBesselFunctionoftheSecondKind.html

    As this function is not usually defined by its integral (but by a more generic differential equation), you might not recognize it at first sight. The one we have actually corresponds to equation 7 in the second link, with n=1 and knowing that
    [tex]\frac{1}{2}!=\frac{\sqrt{pi}}{2}[/tex]

    The exact relation is:
    [tex]\int_0^{\infty}e^{-\sqrt{x^2+z^2}}dx=2zK_1(z)[/tex]
     
  20. Sep 29, 2005 #19

    hotvette

    User Avatar
    Homework Helper

    thanks. Glad you solved it.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?



Similar Discussions: Tricky integral
  1. Tricky Integrals (Replies: 6)

  2. A tricky integral (Replies: 14)

  3. Tricky Mechanics (Replies: 2)

  4. Tricky voltmeter (Replies: 4)

Loading...