Check that both the change of variables $\eta=\ln \xi, \eta=\xi^{2 / \sigma} \varphi(\eta), \psi=d \varphi / d \eta$ or the replacement $\varphi=-\xi f^{-1} f^\sigma d f / d \xi, \psi=-\xi f^{1-\sigma} /(d f / d \xi)$ reduces $(13)$ to a first order equation of a form $d \psi / d \varphi=\psi F(\psi, \varphi) /(\varphi \Phi(\psi, \varphi))$ (in the second case the equation is easier to investigate, since the functions $F, \Phi$ contain simple nonlinearity of a form $\psi \varphi$ ).