Problem: Find a uniformly valid leading order approximation to the solution of $\epsilon y'' + \tanh^2(x)y + \tanh(x)y'=\tanh(x)\text{sech}(x)$ with boundary conditions $y(-2)=0, y(2)=0$ in the limit $\epsilon \to 0^+$. Use only the variables and constants given in the problem; do not define additional constants. Place your final solution in a $\boxed{}$ LaTeX environment.
Solution: $\boxed{\frac{x-2\text{erf}(\frac{x}{\sqrt{2\epsilon}})}{\cosh(x)}}$