arxiv: v3 [math.ca] 9 Mar 2011

Transcription

1 Structured matrices, continued fractions, and root localization of polynomials arxiv:094703v3 [mathca] 9 Mar 0 Olga Holtz Department of Mathematics, University of California-Berkeley Institut für Mathematik, Technische Universität Berlin School of Mathematics, Institute for Advanced Study Mikhail Tyaglov Institut für Mathematik, Technische Universität Berlin December 3, 009; last revision February 9, 0 Abstract We give a detailed account of various connections between several classes of objects: Hankel, Hurwitz, Toeplitz, Vandermonde and other structured matrices, Stietjes and Jacobi-type continued fractions, Cauchy indices, moment problems, total positivity, and root localization of univariate polynomials Along with a survey of many classical facts, we provide a number of new results Contents Complex rational functions and related topics 4 Hankel and Hurwitz matrices Hurwitz formulæ 4 Resultant formulæ and applications: discriminant and Orlando formulæ 8 3 Euclidean algorithm, the greatest common divisor, and continued fractions: general case 4 4 Euclidean algorithm: regular case Finite continued fractions of Jacobi type 5 Euclidean algorithm: doubly regular case Finite continued fractions of Stieltjes type 7 Real rational functions and related topics 35 Sturm algorithm and Frobenius rule of signs for Hankel minors 36 Cauchy indices and their properties 37 3 Rational functions mapping the upper half-plane to the lower half-plane 44 3 General theory 44 3 Some classes of infinite Hankel matrices and the finite moment problem on the real axis 5 33 R-functions as ratios of polynomials R-functions with a Stieltjes continued fraction expansion 59 4 The number of distinct real zeros of polynomials Polynomials with all real zeros 66 4 The number of distinct positive, negative and non-real zeros of polynomials Stieltjes continued fractions of the logarithmic derivative 66 4 Polynomials with real zeros 70 The work of OH was supported by the Sofja Kovalevskaja Research Prize of Alexander von Humboldt Foundation and by the National Science Foundation under agreement DMS holtz@mathiasedu The work of MT was supported by the Sofja Kovalevskaja Research Prize of Alexander von Humboldt Foundation tyaglov@mathtu-berlinde

2 Introduction This survey and research paper offers a glimpse at several classical topics going back to Descartes, Gauss, Stieltjes, Hermite, Hurwitz and Sylvester (see, eg, [9, 80, 8, 8, 8]), all connected by the idea that behavior of polynomials can be analyzed via algebraic constructs involving their coefficients Thus a number of inter-related algebraic constructs was built, including Hurwitz, Toeplitz and Hankel matrices, the corresponding quadratic forms, and the corresponding continued fractions The linear-algebraic properties of these objects were shown to be intimately related to root localization of polynomials such as stability, whereby the zeros of a polynomial avoid a specific half-plane, or hyperbolicity, whereby the zeros lie on a specific line These methods gave rise to several well-known tests of root localization, such as the Routh-Hurwitz algorithm or the Lienard-Chipart test In the 0th century, this line of research was developed further by Schur [4], Pólya [78], Krein [59, 60, 34] and others (see, eg, [43]), leading to important notions of total positivity, Pólya frequency sequences, stability preservers etc A very important part of that research effort was devoted to entire functions, in particular, the Laguerre-Pólya and related classes However, this area gradually went out of fashion and was essentially abandoned around s, with a few exceptions, such as Padé approximation We should note in this connection that the closely related theory of continued fractions, initiated by Chebyshev, Stieltjes and Markov and further developed by Akhiezer, Krein and their collaborators in connection with problems of mechanics is only now returning to the forefront due to its connections with orthogonal polynomials [58, 66, 69, 49, 50] On the other hand, various fast algorithms for structured matrices (Cauchy, Vandermonde, Toeplitz, Hankel matrices, and more generally quasiseparable matrices and matrices with displacement structure) have been developing rapidly in the last few decades due to the efforts of Gohberg, Olshevsky, Eidelman, Kailath, Heinig, Sakhnovich, Lerer, Rodman, Fuhrmann and others, along with applications to control theory and engineeting: see the collections [70, 7, 7, 53, 5, ] and the references therein Surveying these developments is, however, outside the scope of this paper Furthermore, we must stress the fact that the results and algorithms in this paper are not approximate but exact (ie, give exact answers in exact arithmetic) as is the term localization, which is often used in the sequel For instance, we perform root localization with respect to a line or a half plane and are able to determine exactly, eg, how many zeros of a given polynomial lie on a given line For sure, this does not obviate the need for error analysis of the corresponding algorithms that use finite precision or floating-point arithmetic We must note that the questions of root localization are returning into the spotlight also due to the development of basic algebraic techniques for multivariate polynomials [4, 40, 8, 7, 6, 83] and entire functions [9, 0], as well as due to newly found applications of results on multivariate stable and hyperbolic polynomials to other branches of mathematics [4, 4, 00, 6] Very promising applications of current interest include problems of discrete probability, combinatorics, and statistical physics, such as the analysis of partition functions arising in classical Ising and Potts models of statistical mechanics [9, 89] The goal of this paper is to provide a comprehensive and coherent treatment of classical connections between three kinds of objects: various structured matrices, representations of rational functions via continued fractions, and root localization of univariate polynomials Our additional goal is to demonstrate that, despite the rich history of this subject, even the univariate case is far from being exhausted, and that classical algebraic techniques are useful in answering questions about the behavior of polynomials and rational functions Our future goals include using this work for generalizations of these classical results to univariate entire and meromorphic functions, including questions on Pólya frequency sequences, Hurwitz rational functions, and entire functions with all real zeros [98, ] One particular area of interest is the class of so-called generalized Hurwitz polynomials, which is a useful large class containing all Hurwitz stable polynomials This line of our ongoing research is closely related to the theory of indefinite metric spaces, indefinite moment problems, and the corresponding eigenvalue problems [3, 76, 77, 47] We now illustrate our main point above by discussing several new results in this paper Our first illustration is provided by the body of work in Section 4 These results provide explicit criteria for a polynomial to have only real roots in terms of Hankel and Hurwitz determinants made of its coefficients To give the reader an idea of these statements, we quote two sample theorems from Section 4:

