For long-memory time series, we show that the Toeplitz system sect;n(f)x = b can be solved inO(n log5=2 n) operations using a well-known version of the preconditioned conjugate gradient method, where sect;n(f) is the npound;n covariance matrix, f is the spectral density and b is a known vector....