Past exam of the perturbation methods and applied complex variables course of the University of Oxford / 2010 / A / 1 / i / Solution

Put $\delta =\sqrt{ϵ}$ and write the outer asymptotic expansion as $y_0+\delta y_{1/2}+{\delta }^2{y}_1+\cdots$. Equating powers in the differential equation givesThe first two equations give $y_0=A{e}^{x^2/2}$ and $y_{1/2}=B{e}^{x^2/2}$. Since $y_0^{(3)}=A(x^3+3x)e^{x^2/2}$, an integrating factor yieldsHere $A,B,C$ are matching constants; boundary conditions have not yet been imposed.