4
$\begingroup$

How can you compute the following integral assuming $n>0$? $$\int_{x=0}^{\infty}e^{\frac{x -x^3}{3n}}dx $$

Mathematica etc. fail to produce anything useful.

EDIT: I would be happy with an asymptotic result in $n$ if it is too hard to compute exactly. I don't know if it helps but $$\int_{x=0}^{\infty}e^{\frac{-x^3}{3n}}dx \approx 1.29 n^{1/3}.$$

  • 0
    I think numerical approximation may be you're best bet. You're going to end up with some nasty-ish infinite series to evaluate IMHO.2012-07-11

2 Answers 2

5

The asymptotic behavior ($n\gg 1$) is given by (the $x$ contribution disappears and the integral becomes solvable) : $$\int_0^\infty e^{\frac{x -x^3}{3n}}dx \sim (3n)^{1/3}\Gamma\left(\frac 43\right)$$

To get the next terms of the expansion, and give you some confidence in this result, let's expand the $e^{x/(3n)}$ factor in series and use (since this integral may be rewritten as a gamma integral) : $$\int_0^\infty \left(\frac x{3n}\right)^k e^{\frac{-x^3}{3n}}dx=\frac {(3n)^{(1-2k)/3}}{k+1}\Gamma\left(\frac{k+4}3\right)$$

Observe that we will have to divide by $(3n)^{2/3}$ at each step $k\to k+1$ (as you may see in GEdgar's answer) getting with the expansion $e^r=1+r+\frac {r^2}2+\cdots$

$$\int_0^\infty e^{\frac{x -x^3}{3n}}dx \sim (3n)^{1/3}\left(\Gamma\left(\frac 43\right)+\frac {\Gamma\left(\frac 53\right)}{2(3n)^{2/3}}+\frac 1{2\cdot 3 (3n)^{4/3}}+\rm{O}\left(n^{-2}\right)\right)$$

  • 0
    How can you rigorously show that the $e^{x}$ term has no effect on the asymptotics of this integral?2012-07-11
  • 0
    @Arlo: you are welcome!2012-07-13
3

I guess Maple has one of the useless answers: $$\int_{0}^{\infty} \operatorname{e} ^{\frac{-x (-1 + x^{2})}{3 n}} d x = \frac{i}{9} \pi \mathrm{BesselY} \biggl(\frac{1}{3},\frac{2 i \sqrt{3}}{27 n}\biggr) \sqrt{3} - \\ \quad{}\quad{}\frac{i}{9} \pi \mathrm{AngerJ} \biggl(\frac{1}{3},\frac{2 i \sqrt{3}}{27 n}\biggr) + \frac{2 i}{3} \mathrm{BesselK} \biggl(\frac{1}{3},\frac{-2\sqrt{3}}{27 n}\biggr) + \\ \quad{}\quad{}\frac{i}{9} \pi \mathrm{BesselJ} \biggl(\frac{1}{3},\frac{2 i \sqrt{3}}{27 n}\biggr) + \frac{i}{9} \sqrt{3} \pi \mathrm{WeberE} \biggl(\frac{1}{3},\frac{2 i \sqrt{3}}{27 n}\biggr) $$ Maple reports asymptotics for this as $$ \frac{2 \sqrt[3]{n} \pi 3^{\frac{5}{6}}}{9 \Gamma \Bigl(\frac{2}{3}\Bigr)} $$ which differs by factor $3^{5/6}$ from Raymond's answer. Perhaps this is more accurate??

More terms $$ \frac{2 \sqrt[3]{n} \pi 3^{\frac{5}{6}}}{9 \Gamma \Bigl(\frac{2}{3}\Bigr)} + \frac{\Gamma \Bigl(\frac{2}{3}\Bigr) 3^{\frac{2}{3}}}{9 \sqrt[3]{n}} - \frac{1}{36 n} + O(n^{-5/3}) $$

  • 0
    Yes. I can't even tell by looking it if that is a real number or not. In fact, if you do int(exp((x-x^3)/100),x=0..infinity); in Maple you get a number that is not real! I don't think it can cope with this.2012-07-11
  • 0
    you are right about accuracy GEdgar (I forgot the 3 in at the bottom in one of my computations, I'll update that!)2012-07-11
  • 0
    @Arlo: If you evaluate int(exp((x-x^3)/100),x=0..infinity); to 20 places, you get $4.24 + 1.8\times 10^{-15}i$, so I think it is actually real.2012-07-11
  • 0
    @GEdgar: How do you get Maple to report the asymptotics?2012-07-11
  • 0
    @Gedgar: just to be clear $3^{1/3}\Gamma(4/3)=\dfrac{2\pi 3^{5/6}}{9 \Gamma(2/3)}$.2012-07-11
  • 0
    In Maple, asympt(U,n,2); for asymptotics of U up to "order 2"2012-07-11
  • 0
    +1: fine and better than my answer now (except the $3^{5/6}$ that I changed anyway...)2012-07-11
  • 1
    @Raymond: yours now agrees with mine2012-07-11