**LA_GBSV** computes the solution to a real or complex linear system of
equations , where is a square band matrix and and
are rectangular matrices or vectors. The decomposition with row
interchanges is used to factor as , where is a product
of permutation and unit lower triangular matrices with subdiagonals,
and is upper triangular with superdiagonals.
The factored form of is then used to solve the above system.

Susan Blackford 2001-08-19