WKB approximation for a slowly varying oscillator (source code)

= WKB approximation for a slowly varying oscillator

For $\epsilon^2y_{XX}+k^2(X)y=0$ with nonzero slowly varying $k$,
$$
y\sim k(X)^{-1/2}\left[
C_+e^{i\epsilon^{-1}\int^Xk(s)ds}
+C_-e^{-i\epsilon^{-1}\int^Xk(s)ds}
\right].
$$
It fails near turning points and where $\epsilon|k'|/|k|^2$ is not small.