Use LU factorization of double precision positive definite band matrix to compute determinant. (No provision for inverse.)
CALL DPBDI(ABD, LDA, N, M, DET)