Example 5  Test the matrix A of Examples 2 and 4 to see if A is strictly diagonally dominant.

Solution 5.

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


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

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

We are done !

Aside. We can use Mathematica's Eigensystem subroutine and actually find the spectral radius of  [Graphics:../Images/GaussSeidelMod_gr_228.gif]  is less than 1.    This is just for fun.

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

[Graphics:../Images/GaussSeidelMod_gr_230.gif]
[Graphics:../Images/GaussSeidelMod_gr_231.gif]
[Graphics:../Images/GaussSeidelMod_gr_232.gif]

Since the spectral radius is  [Graphics:../Images/GaussSeidelMod_gr_233.gif],  iteration in Example 4 will not converge.

Aside.  Test the matrix of Example 1.

This is just for fun.

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

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

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

Aside.  Test the matrix of Example 1to see if the spectral radius of  [Graphics:../Images/GaussSeidelMod_gr_237.gif]  is less than 1.  

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


[Graphics:../Images/GaussSeidelMod_gr_239.gif]
[Graphics:../Images/GaussSeidelMod_gr_240.gif]
[Graphics:../Images/GaussSeidelMod_gr_241.gif]

Since the spectral radius is  [Graphics:../Images/GaussSeidelMod_gr_242.gif],  iteration in Example 1 will converge.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004