An efficient breakdown-free algorithm for numerically evaluating the determinants of (p, q)-pentadiagonal matrices

(p,q)-Pentadiagonal matrices have received considerable attention in recent years, which are a generalization of pentadiagonal matrices. In this paper, a breakdown-free algorithm is presented for numerically evaluating the determinants of n-by-n (p,q)-pentadiagonal matrices. The algorithm is based on the use of a reliable tridiagonalization process which preserves the banded structure and sparsity of the original matrix. Numerical examples are given in order to illustrate the effectiveness of the proposed algorithm. All of the experiments are performed on a computer with the aid of programs written in MATLAB.