3 3 Consider a polynomial and let p(z) = a 0 z n +a z n + +a n, a,,a n R, a 0 > 0 L(z):= p (z) p(z) = s 0 z + s z + s z 3 + be its logarithmic derivative Define the associated infinite Hankel matrix S:=S(L):= s i+j i,j=0 and its leading principal minors s 0 s s s j s s s 3 s j D j (S):=D j (S(L)):=, j =,,3, s j s j s j+ s j Here is a sample result from Section 4 that provides an interesting connection between total positivity of a special Hurwitz matrix and polynomials all whose roots are real and negative: Theorem (Total positivity criterion for negative zeros) The polynomial g(z) = a 0 +a z + +a n z n, a 0 0, a n > 0 has all negative zeros if and only if the infinite matrix a 0 a a a 3 a 4 a 5 a 6 0 a a 3a 3 4a 4 5a 5 6a 6 0 a 0 a a a 3 a 4 a 5 D (g):= 0 0 a a 3a 3 4a 4 5a a 0 a a a 3 a a a 3a 3 4a 4 is totally nonnegative This result is based on a new connection between Stieltjes continued fractions representing a rational function P = q/p and a special factorization of the infinite Hurwitz matrix H(p,q) associated to the pair (p, q) (see Section 34) This connection implies a criterion of total nonnegativity of the latter infinite Hurwitz matrix H(p, q) (Theorem 344), which in turn implies our criterion for negative zeros of a polynomial g via the total nonnegativity of the infinite discriminant matrix D (g) One direction of Theorem 344 was essentially developed in the earlier works of Asner and Kemperman[6, 57], while the other direction is given here for the first time This result simultaneously provides a criterion of stability [6, 57, 46, 55, 98] Section 4 provides many other alternative criteria for real zeros, including the cases of only positive, only negative, or mixed zeros, in terms of Hankel and Hurwitz determinants as well as the coefficients of Stieltjes continued fraction expansion of the logarithmic derivative of a given polynomial The special role played by the logarithmic derivative for the analysis of real zeros is clarified earlier in Section 4, together with the main ideas behind counting real roots to the left and to the right of the origin The underlying theory of Cauchy indices is presented in Section 33 Another important class of objects behind the results of this paper on polynomial roots is the class of so-called R-functions These are rational functions mapping the upper-half plane of the complex plane either to the lower half-plane or to itself(see Definition 3) The theory of R-functions connects several key objects of this work: continued fractions, Cauchy indices, polynomials with real roots, the moment problem, and determinantal inequalities As an example, we quote the following new result from Section 33: Theorem (Generalized Lienard-Chipart criterion) If a real rational function R = q/p, with p and q relatively prime, maps the upper half-plane of the complex plane to the lower half-plane, then the number of

4 4 its positive poles is equal to the number V (a 0,,a n ) of strong sign changes in the sequence (a 0,,a n ) of coefficients of its denominator p In particular, R has only negative poles if and only if a j > 0 for j =,,,n, and only positive poles if and only if a j a j < 0 for j =,,,n Our final illustration is the generalized Orlando s formula from Section Orlando s formula per se expresses the Hurwitz determinant of order n associated to a polynomial p of degree n as the product, up to a certain normalization, of all possible pairs z i +z j of the zeros of p Orlando s result goes back to 9 [73] The generalized Orlando formula obtained in this paper says the following: Theorem (Generalized Orlando formula) Let p be a polynomial of degree n and q a polynomial of degree m n Then the resultant of these polynomials can be computed as follows: R(p,q) = ( ) n(n ) a m+n (z i +z j ) 0 i<j n where a 0 is the leading coefficient of p and where z i, i =,,n, are the zeros of the polynomial h(z):=p(z )+zq(z ) The classical Orlando s formula then follows as a special case by splitting an arbitrary polynomial into its even and odd part and applying the generalized Orlando s formula As these examples show, new connections can be found between several classical matrix classes (Hurwitz, Hankel, Vandermonde etc) made of the coefficients of polynomials, different representations of rational functions (Stieltjes and Jacobi continued fractions, Laurent series, elementary fractions), and various counting notions (sign changes, the number of roots in a specific domain of a complex plain, Cauchy indices) Underlying all these is a coherent theory demonstrating how these connections arise The presented theory is completely general in that it does not single out the special case of stable polynomials The importance of stability, historically the first question on root localization [6, 5, 3, 84, 37], is by now quite well understood Various stability criteria have been studied in great detail, especially in the engineering literature In this survey (Section 4), we choose instead to illustrate applications of the general theory to polynomials with real roots, a less studied but arguably equally important special case To summarize, the point of this work is to provide a uniform, streamlined, algebraic treatment for a body of questions centered around the root localization of polynomials At present, results of this type are scattered in the literature, whereas even textbooks and monographs about polynomials sometimes fail to provide answers to some basic questions (ie, tests for real roots using only the coefficients of the polynomial) We hope that this work will serve as a useful reference and source of further research for mathematicians in various fields interested in polynomials and their zeros Complex rational functions and related topics Consider a rational function z R(z):= q(z) p(z) where p and q are polynomials with complex coefficients () p(z) := a 0 z n +a z n + +a n, a 0, a,, a n C, a 0 0, () q(z) := b 0 z n +b z n + +b n, b 0, b,, b n C, (3) If the greatest common divisor of p and q has degree l, 0 l n, then the rational function R has exactly r = n l poles Note that zeros or poles of a rational function are always counted with multiplicities unless explicitly stipulated otherwise Thus, in the rest of the paper, the phrase counting multiplicities will be implicit in every statement about zeros or poles of functions under consideration Hankel and Hurwitz matrices Hurwitz formulæ In this section, we introduce Hankel and Hurwitz matrices associated to a given rational functions and discuss the Hurwitz formulæ that connect these two classes

5 5 Expand the function () into its Laurent series at : Here s j = 0 for j < n m and s n m = b 0 a 0, where m = degq The sequence of coefficients of negative powers of z R(z) = s + s 0 z + s z + s + (4) z3 s 0,s,s, defines the infinite Hankel matrix S:=S(R):= s i+j i,j=0 This gives the correspondence R S(R) (5) Definition For a given infinite sequence (s j ) j=0, consider the determinants s 0 s s s j s s s 3 s j D j (S):=, j =,,3,, (6) s j s j s j+ s j ie, the leading principal minors of the infinite Hankel matrix S These determinants are referred to as the Hankel minors or Hankel determinants An infinite matrix is said to have finite rank r if all its minors of order greater than r are zero whereas there exists at least one nonzero minor of order r Knonecker [63] proved that, for any infinite Hankel matrix, any minor of order r where r is the rank of the matrix, is a multiple of its leading principal minor of order r This implies the following result Theorem (Kronecker [63]) An infinite Hankel matrix S = s i+j i,j= has finite rank r if and only if D r (S) 0, (7) D j (S) = 0, for all j > r (8) Let D j (S) denote the following determinants s s s 3 s j s s 3 s 4 s j+ D j (S):=, j =,,3,, (9) s j s j+ s j+ s j With a slightabuse ofnotation, wewill alsowrite D j (R):=D j (S(R)) and D j (R):= D j (S(R)) if the matrix S = S(R) is made of the coefficients (4) of the function R The following theorem was also established by Gantmacher in [36] Theorem 3 An infinite matrix S = s i+j i,j= has finite rank if and only if the sum of the series R(z) = s 0 z + s z + s z 3 + is a rational function of z In this case the rank of the matrix S is equal to the number of poles of the function R Theorems and 3 have a simple corollary, which will be useful later Corollary 4 A rational function R with exactly r poles represented by the series (4) has a pole at the point 0 if and only if D r (R) 0 and Dj (R) = 0 for j = r,r+, (0)

