arxiv: v1 [math.ca] 17 Apr 2016

Transcription

1 On some Algebraic Properties for -Meixner Multiple Orthogonal Polynomials of the First Kind arxiv: v1 [math.ca] 17 Apr 2016 Jorge Arvesú and A.M. Ramírez-Aberasturis Department of Mathematics, Universidad Carlos III de Madrid, Avenida de la Universidad, 30, 28911, Leganés, Spain July 3, 2018 Abstract We study a new family of -Meixner multiple orthogonal polynomials of the first kind. The discrete orthogonality conditions are considered over a non-uniform lattice with respect to different - analogues of Pascal distributions. We address some algebraic properties, namely raising and lowering operators as well as Rodrigues-type. Based on the explicit expressions for the raising and lowering operator a high-order linear -difference euation with polynomial coefficients for the -Meixner multiple orthogonal polynomials of the first kind is obtained. Finally, we obtain the nearest neighbor recurrence relation based on a purely algebraic approach. Keywords: multiple orthogonal polynomials; Hermite Padé approximation; difference euations; classical orthogonal polynomials of a discrete variable; Meixner polynomials; -polynomials 2010 Mathematics Subject Classification: 42C05; 33E30; 33C47; 33C65 1 Introduction We begin by introducing the basic background materials. Let µ = (µ 1,...,µ r) be a vector of r positive discrete measures (with finite moments) N i µ i = ω i,k δ xi,k, ω i,k > 0, x i,k R, N i N {+ }, i = 1,2,...,r, k=0 where δ xi,k denotes the Dirac delta function and x i1,k x i2,k, k = 0,...,N i, whenever i 1 i 2. By = (n 1,...,n r) N r we denote a multi-index, where N stands for the set of all nonnegative integers. A type II discrete multiple orthogonal polynomial P (x), corresponding to the multi-index, is a polynomial of degree = n 1 + +n r which satisfies the orthogonality conditions [3] N i k=0 P (x i,k )x j i,k ω i,k = 0, j = 0,...,n i, i = 1,...,r. (1.1) The orthogonality conditions (1.1) give a linear system of homogeneous euations for the + 1 unknown coefficients of P (x). This polynomial solution P always exists. We restrict our attention to a uniue solution (up to a multiplicative factor) with degp (x) =. If this happen for every multi-index, we saythatisnormal [11]. Iftheabovesystemofmeasures forms anat system[3,11] then every multi-index is normal. In this paper we will deal with such system of discrete measures. In [3] some type II discrete multiple orthogonal polynomials on the linear lattice x(s) = s were considered. Moreover, multiple Meixner polynomials (of first and second kind, respectively) were studied. Indeed, the monic multiple Meixner polynomials of the first kind [3] (x), with multiindex N r and degree, satisfy the following orthogonality conditions with different positive parameters α 1,...,α r (0,1) (indexed by α = (α 1,...,α r)) and β > 0 x=0 (x)( x)jυα i,β (x) = 0, j = 0,...,n i, i = 1,...,r, The research of J. Arvesú was partially supported by the research grant MTM C03-01 of the Ministerio de Educación y Ciencia of Spain and mobility research grant from Fundación Caja Madrid. 1

