An iterative method for the calculation of narrow solitary excitations on atomic chains