nag_dgbsv (f07bac) uses the LU decomposition with partial pivoting and row interchanges to factor A as A=PLU, where P is a permutation matrix, L is a product of permutation and unit lower triangular matrices with kl subdiagonals, and U is upper triangular with kl+ku superdiagonals.