Problem: Find a uniformly valid leading order approximation to the solution of $\epsilon y'' + \tanh(x)y' + tanh^2(x)y=tanh^2(x)$ with boundary conditions $y(-2)=1, y(2)=2$ 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{y(x)=1+\frac{\cosh(2)}{2\cosh(x)}(1+\text{erf}(\frac{x}{\sqrt{2\epsilon}}))}$