We describe a method for obtaining an analytic form for an inverse of a class of symmetric semi-infinite banded matrices, which are, apart from a finite number of terms, of the Toeplitz type. The results are applied to the determination of the spectrum of two-magnon excitations in Heisenberg spin chains with next-nearest-neighbour interactions.