This question already has been solved here: How to evaluate $\int_{0}^{\infty}\exp(-x^2-1/x^2)dx$?
And there is also an answer for How should I calculate $\int_0^\infty e^{-\frac{1}{2}(x^2+a^2/x^2)}\,dx$ But I think the answer is incorrect, if this problem solved, and I can try to replace x by u^2 in my question, even this would still be a little different because there would have another more 2u in the equation, but would it still be some help to solve my problem?