Some Convalescent Methods for the Solution of Systems of Linear Equations
Applied and Computational Mathematics
Volume 4, Issue 3, June 2015, Pages: 207-213
Received: May 16, 2015;
Accepted: May 26, 2015;
Published: Jun. 9, 2015
Views 4258 Downloads 168
M. Rafique, Department of Mathematics, Faculty of Science, HITEC University, Taxila, Pakistan
Sidra Ayub, Department of Mathematics, Faculty of Science, HITEC University, Taxila, Pakistan
Follow on us
In a variety of problems in the fields of physical sciences, engineering, economics, etc., we are led to systems of linear equations, Ax = b, comprising n linear equations in n unknowns x1, x2, …, xn, where A = [aij] is an nxn coefficient matrix, and x = [x1 x2 . . .xn]T, b = [b1 b2 . . .bn]T are the column vectors. There are many analytical as well as numerical methods[1}–  to solve such systems of equations, including Gauss elimination method, and its modifications namely Doolittle’s method, Crout’s method and Cholesky’s method, which employ LU-decomposition method, where L = [iij] and u = [uij] are the lower and upper triangular matrices respectively. The LU-decomposition method was first introduced by the mathematician Alan M. Turing- in 1948. Here, in this paper we have made an effort to modify the existing LU-decomposition methods to solve the above mentioned system Ax = b, with the least possible endeavour. It may be seen that the Gauss elimination method, , ,  needs about 2n3/3 operations, while Doolittle’s and Crout’s methods require n2 operations. Accordingly, in these methods we are required to evaluate n2 number of unknown elements of the L and U matrices. Moreover, Cholesky’s method requires 2n2/3 operations. Accordingly this method requires evaluation of 2n2/3 number of unknown elements of the L and U matrices But, in contrast, the improved Doolittle’s, Crout’s and Cholesky’s methods presented in this paper require evaluation of only (n–1)2 number of unknown elements of the L and U matrices. Moreover, an innovative method is also presented in this paper which requires evaluation of even less number of unknown elements of the L and U matrices. In this method we need to evaluate only (n–2)2 number of the said unknown elements. Thus, by employing these methods, the computational time and effort required for the purpose can substantially be reduced.
System of Equations, Matrix, Column Vector, Decomposition, Doolittle’s Method, Crout’s Method, Cholesky’s Method
To cite this article
Some Convalescent Methods for the Solution of Systems of Linear Equations, Applied and Computational Mathematics.
Vol. 4, No. 3,
2015, pp. 207-213.
E. Kreyszig: Advanced Engineering Mathematics, John Wiley, (2011)
A. M. Turing: Rounding-Off Errors in Matrix Processes”, The Quarterly Journal of Mechanics and Applied Mathematics 1: 287-308.doi:10.1093/qjmam/1.1.287 (1948)
Myrick H. Doolittle:. Method employed in the solution of normal equations and the adjustment of a triangulation, U.S. Coast and Geodetic Survey Report, Appendix 8, Paper No. 3, blz.. 115–120.( 1878)
David Poole: Linear algebra: A Modern Introduction (2nd Edition), Thomson Brooks/Cole ISSN 0-534-99845-3 (2006)
James R. Bunch, and John Hopkins: Triangular factorization and inversion by fast matrix multiplication, Mathematics of Computation, 28, 231-236, doi:10.2307/2005828, ISSN 0025-5718
J H Wilkinson: The Algebraic Eigenvalue Problem, Oxford University Press (1988)
B. Faires, Numerical Analysis, PWS Pub., Boston, (1993).
G.H.Golub and C.F.Van Loan, Matrix Computations, John Hopkins, Baltimore, (1989).
J.H. Mathews and K. D. Fink: Numerical Methods Using MATLAB (4th Edition)
L. V. Fausett: Applied Numerical Analysis Using MATLAB, Pearson, (2009)
S.C.Chapra and R.P.Canale, Numerical Methods for Engineers, Mc-Graw-Hill, New York,(1990).
Petre Teodorescu, Nicolae-Doru Stanescu, Nicolae Pandrea: Numerical Analysis with Applications in Mechanics and Engineering, John Wiley (2013).