A method is proposed for generating multivariate time series which are
required to satisfy a given spectral density function, which extends
previous work on univariate time series. The performance of the method
is assessed in a small simulation exercise for a bivariate long memor
y model and is found to perform well.