Example 1 (c).  Find the A = LU  factorization for the matrix  [Graphics:Images/CholeskyMod_gr_54.gif].  Use the Cholesky method.  

Solution 1 (c).

Enter the matrix.

 

[Graphics:../Images/CholeskyMod_gr_83.gif]


Invoke the subroutine Cholesky.  

[Graphics:../Images/CholeskyMod_gr_84.gif]



[Graphics:../Images/CholeskyMod_gr_85.gif]

[Graphics:../Images/CholeskyMod_gr_86.gif]

[Graphics:../Images/CholeskyMod_gr_87.gif]

[Graphics:../Images/CholeskyMod_gr_88.gif]


Verify the factorization.

[Graphics:../Images/CholeskyMod_gr_89.gif]



[Graphics:../Images/CholeskyMod_gr_90.gif]

[Graphics:../Images/CholeskyMod_gr_91.gif]

[Graphics:../Images/CholeskyMod_gr_92.gif]

[Graphics:../Images/CholeskyMod_gr_93.gif]

[Graphics:../Images/CholeskyMod_gr_94.gif]

This time the elements on the diagonals of  L  and  U  are the same.  

Remark.  The Cholesky method is used only when A is a real, symmetric and positive definite matrix.  

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004