A recursive decoupling algorithm for solving banded linear systems

A recursive bisection method for the solution of banded matrix systems of linear equations which occur frequently in the finite difference/element discretisations of boundary value problems is derived. We analyse the proposed method in detail for tridiagonal systems but the algorithm can be readily extended in block form to matrices of bandwidth (2s +1).