We investigate the bound states of the Yukawa potential V (r) = −λ exp(−αr)/r, using different algorithms: solving the Schrödinger equation numerically and our Monte Carlo Hamiltonian approach. There is a critical α = α C , above which no bound state exists. We study the relation between α C and λ for various angular momentum quantum number l, and find in atomic units, α C (l) = λ[A 1 exp(−l/B 1 ) + A 2 exp(−l/B 2 )], with A 1 = 1.020(18), B 1 = 0.443(14), A 2 = 0.170(17), and B 2 = 2.490(180).