Stability analysis of orbital modes for a generalized Lane-Emden equation

We present a stability analysis of the standard nonautonomous systems type for a recently introduced generalized Lane-Emden equation which is shown to explain the presence of some of the structures observed in the atomic spatial distributions of magnetically-trapped ultracold atomic clouds. A Lyapunov function is defined which helps us to prove that stable spatial structures in the atomic clouds exist only for the adiabatic index $\gamma =1+1\mathrm{/}n$ with even $n$. In the case when $n$ is odd we provide an instability result indicating the divergence of the density function for the atoms. Several numerical solutions, which according to our stability analysis are stable, are also presented.