Problem: Find a uniformly valid leading order approximation to the solution of $\epsilon y'' - x^2y' - (3+x^3) = 0$ with boundary conditions $y(1) = 1; y(2) = 1$ 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_{uniform}(x) = \frac{3}{x} -\frac{x^2}{2} -\frac{3}{2} + 3e^{\frac{-4(2-x)}{\epsilon}}}$