Example 3. Given  [Graphics:Images/LUFactorMod_gr_39.gif].  Can A be factored  A = LU?    

Solution 3.

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

Use the LandU subroutine and see what happens.

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



[Graphics:../Images/LUFactorMod_gr_42.gif]
[Graphics:../Images/LUFactorMod_gr_43.gif]
[Graphics:../Images/LUFactorMod_gr_44.gif]

[Graphics:../Images/LUFactorMod_gr_45.gif]
[Graphics:../Images/LUFactorMod_gr_46.gif]
[Graphics:../Images/LUFactorMod_gr_47.gif]

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


We get a division by zero error, hence there is no LU-decomposition of A.

Remark.  The matrix A is nonsingular.  

We are done.  This shows that nonsingular is not sufficient to guarantee an  LU-decomposition.

Aside.  We can interchange rows and the new matrix B will have an   LU-decomposition.

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




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

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

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

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

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

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




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



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

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

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

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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004