For non-classically damped structures subjected to evolutionary random seismic excitations, the non-stationary random responses are computed by means of a high precision direct (HPD) integration scheme combined with the pseudo excitation method. Only real modes are used, so that the reduced equations of motion remain coupled for such non-classically damped structures. In the given examples, the efficiency of this method is compared with that of the Newmark method.