A structure preserving matrix factorization for solving general periodic pentadiagonal Toeplitz linear systems

In the current article, a new algorithm for solving linear systems with the periodic pentadiagonal Toeplitz (PPT) coefficient matrix is presented. This algorithm is based on a structure preserving matrix factorization and Sherman–Morrison–Woodbury formula. Some numerical examples are given in order to illustrate the effectiveness of the proposed algorithms. All of the experiments are performed on a computer with the aid of programs written in MATLAB.