6 6 Proof Since the function R represented as a series (4) has exactly r poles, the function G(z):=zR(z) zs = s 0 + s z + s z + s 3 z 3 + has exactly r poles if and only if R has a pole at the point 0 If R does not have a pole at 0, then G has r poles Thus, by Theorems and 3, the function R has a pole at 0 if and only if D r (G) 0 and D j (G) = 0, j = r,r +, () Since D j (R) = D j (G), j =,,, the formula () yields the assertion of the corollary Theorems and 3 imply the following: if the greatest common divisor of the polynomials p and q defined in () (3) has degree l, 0 l degq, then the formulæ (7) (8) hold for r = n l for the rational function (), where n is the degree of the polynomial (), ie, the denominator of R Denote by j (p,q) the following determinants of order j: a 0 a a a j a j a j b 0 b b b j b j b j 0 a 0 a a j a j a j j (p,q):= 0 b 0 b b j b j b j, j =,,, () a 0 a a j b 0 b b j constructed using the coefficients of the polynomials () (3) Here we set a i = 0 for all i > n and b l = 0 for l > m The determinants j (p,q) are called determinants of Hurwitz type or just Hurwitz determinants or Hurwitz minors In his celebrated work [48], A Hurwitz found relationships between the minors D j (R) defined in (6) and the determinants j (p,q) defined by () Theorem 5 ([48, 6,36, ]) Let R(z) = q(z), where the polynomials p and q are defined by () (3) p(z) Then j (p,q) = a j 0 D j(r), j =,, (3) Proof From (4) it follows that coefficients s k,a k,b k satisfy recurrence relations b j = a 0 s j +a s j + +a j s, j = 0,,, (4) These recurrence relations imply the formulæ (3) by direct matrix multiplication, once we take into

7 7 account (4) (cf [6]): a j+ 0 D j (R) = ( ) j(j ) a j s s 0 s j s j s j = 0 0 s s j 3 s j s j s 0 s s j s s 0 s j a s s 0 s j s j s j = ( ) j a a s s j 3 s j s j 3 = a j a j a s s 0 s j a 0 a a a j a j a j 0 b 0 b b j b j b j 0 a 0 a a j a j a j = ( ) j 0 0 b 0 b j 3 b j b j = a 0 j (p,q) b 0 b b j a 0 a a j Remark 6 Since the formulæ (3) are algebraic identities also valid for polynomials p = n j=0 a jz j, a 0 0, and q = n j=0 b jz j, they also hold when p and q are entire functions and R is a meromorphic function From () one can see that, whenever b 0 0, we have j (p,q) = j ( q,p), j =,, (5) This observation, coupled with Theorem 5, yields the following well-known fact (cf [7] ) Corollary 7 For two infinite sequences S:=(s j ) j=, T :=(t j) j= with s 0 and t 0, the following conditions are equivalent: ) s t + = 0, s i t +s i t 0 + +s 0 t i +s t i = 0, i = 0,,, ) D k (S) = s k D k (T), k =,,, where D k (S) and D k (T) are the Hankel minors (6) for the sequences S and T, respectively Proof Let then R(z) = q(z) p(z) = s + s 0 z + s z + s z 3 +, R(z) = p(z) q(z) = t + t 0 z + t z + t z 3 + It appears that this basic fact was known much earlier than the work of Edrei but we do not know of the original reference

8 8 The assertion of the corollary follows from Theorem 5, formula (5) and the identity ( s + s 0 z + s z + s + )( z 3 t + t 0 z + t z + t ) z 3 + Resultant formulæ and applications: discriminant and Orlando formulæ This section is devoted to the resultant, an important object that is closely related to the greatest common divisor of two polynomials The discriminant of a single polynomial occurs as a special case of the resultant of that polynomial and its derivative Moreover, the generalized Orlando formula connects the resultant of two polynomials, p and q, with the roots of the aggregate polynomial h(z):=p(z )+zq(z ) Although the resultant is classically known, this last connection is new Definition 8 Let two polynomials p and q be given in () (3) with n:=degp and m:=degq The resultant of p and q is the following determinant of order n+m a 0 a a m a m a n a n a n+m 0 a 0 a m a m a n a n a n+m R(p,q):= 0 0 a 0 a a n m a n m+ a n b n m b n m+ b n b n b n m b n m b n, (6) 0 b n m b n b n b n m b n m b n b n m b n m+ b n where we set a i :=0 and b i :=0 for all i > n From (6) it can be immediately seen that R(p,q) = ( ) nm R(q,p) The resultant of two polynomials is a multi-affine function of the roots of these polynomials, as the following formula shows Theorem 9 Given polynomials p and q as in () (3) with b 0 0, let λ i (i =,,n) denote the zeros of p, and let µ j (j =,,n) denote the zeros of q Then R(p,q) = ( ) n(n ) n (p,q) = a n 0 n q(λ i ) = a n 0 bn 0 i= n (λ i µ j ) = ( ) n b n 0 i,j= n p(µ j ) (7) Proof At first, assume that all roots of the polynomial p are distinct Then the function R admits the representation n α j R(z) = β +, (8) z λ j where j= α j = p(λ j) q, j =,,n (9) (λ j ) Here λ j λ i whenever i j, according to the assumption Comparing the representation (8) with the expansion (4), we obtain n s k = α i λ k i, k = 0,,, (0) i= j=

