221
Views
6
CrossRef citations to date
0
Altmetric
Section B

A fast numerical algorithm for solving nearly penta-diagonal linear systems

, &
Pages 851-860 | Received 14 Oct 2011, Accepted 04 Feb 2012, Published online: 13 Mar 2012
 

Abstract

In this paper, we present a fast numerical algorithm for solving nearly penta-diagonal linear systems and show that the computational cost is less than those of three algorithms in El-Mikkawy and Rahmo, [Symbolic algorithm for inverting cyclic penta-diagonal matrices recursively–Derivation and implementation, Comput. Math. Appl. 59 (2010), pp. 1386–1396], Lv and Le [A note on solving nearly penta-diagonal linear systems, Appl. Math. Comput. 204 (2008), pp. 707–712] and Neossi Nguetchue and Abelman [A computational algorithm for solving nearly penta-diagonal linear systems, Appl. Math. Comput. 203 (2008), pp. 629–634.]. In addition, an efficient way of evaluating the determinant of a nearly penta-diagonal matrix is also discussed. The algorithm is suited for implementation using computer algebra systems (CAS) such as MATLAB, MACSYMA and MAPLE. Some numerical examples are given in order to illustrate the efficiency of our algorithm.

2010 AMS Subject Classification :

Acknowledgements

This work was supported by the Natural Science Foundation of China (NSFC) under grant 11071192 and the International Science and Technology Cooperation Program of China under grant 2010DFA14700. The authors would like to thank the referees who substantially enhanced the quality of the paper.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.