Exact late slope of an exponentially damped oscillator (${\lim }_{t\to \infty }{y}_t=e^{-1/2}M(-1/(2{\epsilon }^2),1,1/2)$)

For $y_{tt}+e^{-2\epsilon t}(\epsilon y_t+y)=0$, $y(0)=0$, $y_t(0)=1$, set $r=e^{-2\epsilon t}/2$. The exact equation becomes $r{y}_{rr}+(1-r)y_r+y/(2{\epsilon }^2)=0$. Let $m(t)=M(-1/(2{\epsilon }^2),1,r)$, using the Kummer function. It is regular at $r=0$ and tends to one. The Wronskian $W=m{y}_t-m_{t}y$ obeys $W_t=-\epsilon e^{-2\epsilon t}W$, with $W(0)=M(-1/(2{\epsilon }^2),1,1/2)$. Integration gives the displayed exact late slope. The other local solution grows at most logarithmically in $r$, so $m_{t}y\to 0$ and $W\to y_t$. This provides a check on a WKB approximation even near phases where its leading predicted slope vanishes.