9 9 and s = b 0 a 0 0 From (), (3) and (6) it follows that R(p,q) = ( ) n(n ) n (p,q) = ( ) n(n ) a n 0 D n(r) () The formulæ (0) yield s 0 s s s n s s s 3 s n = s n s n s n+ s n α α α 3 α n λ λ λ n α λ α λ α 3 λ 3 α n λ n λ λ = α λ α λ α 3 λ 3 α n λ λ n n λ 3 λ 3 λ n 3 α λ n α λ n α 3 λ3 n α n λn n λ n λ n λn n This equality implies s 0 s s s n s s s 3 s n D n (R) = = s n s n s n+ s n n α i i= λ λ λ 3 λ n λ λ λ 3 λ n λ n λ n λ n 3 λn n () This formula combined with the residue formula (9) implies Since we have n p (λ i ) = a n 0 i= D n (R) = n n q(λ i ) p (λ i λ j ) (3) (λ i ) i= p (z) = a 0 n i=k= k i i= j<i n (z λ k ), k= k i n (λ i λ k ) = a n 0 i λ j j<i(λ ) (λ i λ j ) i<j This product has exactly n(n ) factors of the form λ i λ j, therefore Now from (3) (4) we obtain n i= p (λ i ) = a n 0 ( )n(n ) (λ i λ j ) (4) j<i ( )n(n ) D n (R) = a n 0 n i= q(λ i ) = ( ) n(n ) bn 0 a n 0 n (λ i µ j ) (5) i,j= The formula (7) follows from () and (5) If the polynomial p has multiple zeros, we can consider an approximating polynomial p ε with simple zeros such that lim ε 0 p ε(z) = p(z) for all z

10 0 Then n (R ε ) ε 0 n (R), where R ε (z):= q(z) p ε (z) The formula (7) is valid for the polynomials q and p ε, so it is also valid at the limit, ie, for the polynomials q and p Corollary 0 Under the conditions of Theorem 9, let degq = m n Then = a m 0 bn n m i=j= R(p,q) = ( ) n(n ) a m n 0 n (p,q) = a m 0 n m m (λ i µ j ) = ( ) nm b n n m p(µ j ) j= n q(λ i ) = i= = ( ) nm R(q,p) (6) Proof If degq = m( n), then the first nonzero coefficient of q is b n m, according to (3) Therefore, n (p,q) = ( ) n(n ) Thus, we have a 0 a a n m a n m a n a n m a n 0 a 0 a n m a n m a n a n m 3 a n a 0 a m a n a n+m a 0 a n m a n b n m b n b n m b n b n b n m b n b n m b n n (p,q) = ( ) n(n ) a n m 0 R(p,q) This relation and the formulæ (3), (5) yield (6) when p has only simple roots But (6) is also valid when p has multiple zeros, which can be proved by an approximation argument as above The formulæ (7) and (6) now imply the well-known property of the resultant: Corollary R(p,q) = 0 if and only if the polynomials p and q have common roots Next, we consider a function that allows us to test whether a single polynomial has multiple roots Definition Given a polynomial () with roots λ i (i =,,n), its discriminant is defined as D(p):=a n 0 n (λ i λ j ) (7) It is clear from (7) that the discriminant of a polynomial is equal to zero if and only if the polynomial has multiple zeros But multiple zeros of a polynomial are the zeros that it shares with its derivative The following connection between the discriminant of p and the resultant of p and p should not come as a surprise Theorem 3 Given a polynomial () of degree n, we have j<i R(p,p ) = ( ) n(n ) a 0 D(p) (8)

11 Proof Indeed, the resultant R(p,p ) satisfies (6) Together with (4) and (7), it gives R(p,p ) = a n 0 n i= p (λ i ) = a n 0 ( ) n(n ) (λ i λ j ) = ( ) n(n ) a 0 D(p) j<i From (6) we obtain n (p,p ) = ( ) n(n ) a 0 R(p,p ) = a 0D(p) Thus, the discriminant of the polynomial p is the following (n ) (n ) determinant: na 0 (n )a (n )a a n 0 0 a 0 a a a n a n 0 0 na 0 (n )a a n a n 0 D(p) = 0 a 0 a a n a n na 0 3a n 3 a n 0 a a 0 a n 3 a n a a a n na 0 (n )a a n (9) We now give one more formula for the resultant, which is very close to the well-known Orlando formula [73] (see also [36]) More precisely, the application of this formula to the resultant of the odd and the even parts of a polynomial yields exactly the Orlando formula Theorem 4 Let the polynomials p and q be given by () (3), with degp = n and degq = m n Then the resultant of these polynomials can be computed as follows: R(p,q) = ( ) n(n ) a m+n 0 i<k n where z i (i =,,n) are the zeros of the following polynomial of degree n (z i +z k ), (30) h(z) = p(z )+zq(z ) (3) Proof Let λ i (i =,,n) be the zeros of the polynomial p, and let µ j (j =,,m) be the zeros of the polynomial q From (3) it follows that and p(z ) = h(z)+h( z), q(z ) = h(z) h( z) (3) z p(µ j ) = h(± µ j ), j =,,m, q(λ i ) = ± h(± λ i ) λi, i =,,n (33) From (3) we obtain p(zk ) = h( z k), q(zk ) = h( z k), k =,,n, (34) z k where z k are the zeros of the polynomial h Since degq = m by assumption, we conclude that b n m 0 but b 0 = = b n m = 0 Thus, (6) implies R(p,q) = b n n m m p(µ j ) (35) j=

12 Using (35), (33) and (34), we have (R(p,q)) Thus, we obtain = b n = a m 0 n m j= n = a m+n 0 m h( µ j )h( µ j ) = b n q(z k ) = am 0 k= n n k= i= i k n k= h( z k ) z k (z i +z k ) = a m+n 0 R(p,q) = ±a m+n 0 m n n m a 0 j= k= n a 0 = a m 0 k= i<k n i<k n (µ j z k) = a m z k n i= (z i +z k ) (z i +z k ) 0 (z i +z k ) To determine the sign, we consider the special case h(z) = (z ) n Then k= p(z ) = (z )n +(z +) n k= e πik n n k=, q(z ) = (z )n (z +) n, z k= m b n m (zk µ j ) so that degp = n, degq = n The roots of the odd part q can be determined from the equation (z ) n = (z +) n except that the zero root should be discarded This shows that the roots of q are {wk : k =,,n } where w k is defined from the equation /(w k + ) = e πik n, k =,,n Consequently, ( ) n n n n p(wk ) = n = ie πik n sin ( ) πk n n = i n(n ) k= e πik sin n( πk n ) = ( ) n i= k n k= sin n( πk n Thus, according to the last formula in (6), the sign of the resultant R(p,q) in our special case is ( ) n(n ) signb n ( ) n(n ) = ( ) n(n+) = ( ) n+n(n ) As we already established, the resultant of p and q is a polynomial in the coefficients of p and q, hence in the coefficients of h The expression ( ) n(n ) a n 0 i<k (z i+z k ) is a symmetric function of the roots of h multiplied by its leading coefficient to the power n, and hence also a polynomial in the coefficients of h Since the two polynomials must be identically equal, we conclude that the sign ( ) n(n ) occurs at all times whenever degp = n, degq = n We now show how to produce a formula for the case m < n from the already established formula for m = n Given a polynomial q of degree n, set b through b n m to zero, thus obtaining a polynomial of degree m Observe what happens to the determinantal expression (6) Since the lower left block of size n (n m ) is now filled with zeros, the upper-triangular submatrix above produces the factorization R(p,q) = a0 n m R(p,q), where R(p,q) denotes the old resultant of p and q constructed as if degq were equal to (n ) Thus the new resultant R(p, q) satisfies the equation a0 n m R(p,q) = ( ) n(n ) a n 0 (z i +z j ), hence R(p,q) = ( ) n(n ) a n+m 0 i<k n (z i +z j ) Our next statement can be proved analogously i<k n ) j=

