This is homogeneous class D \(y^{\prime }=\frac {y}{x}+g\left ( x\right ) F\left ( \frac {y}{x}\right ) \). Hence from lookup table\begin {align*} \xi & =x^{2}\\ \eta & =xy \end {align*}
From above we found the solution to be\[ S=\int \frac {1}{F\left ( R\right ) }dR+c \] In this case \(F\left ( R\right ) =e^{-R}\). Hence\begin {align*} S & =\int e^{R}dR+c\\ S & =e^{R}+c \end {align*}
Now we just need to find canonical coordinates \(\left ( R,S\right ) \) since \(\xi ,\eta \) are known. From above\begin {align*} R & =\frac {y}{x}\\ S & =-\frac {1}{x} \end {align*}
Hence the solution becomes\begin {align*} -\frac {1}{x} & =e^{\frac {y}{x}}+c\\ e^{\frac {y}{x}} & =c_{2}-\frac {1}{x}\\ \frac {y}{x} & =\ln \left ( c_{2}-\frac {1}{x}\right ) \\ y & =x\ln \left ( c_{2}-\frac {1}{x}\right ) \end {align*}
The nice thing about this method is that once we solve for one pattern of an ode, then the same solution in canonical coordinates is used, the only change need is to plug-in in the RHS of the original ode in the solution and integrate.