In this notebook we provide a code for the calculation of eigenvalues and eigenstates for the bound states of the Yukawa Potential following the algorithm described in the following papers i) Phys.Lett.B 816 (2021) 136218 (arXiv: 2012.12969), ii) arXiv:2102.07160, submitted for publication to Progress of Theoretical and Experimental Physics (PTEP); you can freely use it but we kindly ask you to cite these papers as the original source.