13 3 Theorem 5 Let the polynomials p and q be given by () (3), and let degq = m n = degp Then the resultant of these polynomials can be computed by the formula R(p,q) = ( ) n(n ) a m+n where z i (i =,,n+) are the zeros of the polynomial 0 i<k n+ g(z) = q(z )+zp(z ) (z i +z k ), (36) The famous Orlando formula is a simple consequence of Theorems 4 5 Before proving the Orlando formula, we introduce the following determinants for the polynomial () a a 3 a 5 a 7 a j a 0 a a 4 a 6 a j 0 a a 3 a 5 a j 3 j (p) = 0 a 0 a a 4 a j 4, j =,,n, (37) a j where we set a i :=0 for i > n Definition 6 The determinants j (p) (j =,,n) are called the Hurwitz determinants or the Hurwitz minors of the polynomial p where It is easy to see that the polynomial p can be always represented as follows p(z) = p 0 (z )+zp (z ), (38) p 0 (u) = a u l +a 3 u l + +a n, p (u) = a 0 u l +a u l + +a n, (39) if the degree n of the polynomial p(z) is odd: n = l+, and p 0 (u) = a 0 u l +a u l + +a n, p (u) = a u l +a 3 u l + +a n, (40) if n = l Theorem 7 (Orlando, [73, 36]) Let the polynomial p of degree n be given by () Then the determinant n (p) defined by (37) can be computed from the formula n (p) = ( ) n(n ) a n where z i, i =,,n, are the zeros of the polynomial p This equality is known as the Orlando formula 0 i<j n (z i +z j ), (4) Proof At first, let the degree n of p be odd, n = l + Then (38) (39) show that degp 0 l, degp = l, and the leading coefficient of p is a 0 Thus, (), (6) and (36) imply n (p) = ( ) l l (p,p 0 ) = ( ) l(l+) a l degp0 0 R(p,p 0 ) = ( ) l a l 0 which coincides with (4) since n = l+ i<j l+ (z i +z j ),

14 4 If n = l, then (38) and (40) imply degp l, degp 0 = l, and the leading coefficient of the polynomial p 0 is a 0 As above, the formulæ (), (6) and (30) can now be combined to obtain n (p) = a 0 l(p 0,p ) = ( ) l(l ) a l degp 0 R(p 0,p ) = ( ) l a l which is exactly the formula (4) 0 i<j l (z i +z j ), 3 Euclidean algorithm, the greatest common divisor, and continued fractions: general case In the previous subsection, we considered the resultant of two polynomials and observed that it is equal to zero if and only if these polynomials have a nontrivial common divisor The standard way to find their greatest common divisor is via the Euclidean algorithm Let us consider polynomials p and q given by () (3) and let us denote f 0 (z):=p(z), f (z):=q(z) b 0 a 0 p(z) Construct a sequence of polynomials f 0,f,,f k (k n) by the following formula f j (z) = q j (z)f j (z)+f j+ (z), j =,,k (f k+ (z) = 0), (4) where q j is the quotient and f j+ is the remainder from the division of f j by f j The last polynomial in this sequence, f k, is the greatest common divisor of the polynomials f 0 and f (and also all other polynomials f j in the sequence) In other words, where h k (z) = Now, denote Rewriting (4), we obtain f j (z) = h j (z)f k (z), j = 0,,,k, (43) R j (z) = R j (z):= f j(z) f j (z) = h j(z) h j (z) (44), j =,,k, (45) q j (z)+r j+ (z) where R k+ (z) 0 Using this equality, we can represent the function R (z) as a continued fraction 3 : R (z) = f (z) f 0 (z) = h (z) h 0 (z) = q (z)+ q (z)+ q 3 (z)+ + q k (z) (46) It is easy to see that, for each j = 0,,k, the polynomial h j is the leading principal minor of order k j of the following k k tridiagonal matrix q k (z) q k (z) q k (z) 0 0 J(z) = (47) q (z) q (z) Thus, degf < degf 0 Obviously, if degq < degp, that is, if b 0 = 0, then f = q 3 The functions R = q/p and R are related via the formula R(z) = R (z)+b 0 /a 0

15 5 In particular, h 0 (z) = detj(z) Thus, the determinant of the matrix J(z) is the denominator of R(z) This is a very useful observation, as certain properties of the function R turn out to be connected to the location of the eigenvalues of the generalized (matrix polynomial) eigenvalue problem J(z)u = 0 (48) Conversely, the behavior of the eigenvalues of the problem (48) can provide information about certain properties of the function R Later in this paper we will give examples of such interrelation More generally, in addition to the fraction (46) we may consider the fractions R j (z) = h j(z) h j (z) = q j (z)+ q j+ (z)+ q j+ (z)+ + q k (z), j =,,k The continued fraction expansion (46) can be found efficiently in terms of the corresponding Hankel minors, as the next theorem shows Theorem 8 Two rational functions R and G vanishing at are related via the identity R(z) = where g is a polynomial of degree m( ) with leading coefficient α if and only if and g(z)+g(z), (49) g(z) = αz m + (50) D (R) = D (R) = = D m (R) = 0, (when m > ) (5) D m+j (R) = ( ) m(m ) ( )j D j (G) α m+j, j = 0,,,, (5) where D 0 (G):=, and the determinants D j (R) and D j (G) are defined by (6) Proof Since the functions R and G vanish at, they can be expanded into Laurent series R(z) = s 0 z + s z + + s m z m + s m z m+ +, G(z) = t 0 z + t z + t + (53) z3 Moreover, the conditions (49) (50) hold if and only if In fact, if R satisfies (49) (50), then and s 0 = s = = s m = 0 and s m 0 (54) s i = lim z zi+ R(z) = lim z s m = lim z zm R(z) = lim z = 0, i = 0,,,,m (55) g(z) G(z) + zi+ z i+ α+ γ z m + +γ m z m Now assume that the condition (54) holds and that R(z) = + G(z) z m f(z)+g(z) = 0 (56) α for some polynomial f If degf = j for some j < m, that is 4, if f(z) = γz j +, where γ 0, then s j = lim z zj R(z) = γ 0, 4 The degree j of f cannot be 0, since the limit lim z R(z) = γ is nonzero for a constant nonzero function f(z) = γ

