A computationally convenient method is presented for simulating a sequence of sampled values of a stationary Gaussian process having a specified correlation function or power spectrum.

-transform theory is applied to provide a simple recursive formula for generating the values from a set of independent Gaussian random variables.