The Yukawa potential (or screened Coulomb) describes interactions screened by neighbors in many fields of physics: solid state, atomic, nuclear and etc. Especially, in high energy physics, the Yukawa potential is used for the phenomenological quark-quark interaction which, however, needs to be solved numerically since it has no analytic solution. In this work, we try to find the efficient way of the numerical method based on basic principles of Qunatum mechanics; particle states can be expanded with any complete set of eigenstates. We test the validity of this method by solving the Schroedinger equation for the Coulomb potential which has analytic solutions, and we abstract serveral factors affecting the accuracy of solutions: the resemblance between the reference potential and the potential in question, the number of eigenbases and the range of principal quantum numbers of eigenbases.

To obtain eigenenergies and eigenstates for the Yukawa potential, we choose the Coulomb potential as the reference potential and its eigenstates from ground to 100th states. We compare our results to other published ones: Direct root finding method (Direct), Ritz method with Dirichelet boundary condition (Ritz), Super Symmetric quantum mechanics perturbation method (SUSY), Asymptotic Iteration Method (AIM) and Nikiforov-Uvarov method (NU). The relative errors for eigenenergies are within $2.2\%$ and those for expectation values of the position are within $0.075\%$. Also we check that the stronger screening effect makes eigenenergies deeper and wavefunctions wider. This effect is more significant for higher states.