Example 2.  Solve the 31 by 31 system  [Graphics:Images/Tri-DiagonalMod_gr_90.gif]  in Example 1 using our TriDiagonal[a,d,c,b] subroutine.  Compute the solution using decimal arithmetic and rational arithmetic.  

Solution 2 (b).

[Graphics:../Images/Tri-DiagonalMod_gr_101.gif]

[Graphics:../Images/Tri-DiagonalMod_gr_102.gif]

[Graphics:../Images/Tri-DiagonalMod_gr_103.gif]

[Graphics:../Images/Tri-DiagonalMod_gr_104.gif]

[Graphics:../Images/Tri-DiagonalMod_gr_105.gif]

[Graphics:../Images/Tri-DiagonalMod_gr_106.gif]

[Graphics:../Images/Tri-DiagonalMod_gr_107.gif]


[Graphics:../Images/Tri-DiagonalMod_gr_108.gif]

We can compare this answer with the previous one in Example 1.

[Graphics:../Images/Tri-DiagonalMod_gr_109.gif]


[Graphics:../Images/Tri-DiagonalMod_gr_110.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004