In this paper, for a prescribed real spectrum, using properties of the circ
ulant matrices and of the symmetric persymmetric matrices, we derive a fast
and stable algorithm to construct a symmetric nonnegative matrix which rea
lizes the spectrum. The algorithm is based on the fast Fourier transform. (
C) 2001 Elsevier Science Ltd. All rights reserved.