We present a method for the multiplication of an arbitrary vector by a symm
etric centrosymmetric matrix, requiring 5/4n(2) + O(n) floating-point opera
tions, rather than the 2n(2) operations needed in the case of an arbitrary
matrix. Combining this method with Trench's algorithm for Toeplitz matrix i
nversion yields a method for solving Toeplitz systems with the same complex
ity as Levinson's algorithm. (C) 2000 Elsevier Science Inc. All rights rese
rved. AMS classification: 65F05.