A novel method to find bound states in general tight‐binding Hamiltonians with semi‐infinite leads is reported. The method is based on the recursive S‐matrix method, which allows to compute iteratively the S‐matrix of a general system in terms of the S‐matrices of its subsystems. The condition that the S‐matrices of the subsystems must accomplish to have a bound state at energy E is established. Energies that accomplish this relation, can be determined with high accuracy and efficiency by using the Taylor series of the S‐matrices. The method allows to find bound states energies and wavefunctions in (BIC) and out (BOC) of the continuum, including degenerate ones. Bound states in nanoribbons with wider sections are computed for square and honeycomb lattices. Using this method, the bound states in a graphene nanoribbon with two quantum‐dot‐like structures which are reported to have BICs by using another technique are verified. However, this new analysis reveals that such BICs are double, one with even and the other with odd wavefunction, with slightly separated energies. In this way, the new method can be used to efficiently find new BICs and to improve precision in previously reported ones.