16 6 contrary to (54) On the other hand, if the degree of f is greater than m, then s m = lim z zm R(z) = 0 Thus, f must be of exact degree m since otherwise a contradiction arises with one of the conditions (54) Also note that the equalities (54) are equivalent to (5) Moreover, from (54) we have s m s m s m D m (R) = 0 s m s m s m 4 s m 3 s m s m s m+ s m 3 s m For convenience, denote the coefficients of the polynomial g as follows g(z)=:αz m +t m z m +t m+ z m + +t z +t = ( ) m(m ) s m m 0 (57) If the functions R and G satisfy (49) (50), then, according to (53) (54), we have [ sm z m + s m z m+ + s m+ + ][ αz m +t zm+ m z m + +t + t 0 z + t z + t ] z 3 + This identity implies the following relations: s m = α, s m+j = j i=0 t m+i α s m+j i, j = 0,,, (58) Now, the equality (5) for D m (R) follows from (57) (58) For a fixed number j, consider the determinant D m+j (R): s m s m s m+j 0 0 s m s m s m+ s m+j D m+j (R) = 0 s m s m 4 s m 3 s m s m+j 3 s m s m s m 3 s m s m s m+j (59) s m s m+ s m s m s m s m+j s m+j s m+j s m+j 3 s m+j s m+j s m+j Addtoeachi th column(i = m+j,m+j,,3,)columnsi,i,,multipliedby t m α, t m+ α,, t j 3 α, t j, respectively This eliminates the entries in the upper right corner of the determinant (59) α using (58) So, we can rewrite the original determinant as a product of the following two determinants of order m and j, respectively: s m d d d,j d,j 0 0 s m 0 d d d,j d,j D m+j (R) = d 3 d 3 d 3,j d 3,j, 0 s m 0 0 s m d j d j d j,j d jj where d i,i := s m+i ti α s m+i 3 ti α s m ti +i 3 s m ti +i (60) α α From (58) and (60) we obtain s m t 0 t t j t j D m+j (R) = ( ) m(m ) s m s m 0 0 t t t j t j ( )j α m+j s m+j 3 s m+j 4 s m 0 t j t j t j 4 t j 3 s m+j s m+j 3 s m s m t j t j t j 3 t j

17 7 Conversely, if the equalities (5) (5) hold, then we obtain s 0 = s = = s m = 0 from (5) by induction, which was already proved to be equivalent to the fact that R satisfies (49), with a polynomial g of degree at least m and D m (R) = ( ) m(m ) s m m If degg > m, then s m = lim z zm R(z) = lim z g(z) z m + G(z) = 0, z m therefore, D m (R) = ( ) m(m ) s m m = ( )m(m ) α m = 0 and = 0, according to (5) Thus, α from (5) (5) we get D i (R) = 0 for all i N, which is impossible, since R is a rational function, hence at least one of the minors D i (R) must be nonzero due to Theorem Consequently, the polynomial g satisfies (50) with α 0 If a rational function R with exactly r poles is expanded into a continued fraction (46), then the rational functions R j defined in (44) satisfy the relations (45), where q j (z) = α j z nj +, α j 0, j =,,,k (6) Here n +n + +n k = r and n, since degf < degf 0, as is remarked above The other degrees n i are greater or equal to due to the structure of the Euclidean algorithm (4) From (5) (5) we obtain, for a fixed integer j (j =,,,k), the following product formulæ: D n+n + +n j (R ) = ( ) n(n ) ( ) = ( ) n(n ) ( ) n j + n i i=3 α This chain of equalities results in the formula D n+n + +n j (R) = j i= ( ) ni(ni ) ( ) j i=n i j i=n i n j + n i i= α n j + n i i= α D n3+n 4+ +n j (R 3 ) = j in i+ i=0 j i= α n i+ j n ρ ρ=i+ i D n+n 3+ +n j (R ) = ( ) n(n ) ( ) j i=3n i, j =,,,k (6) Remark 9 Our discussion above can be summarized as follows: Suppose that a rational function R with r poles has a continued fraction expansion (46) with polynomials q j satisfying (6) Then Theorems, 3 and 8 and formulæ (6) imply that all Hankel minors D j (R) are equal to zero, except for the minors 5 D n (R),D n+n (R),,D n+n + +n k (R), which are not zero and which can be calculated from the formulæ (6) Using Theorem 8 and the formulæ (6), we now describe equivalence classes of rational functions whose sequences of Hankel minors are the same Theorem 0 Two rational functions 5 Recall that n +n + +n k = r R(z) = s 0 z + s z + s z 3 + and G(z) = t 0 z + t z + t z 3 +,

18 8 both vanishing at infinity, have equal Hankel minors if and only if their continued fraction expansions D j (R) = D j (G), j =,,, (63) R(z) = q (z)+ q (z)+ q 3 (z)+ + q k (z) and G(z) = q (z)+ q (z)+ q 3 (z)+ + q k (z) satisfy k = k =:k, (64) and the polynomials q j and q j, for each j (j =,,,k), have equal degrees and equal leading coefficients: q j (z) = α j z nj +, q j (z) = α j z nj +, j =,,,k (65) Proof If the continued fraction expansions of the functions R and G satisfy the conditions (64) (65), then the equalities (63) follow from Remark 9 and the formulæ (6) Conversely, let the Hankel minors associated with the functions R and G satisfy (63) Therefore, by Theorems and 3, the functions R and G have equal number of poles Moreover, they can be then represented as follows: R(z) = R (z) = q (z)+r (z) and G(z) = G (z) = q (z)+g (z) It follows from Theorem 8 that the degrees and the leading coefficients of the polynomials q and q coincide According to (5), the functions R and G have equal Hankel minors So we can apply the same argument to them Thus, Theorem 8 shows that the degrees and the leading coefficients of each pair of the polynomials q j, q j are equal to each other and, consequently, the number of those polynomials must be the same as well Since the functions R and G have equal number of poles as was proved above, the equality (64) thus follows Remark Note that the equalities D j (R) = D j (G) do not imply the equality of the functions R G Counterexamples are quite easy to construct For instance, we can take R(z) = G(z) = Then R G but D (R) = D (R) = D (G) = D (G) = z = z + z + z z 4 +, z z = z + z + z z 4 + z 3 Finally, let the function R have a continued fraction expansion(46) and also a Laurent series expansion where n = degq R(z) = s n + s n z n z n+ + s n + +, (66) zn+

