Frequently in signal processing one is faced with situations where a large system of linear equations, with a Toeplitz or a Hankel coefficient matrix, needs to be solved. One efficient way of solving these kinds of equations is by Levinson recursion. The Levinson recursion does not require explicit storage of the Toeplitz (or Hankel) coefficient matrix and the number of multiplies required is proportional to the square of the number of unknowns.