Abstract
This paper describes an algorithm for the approximate solution of vector Toeplitz systems via an iterative-improvement method. The algorithm exploits the special structure of Toeplitz matrices, namely, their similarity to vector circulants, and is particularly well suited for solving large systems. Sufficient convergence conditions and concrete error bounds for the iteration are presented along with an application of the routine to a problem in the design of planar digital filters for image processing.