2 where Γ(β +x) α x i υ αi,β, if x R\(Z { β, β,β 2,...}), (x) = Γ(β) Γ(x+1) 0, otherwise, and (x) j = (x)(x+1) (x+j ), (x) 0 = 1, j 1, denotes the Pochhammer symbol. Notice that the multi-index N r is normal whenever 0 < α i < 1, i = 1,2,...,r, and with all the α i different (see [3]). Furthermore, it was found the following raising operators where R α i,β [ M α,β (x)] = + e i (x), i = 1,...,r, (1.2) R α i,β α i(β ) = (1 α i)υ α i,β (x) υα i,β (x), and f(x) = f(x) f(x) denotes the backward difference operator. As a conseuence of (1.2) the following Rodrigues-type formula ( ) nj ( ) (x) = (β) αj Γ(β)Γ(x+1) Γ(β + +x) M α, (1.3) α j Γ(β +x) Γ(β + )Γ(x+1) where M α = ( α x i n i α x i), was obtained. Moreover, two important algebraic properties were deduced for multiple Meixner polynomials of the first kind [3], namely the (r + 1)-order linear difference euation [9] R α i,β (x)+ n i R α j,β where f(x) = f(x+1) f(x), and the recurrence relation [3, 8] ( ( ) ) x αk n i (x) = M α,β + e k (x)+ (β + ) + 1 α k 1 α i + (x) = 0, (1.4) M αβ (x) α in i(β + ) (α i ) 2 e i (x), (1.5) where the multi-index e i is the standard r dimensional unit vector with the i-th entry euals 1 and 0 otherwise. The multiple Meixner polynomials of the first kind M α,β(x) are common eigenfunctions of the above two linear difference operators of order (r + 1), given by (1.4) and (1.5), respectively. In this paper we will introduce a -analogue of such multiple orthogonal polynomials, i.e. when the component measures of µ are different -Poisson distributions and study the aforementioned algebraic properties. Our goal is to continue the recent investigations in [2, 5, 12] for some families of -multiple orthogonal polynomials regarding their algebraic properties. In [4, 9, 13] an (r + 1)-order difference euation for some discrete multiple orthogonal polynomials was obtained. Furthermore, the explicit expressions for the coefficients of (r + 2)-term recurrence relations are a very important issue for the study of some type of asymptotic behaviors for discrete multiple orthogonal polynomials. In [1] the weak asymptotics was studied for multiple Meixner polynomials of the first and second kind, respectively. The zero distribution of multiple Meixner polynomials was also studied. Another interesting fact involving the knowledge of the (r + 2)-term recurrence relations is the attainment of a Christoffel Darboux kernel [6] among other applications, which plays important role in correlation kernel as in the unitary random matrix model with external source. The content of this paper is as follows. In Section 2 we define the -Meixner multiple orthogonal polynomials of the first kind. Moreover, we will prove that these multiple orthogonal polynomials can be explicitly expressed by means of Rodrigues-type formula. This fact provides the background materials for the next Section 3, in which we obtain the (r + 1)-order -difference euation, with polynomial coefficients on a non-uniform lattice x(s). Finally, in Section 4 the nearest neighbor recurrence relation is obtained (an (r + 2)-term recurrence relations). Explicit expressions for the recurrence coefficients are given. The paper ends by summarising our findings in Section 5. 2

