I am trying to solve a differential equation of the form:
\begin{equation}x^2y''+2xy'+x^2e^{ay}=0\end{equation}
This arises from calculating the electric potential of ions following the Boltzmann distribution attracted to a spherically symmetric electrode for those who are curious.
Anyway, seeing as this is very similar to #33 from here, I made the substitutions $z=x^2e^{ay}$ and $w=xy'$. I then found
\begin{equation}dz=2xe^{ay}dx+ax^2e^{ay}y'dx=(2+w)z\frac{dx}{x}\end{equation}
\begin{equation}dw=y'dx+xy''dx\end{equation}
I therefore got
\begin{equation}\frac{dw}{dz}=w'=\frac{x^2y''+xy'}{(2+aw)z}\end{equation}
Substituting this and the definition of $z$ into the original ODE gives
\begin{equation}z(2+aw)w'+w+z=0\end{equation}
I thought this would be easier to solve, as it is a 1st order ODE. However, it is non-linear and not separable so far as I can see. In the original from the link there is no 2 on the $y'$ term, which makes the transformed ODE separable.
Can anyone confirm that what I did with the substitution was correct? Also does anyone have any other suggestions for solving either the original or transformed ODE analytically?