How to compute $\int_{-1}^{1} e^{-1/(1-x^2)}dx$? As in the title, I would like to compute the integral:
\begin{equation}
\int_{-1}^{1}e^{-1/(1-x^2)}dx
\end{equation}
My hunch tells me that I should try to transform it to the correspoding $\int_{-1}^{1} e^{-y^2}dy$ but I am not sure. 
Thank you for your time!

$$\begin{eqnarray*}\color{red}{I}=\int_{-1}^{1}\exp\left(\frac{1}{x^2-1}\right)\,dx &=& 2\int_{0}^{1}\exp\left(\frac{1}{x^2-1}\right)\,dx=\int_{0}^{1}\exp\left(\frac{1}{z-1}\right)\frac{dz}{\sqrt{z}}\\&=&\int_{0}^{1}\exp\left(-\frac{1}{z}\right)\frac{dz}{\sqrt{1-z}}=\int_{1}^{+\infty}\frac{dt}{t e^t\sqrt{t^2-t}}\\&=&\int_{0}^{+\infty}\frac{e^{-(u+1)}du}{(u+1)\sqrt{u(u+1)}}=\frac{2}{e}\int_{0}^{+\infty}\frac{e^{-\eta^2}\,d\eta}{(1+\eta^2)^{3/2}}\\&=&\color{red}{\frac{\sqrt{\pi}}{e}\,U\left(\frac{1}{2},0,1\right)}\end{eqnarray*}$$
where $U(a,b,z)$ is the Tricomi's confluent hypergeometric function.
If we take the last integral and switch to Fourier transforms, that can be written also as:
$$\frac{2}{e\sqrt{\pi}}\int_{0}^{+\infty} e^{-s^2/4} s\, K_1(s)\,ds $$
where $K_1$ is a modified Bessel function of the second kind. Tight numerical approximations follows from the fact that the last integrand function is smooth and essentially supported on $[0,4]$, since the integral over $[4,+\infty)$ is extremely small. We also have:
$$ \int_{-1}^{1}\exp\left(\frac{1}{x^2-1}\right)\,dx = \frac{2}{e}\int_{0}^{\pi/2}e^{-\tan^2 t}\cos t\,dt.$$
Trivial inequalities are:
$$\color{red}{0.412\ldots}=\sqrt{\frac{2\pi}{5e^2}}=\frac{2}{e}\int_{0}^{+\infty}e^{-5u^2/2}\,du \leq \color{red}{I}\leq \frac{2}{e}\int_{0}^{+\infty}\frac{du}{(1+u^2)^{5/2}}=\frac{4}{3e} =\color{red}{0.490\ldots}$$
