Appendix E. THOMAS ALGORITHM FOR MULTI-LAYERS OF VG MODEL


Thomas algorithm solves a tridiagonal matrix equation:

into a form:

where ai's, bi's, ci's, and fi's are give constants.bi's and yi'$s can be calculated by comparing corresponding elements of the above two equations, and xi's can be expressed:

This equation can be solved by assuming: 

and obtaining ai's, bi's, and yi's incrementally: 

Then xi's can be obtained by a back substitution:


 



Last modified: Oct 15, 1999
VG Model / Samuel Lee / VADOSE.NET