The energy levels of neutral atoms supported by Yukawa potential, $V(r)=-Z
exp(-\alpha r)/r$, are studied, using both dimensional and dimensionless
quantities, via a new analytical methodical proposal (devised to solve for
nonexactly solvable Schrodinger equation). Using dimensionless quantities, by
scaling the radial Hamiltonian through $y=Zr$ and $\alpha^{'}=\alpha/Z$, we
report that the scaled screening parameter $\alpha^{'}$ is restricted to have
values ranging from zero to less than 0.4. On the other hand, working with the
scaled Hamiltonian enhances the accuracy and extremely speeds up the
convergence of the energy eigenvalues. The energy levels of several new
eligible scaled screening parameter $\alpha^{'}$ values are also reported.