A theoretical formalism for calculating the bulk and surface spin modes in Heisenberg semi-infinite lattices is presented on a ferromagnetic cubic network of spins, coupled via nearest and next-nearest neighbors exchange interactions. The magnetic surface can be considered as semi-infinite slabs at the end of the bulk structures. The breakdown of translation symmetry, in the normal direction of the surface, gives rise to localized spin wave modes in its neighborhood. The localized magnon spectrum is derived as elements of a Landauer-type scattering matrix, in the three cubic lattices sc, bcc and fcc. The magnon properties are simulated and determined numerically for each cubic lattice by using the matching technique. The observed fluctuations in the numerical results demonstrate the interference magnon effects between scattered spinwaves and the localized magnon states, generated by the surface region with characteristic Fano resonances. In cubic leads, the localized surface spin states are sensitive to the local magnetic coupling and the incident direction in the surface boundary. In this contribution, the normalized energy of spinwaves arising from the absence of translation symmetry is analyzed for each cubic system as a function of the exchange integral parameters. This addresses the dependence of the surface magnon on the different possibilities of the of the exchange parameters variation from softening to hardening in the neighborhood of the surface region.