<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ALAMT</journal-id><journal-title-group><journal-title>Advances in Linear Algebra &amp; Matrix Theory</journal-title></journal-title-group><issn pub-type="epub">2165-333X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/alamt.2015.52004</article-id><article-id pub-id-type="publisher-id">ALAMT-56990</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Inverse Nonnegativity of Tridiagonal &lt;i&gt;M&lt;/i&gt;-Matrices under Diagonal Element-Wise Perturbation
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ohamed</surname><given-names>A. Ramadan</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mahmoud</surname><given-names>M. Abu Murad</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Mathematics, Faculty of Science, Menoufia University, Shebeen El-Kom, Egypt</addr-line></aff><aff id="aff2"><addr-line>Department of Mathematics, Faculty of Science, Zagazig University, Ash Sharqiyah, Egypt</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>ramadanmohamed13@yahoo.com(OAR)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>09</day><month>06</month><year>2015</year></pub-date><volume>05</volume><issue>02</issue><fpage>37</fpage><lpage>45</lpage><history><date date-type="received"><day>24</day>	<month>March</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>6</month>	<year>June</year>	</date><date date-type="accepted"><day>9</day>	<month>June</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  One of the most important properties of M-matrices is element-wise non-negative of its inverse. In this paper, we consider element-wise perturbations of tridiagonal M-matrices and obtain bounds on the perturbations so that the non-negative inverse persists. The largest interval is given by which the diagonal entries of the inverse of tridiagonal M-matrices can be perturbed without losing the property of total nonnegativity. A numerical example is given to illustrate our findings.
 
