El comportamiento asintótico (n\gg 1) está dada por (el x contribución desaparece y la integral se convierte en solucionable) :
\int_0^\infty e^{\frac{x -x^3}{3n}}dx \sim (3n)^{1/3}\Gamma\left(\frac 43\right)
Para obtener el siguiente términos de la expansión, y darle un poco de confianza en este resultado, vamos a ampliar el e^{x/(3n)} factor en la serie y uso (ya que esta integral puede escribirse como una 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)
Observar que tendremos que dividir por (3n)^{2/3} a cada paso k\to k+1 (como se puede ver en GEdgar la respuesta) llegar con la expansión de la 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)