AbstractSolution of the radial Schrodinger equation for the Woods–Saxon potential together with spin–orbit interaction, coulomb and centrifugal terms by using usual Nikiforov–Uvarov (NU) method is not possible. Here, we have presented a new NU procedure with which we are able to solve this Schrodinger equation and any other one-dimensional ones with any shape of the potential profile. For this purpose, we have combined the NU method with numerical fitting schema. The energy eigenvalues and corresponding eigenfunctions for various values of n, l, and j quantum numbers have been obtained. Good agreement with experimental values is also achieved. We have calculated the 1/2+ state energy with more accuracy (our absolute error = 0.023 MeV and Hagen et al. absolute error = 0.0918 MeV), while Hagen et al. have calculated the 5/2+ state energy with higher accuracy (our absolute error = 0.71 MeV and Hagen et al. absolute error = 0.0337 MeV). Our wave functions are in agreement with Kim et al.’s work, too. پرونده مقاله