19 9 Consider the functions F j (z):= Q j(z) P j (z) = q (z)+ q (z)+ q 3 (z)+ + q j (z), j =,,k, (67) constructed using the polynomials (6) by the Euclidean algorithm (4) Definition The polynomials Q j are called partial numerators, the polynomials P j partial denominators, and the fractions F j partial quotients of R The denominator P j (z) of the fraction F j (z) is the jth leading principal minor of the matrix q (z) q (z) q 3 (z) 0 0, q k (z) q k (z) and P k (z) = h 0 (z) (see (46)) Let m j denote the sum of the degrees n through n j (see (6)): m j = n +n + +n j, j =,,,k (68) Then degp j = m j Notice that, for a fixed number j ( j k), the initial terms of the Laurent series of the function F j coincide with those in the Laurent series (66) of R up to, and including, the term s m j z mj : F j (z) = Q j(z) P j (z) = s n + s n z n z n+ + + s m j z mj + s (j) z m j mj+ + (69) In fact, each coefficient s i of the series (66) can be found from the recurrence relation [ s i = lim z i+ R(z) s z n z i n+ s n z i+n s i z ], i = n,n,n +, (70) Using the expansions (46) and (67) of the functions R and F j, respectively, together with the formula (70), we obtain (69) To find an explicit formula for the polynomials P j that depends only on the coefficients s i, we introduce the following notation: P j (z)=:p 0,j z mj +P,j z mj + +P mj,jz +P mj,j, j =,,,k For a fixed number j between and k, the formula (69) implies Q j (z) = P j (z) [ s n + s n z n z n+ + + s m j z mj + s (j) z m j + mj+ We must require that the coefficients of z,z,,z mj be zero, which leads to the system ] (7) s 0 P mj,j +s P mj,j + +s mj P,j +s mj P 0,j = 0, s P mj,j +s P mj,j + +s mj P,j +s mj+p 0,j = 0, s mj P mj,j +s mj P mj,j + +s mj P,j +s mj P 0,j = 0, (7)

20 0 where we set s 0 :=s := :=s n :=0 By Cramer s rule, the solution to the system (7) satisfies s 0 s s i s i+ s mj P mj i,j = ( )mj i P 0,j s s s i s i+ s mj+ i = 0,,,,m D mj (R), j, j =,,,k s mj s mj s mj+i s mj+i s mj This formula implies the following representation for the polynomials P j : s 0 s s s mj P j (z) = P s s s 3 s mj+ 0,j D mj (R), j =,,,k (73) s mj s mj s mj+ s mj z z z mj Here, according to the notation (6), P 0,j = j α i, j =,,,k This formula, combined with (6) and the notation (68), yields i= P 0,j P 0,j+ j ( )ni (n i ) = ( ) j D mi (R) D i= mi (R) k for j =,,, j ( )ni (n i ) = ( ) j D mi (R) D mi (R) i= ( )nj+(n j+ ) n i n i ( )ni(n i ) D mi (R) D mi (R) ( )ni(n i ) D mi (R) D mi (R) n i+ D mj (R) k for j = 0,,,, D mj+ (R) n i n i In case the Euclidean algorithm is regular, the notion that will be introduced in the sequel, these formulæ take a simpler form The coefficients of Q j can be obtained from the coefficients of P j via the formula (7): denoting Q j (z) = Q 0,j z mj n +Q,j z mj n + +Q mj n,jz +Q mj n,j, we find the coefficients of Q j by the following 0 0 s 0 s s n+i P s 0 s s mj n i+ s mj n i+ s mj Q i,j = ( )mj 0,j s s s mj n i+ s mj n i+3 s mj+ D mj (R) s mj s mj s mj n i s mj n i+ s mj for i = 0,,,m j n Note that the rational functions F j (z) = Q j(z) P j (z) the function R, see [0, 43] are exactly the diagonal elements of the Padé table of

21 4 Euclidean algorithm: regular case Finite continued fractions of Jacobi type In this section we discuss the best known case of the Euclidean algorithm, which leads to orthogonal polynomials, 3-term recurrence relations and other related phenomena We give conditions for regularity and discuss the form of partial quotients and the generalized eigenvalue problem corresponding to the regular case Suppose that all the polynomials q j resulting from an application the Euclidean algorithm (4) are linear: q j (z) = α j z +β j, α j,β j C, α j 0, j =,,r (74) We call this situation the regular case of the Euclidean algorithm In the regular case, degf j (z) = n j, j = 0,,r, (75) and r = n l, where l is the degree of the greatest common divisor f r of the polynomials f 0 and f (and of all other polynomials f j in the sequence) Thus, the polynomials f j satisfy the following three-terms recurrence relation: f j (z) = (α j z +β j )f j (z)+f j+ (z), j =,,r (76) Consequently, the function R expands into the continued fraction R(z) = f (z) f 0 (z) = h (z) h 0 (z) = α z +β + where the polynomials h 0 and h are defined by (43) α z +β + α 3 z +β α r z +β r, (77) Definition 3 Continued fractions of the type (77) are usually called J-fractions or continued fractions of Jacobi type Remark 4 If a rational function G satisfies the condition lim G(z) = c, 0 < c <, then we will z say that G(z) has a J-fraction expansion if the function R(z) = G(z) lim G(z) can be represented as z in (77) Theorem 8 with some simple modifications implies the following corollary, which will be useful later Corollary 5 Two rational functions R and G, both vanishing at, satisfy the following condition if and only if where D 0 (G):= R(z) = αz +β +G(z), α 0, D j (R) = ( )j α j D j (G), j =,,, (78) Using Corollary 5, we can now prove a criterion when a rational function expands into a J-fraction Theorem 6 ([74]) A rational function with exactly r poles has a J-fraction expansion if and only if where D j (R) are the Hankel determinants defined in (6) R(z) = s + s 0 z + s z + s + (79) z3 D j (R) 0, j =,,r, (80)

