  Multi-peaked localized stationary solutions of the discrete nonlinear
Schrodinger (DNLS) equation are presented in one (1D) and two (2D) dimensions.
These are excited states of the discrete spectrum and correspond to
multi-breather solutions. A simple, very fast, and efficient numerical method,
suggested by Aubry, has been used for their calculation. The method involves no
diagonalization, but just iterations of a map, starting from trivial solutions
of the anti-continuous limit. Approximate analytical expressions are presented
and compared with the numerical results. The linear stability of the calculated
stationary states is discussed and the structure of the linear stability
spectrum is analytically obtained for relatively large values of nonlinearity.
