Example 2.  Solve the 31 by 31 system    in Example 1 using our TriDiagonal[a,d,c,b] subroutine.  Compute the solution using decimal.  

Solution 2 (a).

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

Mathematica is able to get a slightly better answer.

(c) John H. Mathews 2004