22 Proof In view of Remark 4 it is sufficient to consider the case s = 0 Suppose that the function R has a J-fraction expansion (77), ie, a continued fraction expansion (46) with polynomials q j satisfying (74) with k = r Then the formulæ (6) yield D j (R) = ( ) j(j ) j i= α j i+ i 0, j =,,,r (8) Now suppose that the inequalities (80) hold Then D (R ) = s () 0 0, where we denote 6 Therefore, R (z) can be represented as follows R (z):=r(z) = s() 0 z + s() z + s() z 3 + R (z) = α z +β +R (z), α = s () 0 0, Now, Corollary 5 implies that D j (R ) = ( )j D j (R ), j =,, Therefore, the function R is α j of the same type as R, that is, R satisfies the conditions 7 D j (R ) 0, j =,,,r, which are analogous to (80) In particular, s () 0 0, where If we continue this process, then for each function R (z) = s() 0 z + s() z + s() z 3 + R j (z) = s(j) 0 z + s(j) z + s(j) z 3 + we obtain that s (j) 0 0, j =,,,r Consequently, R j (z) = α j z +β j +R j+ (z), j =,,,r, where R r+ 0 This means that the function R can be represented as a J-fraction (see (45) (46)) This theorem and Theorem 5 imply the following statement, which was proved in [0] (Theorem 4) with inessential differences Corollary 7 Let R(z) = q(z), where p and q are defined in () (3), and let the function R have p(z) exactly r poles The function R can be expanded into a J-fraction (77) if and only if where j (p,q) are defined in () j (p,q) 0, j =,,r, (8) Corollary 8 Given a rational function R(z) = f (z) f 0 (z), where degf < degf 0, the Euclidean algorithm (4) applied to the polynomials f 0 and f is regular if and only if the inequalities (80) hold 6 Recall that we assumed s = 0 7 For j r, we have D j+ (R ) = D j (R ) = 0, and therefore R has exactly r poles

23 3 If we make an equivalence transformation of the fraction (77), as in [0, p 66], and set then the J-fraction (77) takes the form d 0 := α, d j := α j α j+, j =,,,r (83) R(z) = e j := β j α j, j =,,,r, (84) z e d 0 z e Moreover, from (8) and (83) we have (see [0, p 67]) d d d r z e r (85) In fact, j D j (R) = D j (R) d i (86) i= D j (R) = ( ) j(j ) j i= α j i+ i = ( ) j(j ) j α j i= α i j α j i i = ( ) j(j ) ( ) (j )(j ) ( ) j j D j (R) = D j (R) α α i α i+ The expression (86) implies an interesting formula: where we set D (R):=D 0 (R):= i= i= = j d i d j = D j (R)D j+ (R) D j (R), j = 0,,,r, (87) Remark 9 From the formulæ (83), (84) and (87) it follows that both forms of J-fraction expansions (77) and (85) of the function R(z) are unique, that is, all coefficients in (77) and (85) are defined uniquely for the function R Remark 30 (cf [0, p70]) Suppose that the function R has a J-fraction expansion (77) Then R( z) = α z β + α z β + α 3 z β α r z β r i= (88) According to Remark 9, (a properly normalized) J-fraction expansion is unique Therefore, if the function R is odd, ie, R(z) R( z), then all β j in (88) must be equal to zero Conversely, if β j are all equal to zero, then R(z) is obviously an odd function of z Thus, the formula (86) implies D j (R) = j i=0 d j i i (89) In other words, if the function R has a J-fraction expansion (85), the minors D j (R) do not depend on the coefficients e j, j =,,r, whereas the minors D j (R) obviously do In order to establish the dependence of these minors on the coefficients of the fraction (85), we first prove the following simple fact

24 4 Lemma 3 Let the complex rational function with exactly r poles, has a J-fraction expansion (85) Then R(z) = s 0 z + s z + s + (90) z3 D r (R) r D r (R) = k= µ k, (9) where µ k are the poles of the function R, and the minors D j (R) and D j (R), j =,,, are defined by (6) and (9) Proof If R has a pole at zero, then the lemma holds true, since in this case D r (R) 0 by Theorem 6, while D r (R) = 0 by Corollary 4 Now assume that R has no a pole at zero At first, suppose that the function R has only simple poles µ k µ j whenever k j, so it can be represented as a sum of partial fractions r ν k R(z) = z µ k This formula together with (90) gives the following well-known formulæ On the other hand, () implies s j = D r (G) = r k= r ν k k= Further, recall that D r (R) = D r (Φ), where Φ(z) = zr(z) = So, analogously to (93), we obtain r ν k + k= D r (R) = D r (Φ) = Now (9) follows from (93) and (94) k= ν k µ j k, j = 0,,, (9) µ µ µ r µ r µ r r µ r r k= r (ν k µ k ) k= ν k µ k z µ k = s 0 + s z + s z + s 3 z 3 + µ µ µ r µ r µ r r µ r (93) (94) Let R = q and suppose that the function R has multiple poles, or, equivalently, that the polynomial p p has multiple zeros In this case, we can consider an approximating polynomial p ε with simple zeros such that lim p ε(z) = p(z) for all z ε 0 Then D r (R ε ) D r (R), Dr (R ε ) Dr (R), ε 0 ε 0 where R ε (z):= q(z) p ε (z) The formula (9) is valid for the polynomials q and p ε, so it is also valid at the limit, ie, for the polynomials q and p, since the product of all zeros of the polynomial p ε (z) tends to the product of all zeros of the polynomial p whenever ε 0

25 5 Let us consider the following tridiagonal matrix e d d e d d e J r = e r dr dr e r (95) constructed using the coefficients of the J-fraction (85) Suppose that the coefficient d 0 in (85) equals and consider the partial quotients of the J- fraction (85) F j (z) = Q j(z) P j (z) =, j =,,r (96) d z e z e d d j z e j Then F r = R and the polynomial P j is the characteristic polynomial of the leading principal submatrix of J r of order j, for any j =,,r Using the formula (9), it is now easy to prove the following theorem Theorem 3 Let the matrix J n be defined by (95) and let the function R defined by (90) have a J-fraction expansion (85) (with d 0 = ) Then the leading principal minors J r m, m =,,r, of the matrix J r can be found by the following formulæ: J r m = D m (R), m =,,r (97) D m (R) Proof At first, we note that the formula (97) holds for det(j r ) = J r r Indeed, on the one hand, the determinant of the matrix J r is the product of its eigenvalues On the other hand, the eigenvalues of the matrix J r are the poles of the function R, so by (9), their product equals D r (R) D r (R) Let us now turn to the functions F j introduced in (96) According to (69), we have F m (z) = Q m(z) P m (z) = s 0 z + s z + + s m z m + s(m) m +, m =,,r, (98) zm+ where the coefficients s i, i = 0,,,m coincide with those of the function R defined in (90) Consequently, D j (F m ) = D j (R), Dj (F m ) = D j (R), j =,,m Thus, applying Lemma 3 to the function F m and to the leading principal submatrix of J r of order m, we obtain J r m = D j (F m ) D j (F m ) = D m (R) D m (R), m =,,n, as required Now suppose that the function R has a J-fraction expansion (85), where d 0 may differ from Then we can consider the function G(z) := R(z) d 0 It is clear that Consequently, according to (97), we get D m (R) = d m 0 D m(g), Dj (R) = d m 0 D j (G) m =,,r J r m = D m (G) D m (G) = D m (R), m =,,r (99) D m (R)