**LA_PBSV** computes the solution to a linear system of equations
, where has band form and is
real symmetric or complex Hermitian and, in either case, positive definite,
and where and are rectangular matrices or vectors.
The Cholesky decomposition is used to factor as

where is an upper triangular band matrix and is a lower triangular band matrix, each with the same number of superdiagonals or subdiagonals as . The factored form of is then used to solve the above system.

Susan Blackford 2001-08-19