3 2 -Meixner multiple orthogonal polynomials of the first kind Aimed to define a new family of -multiple orthogonal polynomials [2] let us consider the following r positive discrete measures on R +, µ i = ω i(k)δ(k s), ω i > 0, i = 1,2,...,r. (2.1) Here ω i(s) = υ α i,β (s) x(s/2), x(s) = ( s )/( ), and α s i Γ (β +s) υ α, if s R + {0}, i,β (s) = Γ (s+1) Γ (β) 0, otherwise, where 0 < α i < 1, β > 0, i = 1,2,...,r, and with all the α i different. Recall that the -Gamma function is given by (1 k+1 ) f(s;) = (1 ) 1 s k 0 Γ (s) = (1 s+k ), 0 < < 1, k 0 (s)(s 2) 2 f ( s; ), > 1. See also [7, 10] for the above definition of the -Gamma function. Lemma 2.1. [5] The system of functions α s 1,x(s)α s 1,...,x(s) n 1 α s 1,...,α s r,x(s)α s r,...,x(s) nr α s r, with α i > 0, i = 1,2,...,r, and (α i/α j) k, k Z, i,j = 1,...,r, i j, forms a Chebyshev system on R + for every = (n 1,...,n r) N r. As a conseuence of Lemma 2.1 the system of measures µ 1,µ 2,...,µ r given in (2.1) forms an AT system on R +. For this system of measures we will define a new family of -multiple orthogonal polynomials. Recall that the -analogue of the Stirling polynomials denoted by[s] (k), is a polynomial of degree k in the variable x(s) = ( s )/( ), i.e. k [s] (k) = j=0 s j = x(s)x(s) x(s k+1) for k > 0, and [s] (0) = 1. Moreover, [s] (k) = (k 2) x k (s)+lower terms = O( ks ), where O( ) stands for the big-o notation. Observethat, when goes to1, thesymbol [s] (k) converges to() k ( s) k, where (s) k is the Pochhammer symbol. Definition 2.2. A polynomial, (s), with multi-index Nr and degree that verifies the orthogonality conditions, (s)[s](k) υ α i,β (s) x(s/2) = 0, 0 k n i, i = 1,...,r, (2.2) (with respect to the measures (2.1)) is said to be the -Meixner multiple orthogonal polynomial of the first kind. In the seuel, we will use the following difference operators def = n i = }{{}, where def = n i times x(s/2), (2.3) x(s+1/2), and x 1(s) def = x(s+1/2) = x(s/2) = s/2. Moreover, we will only consider monic -Meixner multiple orthogonal polynomials of the first kind. 3

4 Proposition 2.3. There holds the following -analogue of Rodrigues-type formula ( ) Γ (β)γ (s+1) Γ(β, (s) = K, α,β M, α + +s), (2.4) Γ (β +s) Γ (β + )Γ (s+1) where M α, = M α i,n i, M α i,n i = (α i) s n i (α i n i ) s, (2.5) and K, α,β = () [ β] ( ) 2 n i n i i+β+j (α i +β+j ) α n i i n i j=i n j, where i = n 1 + +n i, 1 = 0. Proof. We first start by finding the raising operators. Thus, we substitute [s] (k) following finite-difference expression in (2.2) by the i.e., [s] (k) = (k)/2 [k +1] (1) [s+1] (k+1),, (s) [s+1](k+1) υ α i (s) x 1(s) = 0, 0 k n i, i = 1,...,r. Then, using summation by parts along with conditions υ α i,β () = υ α i,β ( ) = 0, we get r raising operators D α i,β, (s) = /2 M α i,1/,β,+ e i (s), i = 1,...,r, (2.6) where D α i,β = ( α ix(β ) (1 α i +β )υ α i/,β (s) υαi,β (s) ), α i,1/ = (α 1,...,α i/,...,α r). Indeed, /2 D α i,β 1 [ ] f(s) = α i β (x(s) x(1 β)) x(s) f(s) 1 α i +β α i +βx(s) f(s), for any function f(s) defined on the discrete variable s. Observe that D α i,β raises by 1 the i-th component of the multi-index in (2.6). Finally, using the raising operators (2.6) in a recursive way one obtains the Rodrigues-type formula (2.4). Remark 2.4. Since we are dealing with an AT-system of positive discrete measures (2.1), then the -Meixner multiple orthogonal polynomial of the first kind, (s) has exactly different zeros on R + (see [3, Theorem 2.1, pp ]). 4

5 3 High-order -dif ference euation The strategy that we will follow to deduce the (r + 1)-order -difference euation is the foliowing. First step. Define an r-dimensional subspace V of polynomials on the variable x(s) of degree at most 1 by means of interpolatory conditions. (s) as a linear combination of the basis elements of V. Third step. Combine the lowering and the raising operators (2.6) to get an(r+1)-order -difference euation (in the same fashion that [4], [5], and [9]). These steps represent the general features of the algebraic approach we are using in this section, however some ad hoc computations are needed because of the dependence of the explicit expressions involved in the above steps on the given family of multiple orthogonal polynomials. Second step. Find the lowering operator and express its action on, Lemma 3.1. Let V be the linear subspace of polynomials Q(s) on the lattice x(s) of degree at most 1 defined by the following conditions Q(s)[s] (k) υ α j,β+1 (s) x 1(s) = 0, 0 k n j 2 and j = 1,...,r. Then, the system {M α i,,β+1, e i (s)} r, where α i, = (α 1,...,α i,...,α r), is a basis for V. Proof. Considering the orthogonality relations M α i,,β+1, e i (s)[s] (k) υ α j,β+1 (s) x 1(s) = 0, 0 k n j 2, j = 1,...,r, we have that polynomials M α i,,β+1, e i (s), i = 1,...,r, belong to V. Let proceed by reductio ad absurdum and assume that there exists constants λ i, i = 1,...,r, such that λ im α i,,β+1, e i (s) = 0, where Then, multiplying the previous euation by [s] (n k) on s from 0 to, one gets Hence, from relations λ i M α i,,β+1, e i M α i,,β+1, e i (s)[s] (n k) (s)[s] (n k) υ α k,β λ i > 0. (s) x 1(s) and then taking summation υ α k,β (s) x 1(s) = 0. υ α k,β (s) x 1(s) = cδ i,k, c R\{0}, (3.1) we deduce that λ k = 0 for k = 1,...,r. Here δ i,k represents the Kronecker delta symbol. Therefore, {M α i,,β+1, e i (s)} r is linearly independent in V. Furthermore, we know that any polynomial of V can be determined with coefficients while ( r) linear conditions are imposed on V, conseuently the dimension of V is at most r. Hence, the system {M α i,,β+1, e i (s)} r spans V, which completes the proof. Now we will prove that operator (2.3) is indeed a lowering operator for the seuence of -Meixner multiple orthogonal polynomials of the first kind +1, (s). Lemma 3.2. There holds the following relation, (s) = Proof. Using summation by parts we have, (s)[s](k) υ α j,β+1 n i+1/2 1 α i n i+β 1 α i +β[ni](1) M α i,,β+1, e i (s). (3.2) (s) x 1(s) = =, (s) [ [s] (k) υ α j,β+1 (s) ] x 1(s), (s)ϕ j,k(s)υ α j,β (s) x 1(s), (3.3) 5

6 where ϕ j,k (s) = 1/2 ( β x(s) x(β) ) +1 [s] (k) /2 x(s) α jx(β) [s](k), is a polynomial of degree k + 1 in the variable x(s). Conseuently, from the orthogonality conditions (2.2) we get, (s)[s](k) υ α j,β+1 (s) x 1(s) = 0, 0 k n j 2, j = 1,...,r. Hence, from Lemma 3.1,, combination of polynomials {M α i,,β+1, e i (s) V. Moreover, M α,β, (s) can univocally be expressed as a linear (s)} r, i.e., (s) = β im α i,,β+1, e i (s), Multiplying both sides of the euation (3.4) by [s] (n k) one has, (s)[s](n k) υ α k,β+1 (s) x 1(s) = β i = β k υ α k,β+1 M α i,,β+1, e k β i > 0. (3.4) (s) x 1(s) and using relations (3.1) M α i,,β+1, e i (s)[s] (n k) υ α k,β+1 (s)[s] (n k) υ α k,β+1 (s) x 1(s) (s) x 1(s). (3.5) If we replace [s] (k) by [s] (n k) in the left-hand side of euation (3.3), then left-hand side of euation (3.5) transforms into relation, (s)[s](n k) υ α k,β+1 (s) x 1(s) = = /2( 1 α k n k+β ) α k x(β), (s)ϕ k,n k (s)υ α k,β (s) x 1(s), (s)[s](n k) υ α k,β (s) x 1(s). (3.6) Here we have used that x(s)[s] (n k) = [s] (n k) to get ϕ k,nk (s) = /2 (1 α k n k +β ) α k [s] (n k) x(β) + lower terms. On the other hand, from (2.6) one has that ( ) /2 1 α k +β υ α k,β (s), α k x(β) (s) = /2 [ υ α k,β+1 (s)m α i,,β+1, e k (s) ]. (3.7) Then, by conveniently substituting (3.7) in the right-hand side of euation (3.6) and using once more summation by parts, we get, (s)[s](n k) υ α k,β+1 (s) x 1(s) = 1 α k n k+β 1 α k +β [s] (n k) [ = 1 α k n k+β 1 α k +β υ α k,β+1 Since [s] (n k) = 3/2 n k [n k ] (1) [s] (n k) we finally have, (s)[s](n k) = n k+1/21 α k n k+β 1 α k +β[n k] (1) (s)m α i,,β+1, e k ] (s) x 1(s) M α i,,β+1, e k (s) [ [s] (n ] k) υ α k,β+1 (s) x 1(s). υ α k,β+1 (s) x 1(s) M α i,,β+1, e k (s)[s] (n k) υ α k,β+1 (s) x 1(s). 6

7 Therefore, comparing this euation with (3.5) we obtain the coefficients in the expansion (3.4), i.e. which proves relation (3.2). β k = n k+1/21 α k n k+β 1 α k +β[n k] (1), Theorem 3.3. The -Meixner multiple orthogonal polynomial of the first kind, (s) satisfies the following (r + 1)-order -difference euation D α i,β, (s) = n i+11 α i n i+β 1 α i +β[ni](1) Proof. Since operators (2.6) are commuting, we write D α i,β = D α j,β Dα i,β D α j,β, (s). (3.8), (3.9) and then using (2.6), by acting on euation (3.2) with the product of operators (3.9), we obtain the following relation D α i,β, (s) = n i+1/2 1 α i n i+β 1 α i +β[ni](1) D α j,β Dα i,β M α i,,β+1, e i (s) which proves (3.8). = n i+11 α i n i+β 1 α i +β[ni](1) 4 Nearest neighbor recurrence relation D α j,β, (s), Let us start recalling that for any function f(s) defined on the discrete variable s and a positive integer n i there holds (see Lemma 5.1 from [5]) M α i,n i (x(s)f(s)) = n i+1/2 x(n i)(α i) s n i (α i n i ) s f(s)+ x(s) x(ni) where difference operator M α i,n i is given in (2.5). Now, let us proceed with the nearest neighbor recurrence relation. n i D α i,n i f(s), (4.1) Theorem 4.1. The -Meixner multiple orthogonal polynomials of the first kind satisfy the following (r + 2)-term recurrence relation where x(s) α,β M, (s) = c α,β,k M,+ e k (s)+b M α,β,k α i α j n j + x(n i) α i n i αj n j αi +n i x(β + ) α i +β+n i α,β M, (s) = ( ) M K, α,β α,β, (s) and b,k = α i +β α i +β+n i ( ) r x(n i), (s) α i +β α i n i α i +β+n i α i +β+n i 1 α i +β+n i 2 i x(n αi n j j=i i) α i +β +( ) +β i M α,β, e i (s), (4.2) 1 (α k nk+1 )x(β + ) α j +β + 1 α k +β+n k+1 α i n i α i +β ) 7

