By using the Pekeris approximation, the Schrödinger equation is solved for the nuclear deformed Woods-Saxon potential within the framework of the asymptotic iteration method (AIM). The energy levels are worked out and the corresponding normalized eigenfunctions are obtained in terms of hypergeometric function. The deformed Woods-Saxon (dWS) potential is a short range potential and widely used in nuclear, particle, atomic, condensed matter and chemical physics [1][2][3][4][5][6][7]. This potential is reasonable for nuclear shell models and used to represent the distribution of nuclear densities.The dWS and spin-orbit interaction are important and applicable to deformed nuclei [8] and to strongly deformed nuclides [9]. The dWS potential parameterization at large deformations for plutonium 237,239,241 P u odd isotopes was analyzed [10]. The structure of single-particle states in