Under the action of stationary random ground motion any power spectral density matrix of structural responses can be produced through the multiplication of two response vectors which are easily obtained by means of a harmonic response analysis; it is a CQC algorithm with high efficiency and simple procedure. The non-orthogonal damping matrix can be dealt with conveniently and accurately. Two examples are given to help understand the procedure of the algorithm.