8 Proof. Consider euation α i α j n j x(n ( i) α i n i ) α i β+ +ni α i n i αj n j α i +β+n i α i +β+n i, c,k = /2 α i +β (α k nk+1 )x(β + ) α i +β+n i α k +β+n k+1. (α k ) s nk+1 (α k nk+1 ) s Γ (β + +1+s) Γ (β + +1)Γ (s+1) [ ( )] = (α k ) s n k s+1/2 (α k nk+1 ) s Γ (β + +1+s) Γ (β + +1)Γ (s+1) = 1/2 (α k ) s n k (α k n k ) s Γ (β + +s) Γ (β + )Γ (s+1) + 1/2 α k +β+n k+1 (α k n k+1 )x(β + ) (α k) s n k (α k n k ) s Γ (β + +s) x(s) Γ (β + )Γ, (s+1) which can be rewritten in terms of difference operators (2.5) as follows /2 M α k,n k +1 Γ (β + +1+s) Γ (β + +1)Γ (s+1) = M α k,n k ( Γ(β + +s) Γ (β + )Γ (s+1) ( 1+ α )) k +β+n k+1 (α k n k+1 )x(β + ) x(s). (4.3) Since operators (2.5) are commuting the multiplication of euation (4.3) from the left-hand side by the product M α j yields j k,n j α k +β+n k+1 Γ (β + +s) (α k n k+1 )x(β + ) M α,x(s) Γ (β + )Γ (s+1) = /2 M,+ e α Γ (β + +1+s) Γ (β + +s) k Γ (β + +1)Γ (s+1) M α, Γ (β + )Γ. (4.4) (s+1) Using recursively relation (4.1) involving the product of r difference operators M α 1,n 1,...,M αr,n r acting on the function f(s) = Γ(β+ +s) Γ (β+ )Γ (s+1), we have M,x(s) α Γ (β + +s) Γ (β + )Γ (s+1) = 1/2 α i α j n j x(n ( i) α i n i )( α j +β+n j ) α i n i αj n j ( αν +β ) ν=1 ( M α Γ l (β + +s),n l δ l,i Γ l=1 (β + )Γ + α i +β+n i (s+1) α i +β x(s) i x(n i) αi n j j=i α i +β ( ) +β α i n i α i +β +( ) r i 1 x(n i) α j +β M α, Γ (β + +s) Γ (β + )Γ. (4.5) (s+1) Hence, by using expressions (4.4), (4.5) one gets Γ (β + +s) Γ (β + )Γ (s+1) = /2 α i +β (α k nk+1 )x(β + ) Γ (β + +1+s) α i +β+n i α k +β+n k+1 M α,+ e k Γ (β + +1)Γ (s+1) x(s)m α, 8

9 + α i +β α i +β+n i ( ) r x(n i) M α, i x(n αi n j j=i i) α i +β +( ) +β i Γ (β + +s) Γ (β + )Γ (s+1) /2 M α l,n l δ l,i l=1 1 (α k nk+1 )x(β + ) α j +β + 1 α k +β+n k+1 Γ (β + +s) Γ (β + )Γ. (s+1) As a result of the above calculations one has x(s)m α, ) α i α j n j x(n ( i) α i n i ) α i n i αj n j α i +β+n i α i n i α i +β Γ (β + +s) Γ (β + )Γ (s+1) = /2 α i +β (α k nk+1 )x(β + ) Γ (β + +1+s) α i +β+n i α k +β+n k+1 M α,+ e k Γ (β + +1)Γ (s+1) +b,k M, α Γ (β + +s) α i α j n j x(n ( i) α i n i ) Γ (β + )Γ (s+1) /2 α i n i αj n j α i +β+n i 1 α i +β+n i M α l,n l δ l,i l=1 Γ (β + +s) Γ (β + )Γ. (s+1) Finally, multiplying from the left both sides of the previous expression by Γ(β)Γ(s+1) Γ (β+s) Rodrigues-type formula (2.4) we obtain (4.2), which completes the proof. and using 5 Concluding remarks In closing, we summarize our findings. We have defined a new family of -Meixner multiple orthogonal polynomials of the first kind and obtained their explicit expression in terms of Rodrigues-type formula (2.4). We have shown that these multiple orthogonal polynomials are common eigenfunctions of two different (r + 1)-order difference operators given in (3.8) and (4.2). By taking limit 1 one recovers the corresponding algebraic relations for multiple Meixner polynomials of the first kind [3]. The expressions (2.4), (3.8), and (4.2) transform into (1.3), (1.4), and (1.5), respectively. Our algebraic approach for the nearest neighbor recurrence relation (4.2) is purely algebraic and it neither reuire to introduce type I multiple orthogonality [14] nor an algebraic Riemann-Hilbert approach. Indeed, the -difference operators involved in the Rodrigues-type formula are the base of the discussed approach. Acknowledgements The research of J. Arvesú was partially supported by the research grant MTM C03-01 (Ministerio de Economía y Competitividad) of Spain. References [1] Aptekarev A., Arvesú J., Asymptotics for multiple Meixner polynomials, J. Math. Anal. Appl. 411 (2014), [2] Arvesú J., On some properties of -Hahn multiple orthogonal polynomials, J. Comput. Appl. Math. 233 (2010), [3] Arvesú J., Coussement J., Van Assche W., Some discrete multiple orthogonal polynomials, J. Comput. Appl. Math. 153 (2003), [4] Arvesú J., Esposito C., A high-order -difference euation for -Hahn multiple orthogonal polynomials, J. Difference Eu. Appl. 18 (2012), [5] Arvesú J., Ramírez-Aberasturis A.M., On the -Charlier multiple orthogonal polynomials, SIGMA 11 (2015), 026. [6] Daems E., Kuijlaars A.B.J., A Christoffel Darboux formula for multiple orthogonal polynomials, J. Approx. Theory 130 (2004),