</p></abstract><kwd-group><kwd>Totally Positive Matrix</kwd><kwd> Totally Nonnegative Matrix</kwd><kwd> Tridiagonal Matrices</kwd><kwd> Compound Matrix</kwd><kwd>  Element-Wise Perturbations</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>In many mathematical problems, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x5.png" xlink:type="simple"/></inline-formula>-matrices and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x6.png" xlink:type="simple"/></inline-formula>-matrices play an important role. It is often useful to know the properties of their inverses, especially when the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x7.png" xlink:type="simple"/></inline-formula>-matrices and the M-matrices have a special combinatorial structure, for more details we refer the reader [<xref ref-type="bibr" rid="scirp.56990-ref1">1</xref>] . M-matrices have important applications, for instance, in iterative methods, in numerical analysis, in the analysis of dynamical systems, in economics, and in mathematical programming. One of the most important properties of some kinds of M-matrices is the nonegativity of their inverses, which plays central role in many of mathematical problems.</p><p>An <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x8.png" xlink:type="simple"/></inline-formula> real matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x9.png" xlink:type="simple"/></inline-formula> is called M-matrix if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x10.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x11.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x12.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x13.png" xlink:type="simple"/></inline-formula>, over the years, M-matrices have considerable attention, in large part because they arise in many applications [<xref ref-type="bibr" rid="scirp.56990-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.56990-ref3">3</xref>] . Recently, a noticeable amount of attention has turned to the inverse of tridiagonal M-matrices (those matrices which happen to be inverses of special form of M-matrices with property <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x14.png" xlink:type="simple"/></inline-formula> whenever<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x15.png" xlink:type="simple"/></inline-formula>) and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x16.png" xlink:type="simple"/></inline-formula> is generalized strictly diagonally dominant. A matrix is said to be generalized (strictly) diagonally dominant</p><p>if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x17.png" xlink:type="simple"/></inline-formula>. Of particular importance to us is the fact that since <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x18.png" xlink:type="simple"/></inline-formula> is an M-matrix it is non-singular and</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x19.png" xlink:type="simple"/></inline-formula>&gt;0where the inequality is satisfied element-wise. A rich class of M-matrices were introduced by Ostrowski in 1937 [<xref ref-type="bibr" rid="scirp.56990-ref4">4</xref>] , with reference to the work of Minkowski [<xref ref-type="bibr" rid="scirp.56990-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.56990-ref6">6</xref>] . A condition which is easy to check is that a matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x20.png" xlink:type="simple"/></inline-formula> is an M-matrix if and only if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x21.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x22.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x23.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x24.png" xlink:type="simple"/></inline-formula> is generalized strictly diagonally dominant.</p><p>In this paper, we consider the inverse of perturbed M-matrix. Specifically we consider the effect of changing single elements inside the diagonal of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x25.png" xlink:type="simple"/></inline-formula>. We are interested in the large amount by which the single diagonal element of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x26.png" xlink:type="simple"/></inline-formula> can be varied without losing the property of total nonnegativity.</p><p>The reminder of the paper is organized as follows. In section 2, we explain our notations and some needed important definitions are presented. In section 3, some auxiliary results and important prepositions and lemmas are stated. In section 4, we present our results.</p></sec><sec id="s2"><title>2. Notations</title><p>In this section we introduce the notation that will be used in developing the paper. For <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula> we denote by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula> the set of all strictly increasing sequences of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula> integers chosen from<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula>. For<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula>, we denote by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula> the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula> submatrix of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula> contained in the rows indexed by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula> and columns indexed by<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x37.png" xlink:type="simple"/></inline-formula>. A matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x38.png" xlink:type="simple"/></inline-formula> is called totally positive (abbreviated TP henceforth) if <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x39.png" xlink:type="simple"/></inline-formula> and totally nonnegative (abbreviated TN) if <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x40.png" xlink:type="simple"/></inline-formula> for all<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x41.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x42.png" xlink:type="simple"/></inline-formula>. For a given index<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x43.png" xlink:type="simple"/></inline-formula>, with property<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x44.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x45.png" xlink:type="simple"/></inline-formula>, the dispersion of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x46.png" xlink:type="simple"/></inline-formula>,</p><p>denoted by<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x47.png" xlink:type="simple"/></inline-formula>, is defined to be<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x48.png" xlink:type="simple"/></inline-formula>.</p><p>Throughout this paper we use the following notation for general tridiagonal M-matrix:</p><disp-formula id="scirp.56990-formula49"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x49.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x50.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x51.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x52.png" xlink:type="simple"/></inline-formula>, and each <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x53.png" xlink:type="simple"/></inline-formula> is large enough that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x54.png" xlink:type="simple"/></inline-formula> is strictly diagonally dominant.</p><p>We let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x55.png" xlink:type="simple"/></inline-formula> to be the square standard basis matrix whose only nonzero entry is 1 that occurs in the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x56.png" xlink:type="simple"/></inline-formula> position.</p><p>Definition 2.1 Compound Matrices ([<xref ref-type="bibr" rid="scirp.56990-ref7">7</xref>] , p. 19).</p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x57.png" xlink:type="simple"/></inline-formula> be a square matrix of order<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x58.png" xlink:type="simple"/></inline-formula>. Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x59.png" xlink:type="simple"/></inline-formula> be the index set of cardinality<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x60.png" xlink:type="simple"/></inline-formula>, defining<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x60.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x61.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x60.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x62.png" xlink:type="simple"/></inline-formula>are the index sets of cardinality<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x60.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x63.png" xlink:type="simple"/></inline-formula>.</p><p>Construct the following table which depends on<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x64.png" xlink:type="simple"/></inline-formula>.</p><p>The created matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x89.png" xlink:type="simple"/></inline-formula></p><p>is called<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x90.png" xlink:type="simple"/></inline-formula>, compound matrix of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x90.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x91.png" xlink:type="simple"/></inline-formula>.</p><p>For example, if</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x92.png" xlink:type="simple"/></inline-formula>with indexed sets<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x93.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x94.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x95.png" xlink:type="simple"/></inline-formula>.</p><p>Then<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x96.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s3"><title>3. Auxiliary Results</title><p>We start with some basic facts on tridiagonal M-matrices. We can find the determinant of any <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x97.png" xlink:type="simple"/></inline-formula> tridiagonal M-matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x98.png" xlink:type="simple"/></inline-formula> by using the following recursion equation [<xref ref-type="bibr" rid="scirp.56990-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.56990-ref9">9</xref>] .</p><disp-formula id="scirp.56990-formula50"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x99.png"  xlink:type="simple"/></disp-formula><p>And we have the following proposition for finding the determinant of a <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x100.png" xlink:type="simple"/></inline-formula> tridiagonal M-matrix<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x100.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x101.png" xlink:type="simple"/></inline-formula>.</p><p>Proposition 3.1 ([<xref ref-type="bibr" rid="scirp.56990-ref10">10</xref>] , formula 4.1) For any <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x102.png" xlink:type="simple"/></inline-formula> tridiagonal M-matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x102.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x103.png" xlink:type="simple"/></inline-formula> the following relation is true</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x104.png" xlink:type="simple"/></inline-formula>.</p><p>We will present now some of propositions of nonsingular totally nonnegative matrices which important for our work.</p><p>Proposition 3.2 [<xref ref-type="bibr" rid="scirp.56990-ref10">10</xref>] [<xref ref-type="bibr" rid="scirp.56990-ref11">11</xref>]</p><p>For any nonsingular totally nonnegative matrix<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x105.png" xlink:type="simple"/></inline-formula>, all principle minors are positive.</p><p>That is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x106.png" xlink:type="simple"/></inline-formula>, for all <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x107.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x107.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x108.png" xlink:type="simple"/></inline-formula>.</p><p>Proposition 3.3 ([<xref ref-type="bibr" rid="scirp.56990-ref7">7</xref>] , p. 21)</p><p>Let M be a nonsingular tridiagonal M-matrix, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x109.png" xlink:type="simple"/></inline-formula> be the inverse of the matrix M then</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x110.png" xlink:type="simple"/></inline-formula>, when.</p><p>In the sequel we will make use the following lemma, see, e.g. [<xref ref-type="bibr" rid="scirp.56990-ref12">12</xref>] .</p><p>Lemma 3.4 (Sylvester Identity)</p><p>Partition square matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x112.png" xlink:type="simple"/></inline-formula> of order n, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x112.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x113.png" xlink:type="simple"/></inline-formula>, as:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x114.png" xlink:type="simple"/></inline-formula>,</p><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x115.png" xlink:type="simple"/></inline-formula> square matrix of order <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x115.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x116.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x115.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x116.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x117.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x115.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x116.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x117.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x118.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x115.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x116.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x117.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x118.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x119.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x115.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x116.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x117.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x118.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x119.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x120.png" xlink:type="simple"/></inline-formula> are scalars.</p><p>Define the submatrices</p><disp-formula id="scirp.56990-formula51"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x121.png"  xlink:type="simple"/></disp-formula><p>If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x122.png" xlink:type="simple"/></inline-formula> is nonsingular, then</p><disp-formula id="scirp.56990-formula52"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x123.png"  xlink:type="simple"/></disp-formula><p>Lemma 3.5 ([<xref ref-type="bibr" rid="scirp.56990-ref11">11</xref>] , p.199) Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x124.png" xlink:type="simple"/></inline-formula> be a square matrix of order<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x125.png" xlink:type="simple"/></inline-formula>, with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x126.png" xlink:type="simple"/></inline-formula>. Then</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x127.png" xlink:type="simple"/></inline-formula>is totally nonnegative<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x127.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x128.png" xlink:type="simple"/></inline-formula>.</p><p>We now state an important result which links the determinant of M-matrix with the value of the elements of its inverse.</p><p>Lemma 3.6 [<xref ref-type="bibr" rid="scirp.56990-ref10">10</xref>] Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x129.png" xlink:type="simple"/></inline-formula> be a tridiagonal matrix of order n, then we can find the elements of in-</p><p>verse matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x130.png" xlink:type="simple"/></inline-formula> by using the following formula</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x131.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s4"><title>4. Main Results</title><p>In this section, we present our results based on the inverse of tridiagonal M-matrices. Firstly we begin with the following theorem.</p><p>Theorem 4.1</p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x132.png" xlink:type="simple"/></inline-formula> be strictly diagonally dominant M-matrix.</p><p>If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x133.png" xlink:type="simple"/></inline-formula> is the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x134.png" xlink:type="simple"/></inline-formula> compound matrix of M then the matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x133.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x134.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x135.png" xlink:type="simple"/></inline-formula> is totally nonnega-</p><p>tive matrix. Moreover, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x136.png" xlink:type="simple"/></inline-formula>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x136.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x137.png" xlink:type="simple"/></inline-formula></p><p>Proof: Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x138.png" xlink:type="simple"/></inline-formula> be strictly diagonally dominant M-matrix.</p><p>Then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x139.png" xlink:type="simple"/></inline-formula> is totally nonnegative matrix. So is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x139.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x140.png" xlink:type="simple"/></inline-formula>.</p><p>You can find this formula in ([<xref ref-type="bibr" rid="scirp.56990-ref7">7</xref>] , p. 21).</p><p>There is an explicit formula for the determinant of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x141.png" xlink:type="simple"/></inline-formula> given as</p><disp-formula id="scirp.56990-formula53"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x142.png"  xlink:type="simple"/></disp-formula><p>Multiply the first row by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x143.png" xlink:type="simple"/></inline-formula> and add it to the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x143.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x144.png" xlink:type="simple"/></inline-formula> row to obtain</p><disp-formula id="scirp.56990-formula54"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x145.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x146.png" xlink:type="simple"/></inline-formula>,</p><p>And now apply an induction argument to get the result.</p><p>Numerical Example: Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x147.png" xlink:type="simple"/></inline-formula> be strictly diagonally dominant M-matrix, then</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x148.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x148.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x149.png" xlink:type="simple"/></inline-formula> is totally nonnegative.</p><p>Note that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x150.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x150.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x152.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x150.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x152.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x151.png" xlink:type="simple"/></inline-formula></p><p>Numerically we can conclude the following fact.</p><p>Fact: For any <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x153.png" xlink:type="simple"/></inline-formula> tridiagonal M-matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x153.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x154.png" xlink:type="simple"/></inline-formula> the following formula is true.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x155.png" xlink:type="simple"/></inline-formula>for</p><p>Moreover, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x157.png" xlink:type="simple"/></inline-formula></p><p>To prove this result we use Theorem 4.1.</p><p>Suppose M is nonsingular then<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x158.png" xlink:type="simple"/></inline-formula>, so</p><disp-formula id="scirp.56990-formula55"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x159.png"  xlink:type="simple"/></disp-formula><p>For example, when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x160.png" xlink:type="simple"/></inline-formula>, the M-matrix of our form</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x161.png" xlink:type="simple"/></inline-formula>has an inverse given as</p><disp-formula id="scirp.56990-formula56"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x162.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56990-formula57"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x163.png"  xlink:type="simple"/></disp-formula><p>Similarly we can find<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x164.png" xlink:type="simple"/></inline-formula>.</p><p>Illustrative Example: Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x165.png" xlink:type="simple"/></inline-formula> be a tridiagonal M-matrix and</p><disp-formula id="scirp.56990-formula58"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x166.png"  xlink:type="simple"/></disp-formula><p>Note that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x167.png" xlink:type="simple"/></inline-formula></p><p>Observe that the error came from the rounded to the nearest part of 10,000.</p><p>Theorem 4.2 Let M be a strictly diagonally dominant M-matrix, if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x168.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x168.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x169.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x168.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x169.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x170.png" xlink:type="simple"/></inline-formula>then</p><disp-formula id="scirp.56990-formula59"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x171.png"  xlink:type="simple"/></disp-formula><p>Proof:</p><p>Assume <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x172.png" xlink:type="simple"/></inline-formula> and</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x173.png" xlink:type="simple"/></inline-formula>, , ,.</p><p>Note that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x177.png" xlink:type="simple"/></inline-formula> by previous fact, and by using Sylvester's identity, we have</p><disp-formula id="scirp.56990-formula60"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x178.png"  xlink:type="simple"/></disp-formula><p>Moreover we conclude the following theorem.</p><p>Theorem 4.3 Let M be the M-matrix defined above then</p><disp-formula id="scirp.56990-formula61"><graphic  xlink:href="http://html.scirp.org/file/1-2230076x179.png"  xlink:type="simple"/></disp-formula><p>For example</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x180.png" xlink:type="simple"/></inline-formula>,</p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x182.png" xlink:type="simple"/></inline-formula> then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x182.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x183.png" xlink:type="simple"/></inline-formula></p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x184.png" xlink:type="simple"/></inline-formula>and</p><p>Now, we will perturb elements inside the diagonal band of the inverse of M-matrix without losing the nonnegativity property. We begin with the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x186.png" xlink:type="simple"/></inline-formula> element then generalize to other elements.</p><p>Theorem 4.4 Let M be a strictly diagonally dominant tridiagonal <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x187.png" xlink:type="simple"/></inline-formula>-matrix. Then the matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x187.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x188.png" xlink:type="simple"/></inline-formula></p><p>is totally nonnegative for all<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x189.png" xlink:type="simple"/></inline-formula>.</p><p>Proof:</p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x190.png" xlink:type="simple"/></inline-formula></p><p>Be a nonsingular strictly diagonally dominant tridiagonal M-matrix then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x191.png" xlink:type="simple"/></inline-formula> is totally nonnegative.</p><p>By Lemma 3.5 and Proposition 3.2, we have</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x192.png" xlink:type="simple"/></inline-formula>is totally nonnegative<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x192.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x193.png" xlink:type="simple"/></inline-formula>.</p><p>By using the formula in Proposition 3.3</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x194.png" xlink:type="simple"/></inline-formula>.</p><p>Note that a similar result holds for decreasing the element <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x195.png" xlink:type="simple"/></inline-formula> by considering the matrix<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x195.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x196.png" xlink:type="simple"/></inline-formula>, which reverses the matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x195.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x196.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x197.png" xlink:type="simple"/></inline-formula> as the relation<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x195.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x196.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x197.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x198.png" xlink:type="simple"/></inline-formula>.</p><p>We can generalize this result for the other elements of diagonal.</p><p>Theorem 4.5 Assume M is a strictly diagonally dominant tridiagonal M-matrix. Then the matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x199.png" xlink:type="simple"/></inline-formula></p><p>is totally nonnegative for all<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x200.png" xlink:type="simple"/></inline-formula>.</p><p>Proof: Suppose that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x201.png" xlink:type="simple"/></inline-formula> is not totally nonnegative for all<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x201.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x202.png" xlink:type="simple"/></inline-formula>, then there exist</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x203.png" xlink:type="simple"/></inline-formula>both contain <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x203.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x204.png" xlink:type="simple"/></inline-formula> such that<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x203.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x204.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x205.png" xlink:type="simple"/></inline-formula>.</p><p>To compute <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x206.png" xlink:type="simple"/></inline-formula> expand the determinant along the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x206.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x207.png" xlink:type="simple"/></inline-formula> row of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x206.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x207.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x208.png" xlink:type="simple"/></inline-formula> then</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x209.png" xlink:type="simple"/></inline-formula>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x209.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x210.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x209.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x210.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x211.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x209.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x210.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x211.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x212.png" xlink:type="simple"/></inline-formula> is some minor of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x209.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x210.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x211.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x212.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x213.png" xlink:type="simple"/></inline-formula>.</p><p>Take the case when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x214.png" xlink:type="simple"/></inline-formula> odd. Thus, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x214.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x215.png" xlink:type="simple"/></inline-formula>is a positive linear compination of minors of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x214.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x215.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x216.png" xlink:type="simple"/></inline-formula> and hence is positive, which contradicts the assumption.</p><p>Now suppose <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x217.png" xlink:type="simple"/></inline-formula> such that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x217.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x218.png" xlink:type="simple"/></inline-formula> is totally nonnegative matrix.</p><p>Suppose that<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x219.png" xlink:type="simple"/></inline-formula>, then by Theorem 4.3.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x220.png" xlink:type="simple"/></inline-formula>, since</p><p>which contradicts the nonnegativity of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x222.png" xlink:type="simple"/></inline-formula>.</p><p>Numerical Example: Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x223.png" xlink:type="simple"/></inline-formula> is strictly diagonally dominant tridiagonal M-matrix</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x224.png" xlink:type="simple"/></inline-formula>so</p><p>The matrices</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x225.png" xlink:type="simple"/></inline-formula>, ,</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x227.png" xlink:type="simple"/></inline-formula>, ,</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x229.png" xlink:type="simple"/></inline-formula>, , and</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x231.png" xlink:type="simple"/></inline-formula>,.</p><p>are TNN matrices.</p><p>Note that<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-2230076x233.png" xlink:type="simple"/></inline-formula>.</p></sec></body><back><ref-list><title>References</title><ref id="scirp.56990-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">McDonald, J.J., Nabben, R., Neumannand, M., Schneider, H. and Tsatsomeros, M.J. (1998) Inverse Tridiagonal Z-Matrices. Linear and Multilinear Algebra, 45, 75-97. http://dx.doi.org/10.1080/03081089808818578</mixed-citation></ref><ref id="scirp.56990-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Berman, A. and Plemmons, R. (1979) Nonnegative Matrices in the Mathenlatical Sciences. Academic, New York.</mixed-citation></ref><ref id="scirp.56990-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Pena, J.M. (1995) M-Matrices Whose Inverses Are Totally Positive. Linear Algebra and Its Applications, 221, 189-193. http://dx.doi.org/10.1016/0024-3795(93)00244-T</mixed-citation></ref><ref id="scirp.56990-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Ostrowski, A. (1937) über die Determinanten mit überwiegender Hauptdiagonale. Commentarii Mathematici Helvetici, 10, 69-96. http://dx.doi.org/10.1007/BF01214284</mixed-citation></ref><ref id="scirp.56990-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Minkowski, H. (1900) Zur Theorie der Einheiten in den algebraischen Zahlkorper. Nachrichten von der Konigl. Gesellschaft der Wissenschaften und der Georg-Augusts-Universitat zu Gottingen Mathematisch-Physikalische Klasse: Fachgruppe II, Nachrichten aus der Physik, Astronomie, Geophysik, Technik, 90-93.</mixed-citation></ref><ref id="scirp.56990-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Minkowski, H. (1907) Diophantische Approximationen. Tuebner, Leipzig. http://dx.doi.org/10.1007/978-3-663-16055-7</mixed-citation></ref><ref id="scirp.56990-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Horn, R. and Johnson, C. (1991) Topics in Matrix Analysis. Cambridge University Press, Cambridge.http://dx.doi.org/10.1017/CBO9780511840371</mixed-citation></ref><ref id="scirp.56990-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Adm, M. and Garloff, J. (2014) Invariance of Total Nonnegativity of a Tridiagonal Matrix under Element-Wise Perturbation. Operators and Matrices, 8, 129-137. http://dx.doi.org/10.7153/oam-08-06</mixed-citation></ref><ref id="scirp.56990-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Ando, T. (1987) Totally Positive Matrices. Linear Algebra and Its Applications, 90, 165-219.http://dx.doi.org/10.1016/0024-3795(87)90313-2</mixed-citation></ref><ref id="scirp.56990-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Pinkus, A. (2010) Totally Positive Matrices. Cambridge Tracts in Mathematics	 (No. 181). Cambridge University Press, Cambridge.</mixed-citation></ref><ref id="scirp.56990-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Fallat, S.M. and Johnson, C.R. (2011) Totally Nonnegative Matrices. Princeton University Press, Princeton, Oxford. http://dx.doi.org/10.1515/9781400839018</mixed-citation></ref><ref id="scirp.56990-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Adam, M. and Garloff, J. (2013) Interval of Totally Nonnegative Matrices. Linear Algebra and Its Applications, 439, 3796-3806. http://dx.doi.org/10.1016/j.laa.2013.10.021</mixed-citation></ref></ref-list></back></article>