NOTA:
Quería dar un agradecimiento especial a @robjon por sus acertados comentarios.
En primer lugar, observamos que $\lim_{\epsilon\to 0}e^{-\tan z/\epsilon}=0$ a menos que $z=\ell \pi$ , $\ell$ un número entero. Por lo tanto, toda la "acción" de la integración tendrá lugar en intervalos alrededor de $\ell \pi$ . Así pues, veamos primero lo que ocurre para $0<x<\pi/2$ .
En el espíritu de Método de Laplace tenemos para $0<z<\pi/2$ , $\tan^z =z^2+O(z^4)$ y por lo tanto para $0<x<\pi/2$
$$\begin{align} \epsilon^{-1/2}\int_0^xze^{-\tan^2z/\epsilon}dz&\sim\epsilon^{-1/2}\int_0^xze^{-z^2/\epsilon}dz\\\\ &=\epsilon^{-1/2}\left.\left(-\epsilon^{-z^2/\epsilon}\right)\right|_{z=0}^{z=x}\\\\ &=\epsilon^{1/2}\left(1-e^{-x^2/\epsilon}\right) \end{align}$$
que claramente va a cero como $\epsilon\to 0$ .
A continuación, observamos que la integración alrededor de las singularidades de la función tangente no plantea ningún problema. Así, para una $(L-1)\pi<x<L\pi$ y $\delta >0$ podemos escribir
$$\begin{align} \epsilon^{-1/2}\int_0^x ze^{-\tan^2z/\epsilon}dz&=\epsilon^{-1/2}\sum_{\ell=0}^{L-2}\left(\int_{\ell \pi+\delta}^{(\ell+1)\pi-\delta}ze^{-\tan^2z/\epsilon}dz+\int_{(\ell+1)\pi-\delta}^{(\ell+1)\pi+\delta}ze^{-\tan^2z/\epsilon}dz\right)\\\\ &+\epsilon^{-1/2}\int_{(L-1)\pi+\delta}^{x}ze^{-\tan^2z/\epsilon}dz \tag 1\\\\ \end{align}$$
Observamos que en $(1)$ las únicas integrales que contribuirán en el límite como $\epsilon \to 0$ son aquellos en torno a múltiplos enteros de $\pi$ . Así, tenemos para $(L-1)\pi<x<L\pi$ y $\delta>0$
$$\begin{align} \lim_{\epsilon \to 0}\epsilon^{-1/2}\int_0^x ze^{-\tan^2z/\epsilon}dz&=\lim_{\epsilon \to 0} \epsilon^{-1/2}\sum_{\ell=0}^{L-2}\left(\int_{(\ell+1)\pi-\delta}^{(\ell+1)\pi+\delta}ze^{-\tan^2z/\epsilon}dz\right) \tag 2\\\\ \end{align}$$
Procedemos a evaluar las integrales en $(2)$ . Para ello tenemos
$$\begin{align} \epsilon^{-1/2}\int_{(\ell+1)\pi-\delta}^{(\ell+1)\pi+\delta}ze^{-\tan^2z/\epsilon}dz &=\epsilon^{-1/2}\left(\int_{-\delta}^{\delta}ze^{-\tan^2z/\epsilon}dz+(\ell +1)\pi\int_{-\delta}^{\delta}e^{-\tan^2z/\epsilon}dz\right)\\\\ &=(\ell +1)\pi\epsilon^{-1/2}\int_{-\delta}^{\delta}e^{-\tan^2z/\epsilon}dz\\\\ &\sim (\ell +1)\pi\epsilon^{-1/2}\int_{-\delta}^{\delta}e^{-z^2/\epsilon}dz\\\\ &= (\ell +1)\pi\int_{-\delta/\epsilon^{1/2}}^{\delta/\epsilon^{1/2}}e^{-z^2}dz\\\\ &\to (\ell +1)\pi^{3/2} \end{align}$$
Suma de $\ell$ encontramos para $(L-1)\pi<x<L\pi$
$$\lim_{\epsilon \to 0}\epsilon^{-1/2}\int_0^xze^{-\tan^2z/\epsilon}dz=\frac{L(L-1)\pi^{3/2}}{2}$$
Una última nota se refiere al caso en que $x=L\pi$ . Para este caso, vemos que hay que añadir una integral más, a saber
$$\begin{align} \lim_{\epsilon\to 0}\epsilon^{-1/2}\int_{L\pi-\delta}^{L\pi}ze^{-\tan^z/\epsilon}&=L\pi\int_{-\infty}^0e^{-z^2}dz\\\\ &=\frac12 L\pi^{3/2} \end{align}$$
Así, para $x=L\pi$ tenemos
$$\lim_{\epsilon \to 0}\epsilon^{-1/2}\int_0^xze^{-\tan^2z/\epsilon}dz=\frac{L^2\pi^{3/2}}{2}$$
Si lo juntamos todo, tenemos
$$\lim_{\epsilon \to 0}\epsilon^{-1/2}\int_0^xze^{-\tan^2z/\epsilon}dz= \begin{cases} \frac{L(L-1)\pi^{3/2}}{2},&(L-1)\pi<x<L\pi\\\\ \frac{L^2\pi^{3/2}}{2},&x=L\pi \end{cases} $$