10 [7] Gasper G., Rahman M., Basic hypergeometric series, Encyclopedia of Mathematics and its Applications, Vol. 96, 2nd ed., Cambridge University Press, Cambridge, [8] Haneczok M., Van Assche W., Interlacing properties of zeros of multiple orthogonal polynomials,j. Math. Anal. Appl. 389 (2012), [9] Lee D.W., Difference euations for discrete classical multiple orthogonal polynomials, J. Approx. Theory 150 (2008), [10] Nikiforov A.F., Suslov S.K., Uvarov V.B., Classical orthogonal polynomials of a discrete variable, Springer Series in Computational Physics, Springer-Verlag, Berlin, [11] Nikishin E.M., Sorokin V.N., Rational approximations and orthogonality, Translations of Mathematical Monographs, Vol. 92, Amer. Math. Soc., Providence, RI, [12] Postelmans K., Van Assche W., Multiple little -Jacobi polynomials, J. Comput. Appl. Math. 178(2005), ,. [13] Van Assche W., Difference euations for multiple Charlier and Meixner polynomials, in Proceedings of the Sixth International Conference on Difference Euations, CRC, Boca Raton, FL, 2004, [14] Van Assche W., Nearest neighbor recurrence relations for multiple orthogonal polynomials, J. Approx. Theory 163 (2011),