In this paper, we introduce a new method for calculation of energy levels in detail and give our results for several iso-spectrum-level series as examples: [He] where I(Z), T lim (Z), and T (Z, n) denote ionization potential, series limit, and energy level of a given member, respectively. The expression of non-relativistic part of I(Z) is derived from weakest bound election potential model theory and relativistic effects of I(Z) are included by using a six-order polynomial in Z. Our results are compared with the experimental data and with those obtained by other theoretical method.