Multiple orthogonal polynomials associated with an exponential cubic weight
|
|
- Ashlynn Beasley
- 5 years ago
- Views:
Transcription
1 Multiple orthogonal polynomials associated with an exponential cubic weight Walter Van Assche Galina Filipu Lun Zhang Abstract We consider multiple orthogonal polynomials associated with the exponential cubic weight e x over two contours in the complex plane. We study the basic properties of these polynomials, including the Rodrigues formula and nearest-neighbor recurrence relations. It turns out that the recurrence coefficients are related to a discrete Painlevé equation. The asymptotics of the recurrence coefficients, the ratio of the diagonal multiple orthogonal polynomials and the (scaled) zeros of these polynomials are also investigated. Keywords: multiple orthogonal polynomials, exponential cubic weight, Rodrigues formula, nearest-neighbor recurrence relations, string equations, discrete Painlevé equation, zeros, asymptotics Introduction and statement of the results. Orthogonal polynomials associated with an exponential cubic weight A sequence of non-constant monic polynomials {p n } with deg p n n is said to be orthogonal with respect to the exponential cubic weight e x if p n (x)x e x dx = 0, = 0,,..., n, (.) Γ where the contour Γ is chosen such that the above integral converges. These polynomials satisfy the three-term recurrence relation where β n = xp n (x) = p n+ (x) + β n p n (x) + γ np n (x), (.) Γ xp n(x)e x dx Γ p n(x)e x dx, γ n = xp Γ n(x)p n (x)e x dx, (.) Γ p n (x)e x dx Department of Mathematics, KU Leuven, Celestijnenlaan 00B box 400, BE-00 Leuven, Belgium. Walter.VanAssche@wis.uleuven.be Faculty of Mathematics, Informatics and Mechanics, University of Warsaw, Banacha, Warsaw, 0-097, Poland. filipu@mimuw.edu.pl School of Mathematical Sciences and Shanghai Key Laboratory for Contemporary Applied Mathematics, Fudan University, Shanghai 004, People s Republic of China. lunzhang@fudan.edu.cn
2 and the initial condition is taen to be γ0 p = 0. It is shown by A. Magnus [9] that the recurrence coefficients β n and γn satisfy the string equations γn+ + β n + γ n = 0, (.4) γn(β n + β n ) = n. (.5) For the convenience of the reader, we derive the string equations using ladder operators for orthogonal polynomials in the Appendix. Some variants of orthogonal polynomials associated with the exponential cubic weight have recently been studied in the context of numerical analysis [8] and random matrix theory [4]. Γ Γ 0 Γ Figure : The three rays Γ 0, Γ, Γ For our purpose, we are concerned with the polynomials for specific contours Γ. Consider the three rays (see Figure ) Γ = {z C : arg z = ω }, = 0,,, (.6) where ω = e πi/ is the primitive third root of unity and the orientations are all taen from left to right. Clearly, the integral (.) is well-defined for each Γ. We shall denote by p () n the polynomials satisfying (.) with Γ = Γ 0 Γ. The corresponding recurrence coefficients will be accordingly denoted by β n () and (γ n () ). Hence, we have p () n (x)x e x dx = 0, = 0,,..., n, (.7) Γ 0 Γ and From (.), it is readily seen that xp () n (x) = p () n+(x) + β () n p () n (x) + (γ () n ) p () n (x). (.8) β () 0 = Γ 0 Γ xe x dx Γ 0 Γ e x dx = Γ(/) Γ(/) eπi/. (.9) Thus, one can determine (β n (), (γ n () ) ) recursively from the string equations (.4) (.5) with initial condition γ () 0 = 0 and β () 0 = Γ(/) Γ(/) eπi/.
3 In a similar manner, we let p () n be the polynomials satisfying (.) with Γ = Γ 0 Γ, and denote by β n () and (γ n () ) the corresponding recurrence coefficients. To this end, one has β () Γ 0 = 0 Γ xe x dx Γ 0 Γ e x dx = β() 0 = Γ(/) Γ(/) e πi/. (.0) For the recurrence coefficients β (i) n and (γ (i) n ), i =,, the following proposition holds. Proposition.. There exist two real sequences a n and b n, n N = {0,,,,...} such that and a n, b n satisfy the coupled difference relations with initial conditions β () n = b n e πi/, (γ () n ) = a n e πi/, (.) a n + a n+ = b n, (.) a n (b n + b n ) = n, (.) a 0 = 0, b 0 = Γ(/) Γ(/). (.4) Similarly, we have with the same sequences a n and b n. β () n = b n e πi/, (γ () n ) = a n e πi/, (.5) From (.), one can easily eliminate a n in (.) and obtain n b n + b n + n + b n + b n+ = b n. (.6) This difference equation belongs to A c -type equation on the list of discrete Painlevé equations by Grammaticos and Ramani [6, 7], which has a connection with the second Painlevé equation. It is also an alternative discrete Painlevé I equation in Clarson s list [6, Appendix A.4], see also [], [0]. We give a short derivation of the string equations (.4) (.5) in the Appendix, where we also deal with the more general weight e x +tx.. Multiple orthogonal polynomials with an exponential cubic weight Multiple orthogonal polynomials are polynomials of one variable which are defined by orthogonality relations with respect to r different measures µ, µ,..., µ r, where r is a positive integer. As a generalization of orthogonal polynomials, multiple orthogonal polynomials originated from Hermite-Padé approximation in the context of irrationality and transcendence proofs in number theory. They were further developed in approximation theory, we refer to Aptearev et al. [, ], Coussement and Van Assche [8], Niishin and Soroin [, Chapter 4, ], and Ismail [8, Chapter ] for more information.
4 We tae r = and for (, l) N, we are interested in the monic polynomials P,l of degree + l which satisfy the orthogonality conditions x i P,l (x)e x dx = 0, i = 0,,...,, (.7) Γ 0 Γ x i P,l (x)e x dx = 0, i = 0,,..., l. (.8) Γ 0 Γ We call P,l the (type II) multiple orthogonal polynomial for the exponential cubic weight. If one of and l is equal to zero, then P,l reduce to the usual orthogonal polynomials with respect to the exponential cubic weight e x, i.e., P,0 (x) = p () (x), P 0,(x) = p () (x), (.9) where p (i), i =, are defined in Section.. It is the aim of this paper to derive some basic properties of P,l. Our main results are Theorem. (Rodrigues formula). Let n, m N = {0,,,,...}, then e x P n,n+m (x) = ( )n d n ( ) e x P n dx n 0,m (x), (.0) e x P n+m,n (x) = ( )n d n ( ) e x P n dx n m,0 (x). (.) where P 0,m (x) and P m,0 (x) are given in (.9). The polynomials P n,n (x) were already mentioned by Pólya and Szegö in their problem boo [4, Part V, Chapter, Problem 59] and Pólya investigated their zeros in [, Satz IV]. They are also a special case of polynomials introduced by Gould and Hopper [5] and were investigated, among others, by Dominici [9] and Paris []. Their multiple orthogonality (or d-orthogonality, if one only considers the diagonal polynomials) was already noted earlier, see e.g., [] and references there. In this paper we are investigating the full range of polynomials P n,m (x) and not only the diagonal polynomials, but we obtain ratio asymptotics and the distribution of the zeros for the diagonal polynomials in Section 4. For asymptotic approximations and an asymptotic expansion of P n,n (x) we refer to [9] and []. Multiple orthogonal polynomials satisfy a system of nearest-neighbor recurrence relations [8, Theorem.7]. For P,l defined in (.7) (.8) we can represent the recurrence coefficients explicitly in terms of the sequences a n and b n in Proposition., as stated in the following theorem. Theorem. (the nearest-neighbor recurrence relations). Let n, m N, then xp n,n+m (x) = P n+,n+m (x) + c n,n+m P n,n+m (x) xp n,n+m (x) = P n,n+m+ (x) + d n,n+m P n,n+m (x) + a n,n+m P n,n+m (x) + b n,n+m P n,n+m (x), (.) + a n,n+m P n,n+m (x) + b n,n+m P n,n+m (x), (.) 4
5 where c n,n+m = { Γ(/) Γ(/) eπi/, m = 0, b m e πi/, m > 0, (.4) d n,n+m = b m e πi/, (.5) { n a n,n+m = Γ(/) i, m = 0, Γ(/) (.6) nam m eπi/, m > 0, b n,n+m = { n Γ(/) Γ(/) i, m = 0, (n+m)a m m e πi/, m > 0. (.7) Similarly, where xp n+m,n (x) = P n+m+,n (x) + c n+m,n P n+m,n (x) xp n+m,n (x) = P n+m,n+ (x) + d n+m,n P n+m,n (x) + a n+m,n P n+m,n (x) + b n+m,n P n+m,n (x), (.8) + a n+m,n P n+m,n (x) + b n+m,n P n+m,n (x), (.9) c n+m,n = b m e πi/, (.0) { Γ(/) d n+m,n = Γ(/) e πi/, m = 0, (.) b m e πi/, m > 0, { n a n+m,n = Γ(/) i, m = 0, Γ(/) (n+m)a m (.) e πi/, m > 0, m { n b n+m,n = Γ(/) i, m = 0, Γ(/) (.) nam m e πi/, m > 0. Here, a n and b n are the two real sequences generated from (.) (.4). It is also easy to chec that the recurrence coefficients derived in Theorem. satisfy the partial difference equations obtained in [7, Theorem.]. The rest of this paper is organized as follows. Theorems. and. will be proved in Section. The string equation (.4) plays a particular role in the derivation of the coefficients in the nearest-neighbor recurrence relations. We then perform a numerical study of the coefficients a n, b n in Section. The study suggests that a n+ and b n, n N are all strictly positive, and the limits of a n /n / and b n /n / exist as n, and we can identify these limits explicitly. Section 4 deals with the zeros of P,l. We will give precise location and interlacing results for the zeros of the diagonal multiple orthogonal polynomials P n,n and asymptotic results for the ratio of diagonal multiple orthogonal polynomials. The latter allows us to find the asymptotic distribution of the scaled zeros for these diagonal multiple orthogonal polynomials. The zeros of P,l, with l, have a more interesting structure, which depends on the limit of the ratio /l. We investigate these zeros numerically and end this paper with some conclusions and outloo. 5
6 Proofs. Proof of Proposition. This proposition can be proved by induction on the index n. When n = 0, the relation (.) is obvious, which also gives the initial conditions (.4). Suppose we have β () = b e πi/, (γ () ) = a e πi/, (.) and (a, b ) R for n. From (.4), it follows that thus, (γ () n+ ) = ((γ () n ) + (β () n ) ) On the other hand, the equation (.5) implies that = a n e πi/ b n eπi/ = (b n a n)e πi/, (.) β () n+ = n + (γ n+) () β() n = a n+ = b n a n R. (.) ( n + a n+ b n ) e πi/, (.4) thus b n+ = n + b n R. (.5) a n+ The coupled difference equations (.) (.) are immediate from (.) and (.5). The claim for β n () and (γ n () ) can be proved similarly, we omit the details here.. Proof of Theorem. We shall only prove (.0) since the proof of (.) is similar. We first show that P n,n+m defined in (.0) is a monic polynomial of degree n + m. Observe that ( ) n d n ( ) e x P n dx n 0,m = ( )n n ( d n dx n ( e x P 0,m ) ) = (e x P n,n+m (x)), we then obtain from (.0) the following difference-differential equation for P n,n+m : P n,n+m (x) = x P n,n+m (x) P n,n+m (x). (.6) We can now use induction on n. Clearly P 0,m = p () m is a monic polynomial of degree m. Suppose that P n,m+n is a monic polynomial of degree n + m, then (.6) implies that P n,n+m is a monic polynomial of degree n + m. Next, we show that P n,n+m satisfies the orthogonality conditions (.7) (.8). With Γ 0 defined in (.6), it follows from (.0) and integration by parts times that Γ 0 x P n,n+m (x)e x dx = ( )n n Γ0 x dn dx n ( e x P 0,m (x) ) dx = ( )n+! d n ( e x P n dx n 0,m (x) ) x=0 =! +P n,n+m (0), (.7) 6
7 for = 0,,..., n. Similarly, it is easily seen that x P n,n+m (x)e x dx = x P n,n+m (x)e x dx =! Γ Γ P n,n+m (0). (.8) + Combining (.7) and (.8) gives x P n,n+m (x)e x dx = 0, = 0,,..., n, Γ 0 Γ x P n,n+m (x)e x dx = 0, = 0,,..., n. Γ 0 Γ We still need m more orthogonality conditions to complete (.8), but these follow from ( ) x n+ P n,n+m (x)e x dx = ( )n n+ dn x e Γ 0 Γ Γ0 Γ x P n dx n 0,m (x) dx (n + )! = x e x P! n 0,m (x)dx = 0 Γ 0 Γ for = 0,,..., m, where we used the fact that P 0,m = p () m is the orthogonal polynomial for the cubic exponential weight on Γ 0 Γ.. Proof of Theorem. We will present the proof of (.) (.7), the remaining part of the theorem can be proved in a similar manner. Let us denote the coefficients of x +l and x +l in P,l by δ,l and ε,l, respectively, i.e., P,l (x) = x +l + δ,l x +l + ε,l x +l +. (.9) Substituting the above formula into (.6) and comparing the coefficients of x n+m and x n+m on both sides leads to thus, for m N. Similarly, we have which implies δ n,n+m = δ n,n+m, ε n,n+m = ε n,n+m, δ n,n+m = δ 0,m, ε n,n+m = ε 0,m, (.0) P n+m,n (x) = x P n+m,n (x) P n+m,n (x), (.) δ n+m,n = δ m,0, ε n+m,n = ε m,0. (.) If we insert (.9) into (.) (.), then the coefficients of second leading term x n+m give { δ0,0 δ c n,n+m = δ n,n+m δ n+,n+m =,0, m = 0, (.) δ 0,m δ 0,m, m >, d n,n+m = δ n,n+m δ n,n+m+ = δ 0,m δ 0,m+, (.4) 7
8 where we have also made use of the first equalities in (.0) and (.). On account of the facts that xp m,0 (x) = P m+,0 (x) + β m () P m,0 (x) + (γ m () ) P m,0 (x), (.5) xp 0,m (x) = P 0,m+ (x) + β m () 0,m(x) + (γ m () P 0,m (x), (.6) (see (.9) and (.8)), it is immediate that δ m,0 = δ m+,0 + β m () = δ m+,0 + b m e πi/, (.7) δ 0,m = δ 0,m+ + β m () 0,m+ + b m e πi/, (.8) in view of (.) and (.5). The values for c n,n+m, d n,n+m in (.4) (.5) then follow from combining (.), (.4) and (.7) (.8). We now establish the equalities (.6) (.7) for a n,n+m and b n,n+m. Multiplying both sides of (.) by x n+m e x and integrating the equality over Γ 0 Γ, the orthogonality condition (.8) implies x n+m P n,n+m (x)e x dx Γ b n,n+m = 0 Γ. (.9) x n+m P n,n+m (x)e x dx Γ 0 Γ By (.0), (.) and integrating by parts, we find that x n+m P n,n+m (x)e x dx = ( )n Γ 0 Γ n and Γ 0 Γ x n P n,n (x)e x = = (n + m)! n m! (n + m)! n m! dx = ( )n n = (n )! n ( ) e x P dx n 0,m (x) dx x m P 0,m (x)e x dx Γ 0 Γ P0,m dx, (.0) (x)e x Γ 0 Γ n+m dn x Γ0 Γ ( ) n dn x e dx Γ0 Γ x P n,0 (x) dx P,0 (x)e x dx. (.) Γ 0 Γ Hence, we can simplify (.9) as n Γ 0 Γ e x dx Γ 0 Γ P,0 (x)e x dx, m = 0, b n,n+m = n + m Γ 0 Γ P0,m(x)e x dx m Γ 0 Γ P0,m (x)e x dx, m > 0. (.) Note that Γ 0 Γ P 0,m (x)e x dx = (γ () γ() γ m () ) e x dx, m > 0, (.) Γ 0 Γ 8
9 and straightforward calculations give us dx = Γ(/) ( ω ), (.4) Γ 0 Γ e x P,0 (x)e x dx = (x β () 0 )e x dx Γ 0 Γ Γ 0 Γ = Γ(/) ( ω) Γ(/)β() 0 ( ω ) = Γ(/) ( ω)( ( + ω)e πi/ ). (.5) See (.9) for the value of β () 0. Inserting (.) (.5) into (.), we arrive at n Γ(/) i, m = 0, Γ(/) b n,n+m = n+m m (γ() m ), m > 0, (.6) which is (.7) by (.5). We can also represent a n,n+m as a ratio of two integrals. Indeed, by performing similar strategies above, it is easily seen that x n P n,n+m (x)e x dx Γ a n,n+m = 0 Γ = n P 0,m (x)e x dx Γ 0 Γ. (.7) x n P n,n+m (x)e x dx P 0,m+ (x)e x dx Γ 0 Γ Γ 0 Γ Unfortunately, this representation is not suitable for direct calculation except for m = 0, which gives a n,n = n P 0,0 (x)e x dx Γ 0 Γ P 0, (x)e x dx Γ 0 Γ = n e x dx Γ 0 Γ (x β () 0 )e x dx Γ 0 Γ = b n,n = n Γ(/) i. (.8) Γ(/) For m > 0 the integrals of P 0,m over Γ 0 Γ are involving polynomials orthogonal on the contour Γ 0 Γ, hence it is then difficult to deal with them. So we proceed in another way and we calculate the sum a n,n+m + b n,n+m. Recall the notation δ,l and ε,l in (.9). By comparing the coefficient of x n+m on both sides of (.), it follows from (.0) that a n,n+m + b n,n+m = ε n,n+m ε n+,n+m c n,n+m δ n,n+m = ε 0,m ε 0,m c n,n+m δ 0,m (.9) 9
10 for m > 0. From (.6), we have This, together with (.5), (.4) and (.9), implies ε 0,m = ε 0,m+ + β () m δ 0,m + (γ () m ). (.0) a n,n+m + b n,n+m = ε 0,m ε 0,m c n,n+m δ 0,m = β () m δ 0,m (γ () m ) + β () m δ 0,m = β () m (δ 0,m δ 0,m ) (γ () m ) = (β () m ) (γ () m ) = (γ () m ), m > 0, (.) where we have made use of (.8) in the fourth equality and the string equation (.4) in the last step. A combination of (.6), (.) and (.8) finally gives n Γ(/) i, m = 0, Γ(/) a n,n+m = (.) n m (γ() m ), m > 0, which is (.6), on account of (.5). Asymptotics of a n and b n From Theorem., it is clear that the coefficients in the nearest-neighbor recurrence relations are determined by a n and b n generated from (.) (.4). It is then interesting to study their large n behavior. In Figure we have plotted the values of a n /n / and b n /n / for n from 0 to 70, from which we see that a n+ and b n are all strictly positive for n N. We actually have the following conjecture concerning this observation Figure : The values of a n /n / (left) and b n /n / (right) for n from 0 to 70. Conjecture.. There is a unique positive solution of the recurrence relations (.) (.) with a 0 = 0 and a n+ > 0, b n > 0 for n N. This solution corresponds to the initial condition b 0 = Γ(/)/Γ(/). 0
11 The numerical study further suggests that the limits of a n /n / and b n /n / exist as n, which we can identify in the proposition below. Proposition.. Every positive solution of (.) (.) has the property that lim a n/n / = n /, lim b n/n / = n /. Proof. This can be proved by an argument which was already used by Freud in [, ]. First we show that (a n /n / ) n N is a bounded sequence. From (.) and the positivity of a n+ we find that a n b n. From (.) and the positivity of b n we find a n b n n and thus also 9a n b n n. Together this gives 9a n n, so that 0 a n /n / /9 /. Let a = lim inf n a n /n / and A = lim sup n a n /n /, then 0 a A <. From (.) and the positivity of b n we find b n = a n + a n+. Insert this in (.) to find a n ( an + a n+ + a n + a n ) = n. (.) Let n in (.) through a subsequence for which a n /n / a, then one finds 6a a + A. If n through a subsequence for which a n /n / A, then 6A a + A. Together this gives 6A a + A 6a a + A. If a + A = 0 then one automatically has a = A = 0 so that the limit exists (further on we will see that a 0 so that this case does not happen). If a + A > 0 then one finds A a, and together with a A we see that also in this case a = A and the limit exists. From (.) we then find that lim n b n /n / = a. If we use that information in (.), then 6a a =, so that a = (/6) / = /( / ). The limit for b n /n / follows immediately from this Figure : Zeros of P 45,0 (after scaling) 4 Zeros The formulas in Theorems. and. can be used to generate the multiple orthogonal polynomials P,l defined in (.7) and (.8). We investigate the distribution of their zeros numerically. If one of and l is zero, the polynomials are orthogonal for the exponential
12 Figure 4: Zeros of P 5,5 (after scaling) cubic weight on the curve (Γ 0 Γ or Γ 0 Γ ) in the complex plane. The zeros of P 45,0 (45 / x) are plotted in Figure. It is nown that, in this case, the zeros of the polynomials, after proper scaling, will accumulate on an analytic contour in the complex plane that possesses the so-called S property; cf. [4, 5]. The zero distribution was investigated earlier by Deaño, Huybrechs and Kuijlaars [8], who in fact used the weight e ix. However, a simple rotation x = ye πi/6 is enough to transform their results to the exponential cubic e y which we are using. Suppose that = l = n. It follows from Theorem. that ( P n,n (x) = ( )n dn x e e x). (4.) n dx n We can describe the asymptotic distribution of the zeros of the diagonal polynomials P n,n in more detail. The main reason is that the zeros of P n,n are all located on the three rays Γ 0 Γ Γ, which simplifies matters considerably (see Figure 4). We have the following result for the diagonal polynomials. Observe that this result is the solution of Problem 59 [4, Part V, Chapter ] for the polynomial R n and q =. Proposition 4.. The polynomials P n,n (x) satisfy the symmetry property P n,n (ωx) = ω n P n,n (x), where ω = e πi/ is the primitive third root of unity. In particular n/ j=0 a jx j, n 0 mod, P n,n (x) = x (n )/ j=0 b j x j, n mod, x (4.) (n )/+ j=0 c j x j, n mod, where (a j ) j, (b j ) j, (c j ) j are real sequences depending on n. Furthermore the number of strictly positive real zeros of P n,n is n, if n 0 mod, (n ), if n mod, (n ) +, if n mod, and P n,n (x) has a zero of multiplicity one at x = 0 when n mod and a zero of multiplicity two at x = 0 when n mod.
13 Proof. We use induction on n. The symmetry property follows easily from the Rodrigues formula, so we only need to prove the result about the positive real zeros. Observe that P 0,0 (x) =, P, (x) = x, P, (x) = x(x /), so that the result is true for n = 0,,. Suppose that the result is true for n and let x > x > > x > 0 be the positive real zeros of P n,n. Clearly the sign of P n,n (x j ) is ( ) j+ for j, hence from P n,n (x) = x P n,n (x) P n,n (x) (4.) we find that the sign of P n,n (x j ) is ( ) j, hence P n,n changes sign times and Rolle s theorem guarantees that there are at least zeros y > y > > y with x j < y j < x j, where x 0 = +. If n 0 mod then n mod and the induction hypothesis says that = (n )/ + is odd and P n,n (x) has a zero of multiplicity one at x = 0. The sign of P n,n (0) is ( ) = so that the sign of P n,n (0) is ( ) + =, hence there is also a zero y + of P n,n between 0 and x, giving a total of + = n/ positive real zeros. The ω-symmetry gives another n/ zeros on Γ and n/ zeros on Γ, which is a total of n zeros. Hence there are no other zeros of P n,n. If n mod then n 0 mod and the induction hypothesis gives = (n )/ and P n,n (x) has no zero at x = 0. Hence there will not be an additional zero between 0 and x so that there are = (n )/ positive real zeros for P n,n. There is double zero of P n,n (x) at x = 0. The ω-symmetry gives another (n )/ zeros on Γ and (n )/ zeros on Γ, hence the total number of zeros is (n ) + = n so that there are no other zeros. If n mod then n mod and the induction hypothesis gives = (n )/ is even and a double zero for P n,n (x) at x = 0. Then (4.) implies that P n,n (x) has a single zero at 0. The polynomial P n,n (x)/x of degree n 4 has positive zeros and the sign of this polynomial as x 0 is ( ) =, so that P n,n (x)/x has sign ( ) + = as x 0. Hence P n,n (x)/x has a zero y + between 0 and x, giving a total of + = (n )/+ positive real zeros. The ω-symmetry gives another (n )/ + zeros on Γ and another (n )/ + zeros on Γ, hence together with the single zero at x = 0 this gives a total of (n ) + 4 = n zeros for P n,n so that there are no other zeros. Observe that the proof also shows that the zeros of P n,n and P n,n interlace in the sense that x < y <, x j < y j < x j for j =,...,, x < y < and 0 < y + < x (the latter only when n 0 mod and n mod ). We can now prove the following results Theorem 4.. Let K be a compact set in C \ (Γ 0 Γ Γ ), then lim n P n,n (n / x) n / P n,n (n / x) = Φ(x),
14 holds uniformly for x K, where Furthermore Φ(x) = e πi/ ( holds uniformly for x K. Proof. Consider the ratio N lim n + 9 4x ) / ( + e πi/ P n,n (n/ x) n / P n,n (n / x) = x Φ(x), d P dx n,n(n / x) = P n,n(n / x) P n,n (N / x) N / P n,n (N / x) = N ) /. 9 4x + x n j= x x j,n /N /, where {x j,n, j n} are the zeros of P n,n which are all on the set Γ 0 Γ Γ, then if x K we have P n,n (N/ x) N / P n,n (N / x) n N x x j,n /N / n Nδ, j= where δ = inf{ x y : x K, y Γ 0 Γ Γ } > 0 is the minimal distance between K and Γ 0 Γ Γ. If n/n we then see that the family of analytic functions P n,n (N/ x) N / P n,n (N / x), is uniformly bounded on K. By Montel s theorem there exists a subsequence (n ) N such that P n lim,n (n / x) n n / = F(x), (4.4) P n,n (n / x) uniformly for x K, where F is an analytic function on K for which F(x) = /x+o(/x ) as x. This function F may depend on the selected subsequence, so our aim is to prove that it is independent of the subsequence. Now consider (4.) for P n,n, then N / P n,n (N / x) P n,n (N / x) = x N / P n,n (N/ x) P n,n (N / x), hence (4.4) implies (with N = n ) lim n n / P n,n (n / x) = Φ(x), (4.5) P n,n (n / x) uniformly on K, where Φ(x) = x F(x)/. This uniform convergence of analytic functions implies also the uniform convergence of the derivatives, hence ( ) Φ P n,n (x) = lim (n / x) n n / P n,n (n / ( x) P n,n = lim (n / x) / n P n,n (n / x) n n / P n,n (n / x) P n,n (n / x) n / P n,n (n / x) ) P n,n (n / x), 4
15 but this means that lim n n / P n,n (n / x) = F(x), (4.6) P n,n (n / x) uniformly on K, with the same limit as in (4.4). But then (4.) implies that lim n n / P n,n (n / x) = Φ(x), (4.7) P n,n (n / x) uniformly on K, with the same limit as in (4.5). We can repeat this reasoning once more and find that P n lim,n (n / x) = F(x), (4.8) n and lim n n / n / P n,n (n / x) P n +,n +(n / x) P n,n (n / x) = Φ(x), (4.9) uniformly on K. We will show that the function Φ satisfies a cubic equation, from which we can determine Φ and hence also F uniquely, so that Φ and F do not depend on the selected subsequence (n ) N. Consider the nearest neighbor recurrence relations for the diagonal n = m xp n,n (x) = P n+,n (x) + c n,n P n,n (x) + a n,n P n,n (x) + b n,n P n,n (x) (4.0) xp n,n (x) = P n,n+ (x) + d n,n P n,n (x) + a n,n P n,n (x) + b n,n P n,n (x). (4.) Subtracting (4.0) and (4.) gives P n+,n (x) P n,n+ (x) = (d n,n c n,n )P n,n (x). Use this for n n to eliminate P n,n (x) in (4.0) to find xp n,n (x) = P n+,n (x) + c n,n P n,n (x) + (a n,n + b n,n )P n,n (x) From Theorem. we have + a n,n (c n,n d n,n )P n,n (x). c n,n = Γ(/) Γ(/) eπi/, d n,n = Γ(/) Γ(/) e πi/, so that c n,n d n,n = i Γ(/)/Γ(/). Furthermore a n,n = ni Γ(/) Γ(/), b n,n = ni Γ(/) Γ(/), so that the recurrence relation becomes xp n,n (x) = P n+,n (x) + Γ(/) Γ(/) eπi/ P n,n (x) + n P n,n (x). (4.) 5
16 Use this for n / x and divide by P n,n (n / x), then x = P n+,n (n / x) + Γ(/) n / P n,n (n / x) n / Γ(/) eπi/ + and by using (4.7) we find lim n n / n P n,n (n / x), n / P n,n (n / x) P n +,n (n / x) P n,n (n / x) = x Φ(x), (4.) uniformly on K. We can repeat the reasoning for n n and use (4.5) to find lim n n / P n,n (n / x) P n,n (n / x) = x Φ(x), (4.4) uniformly on K. But then the uniform convergence also holds for the derivative, and as before (4.6) then implies that lim n n / P n,n (n / x) = F(x). P n,n (n / x) Use (4.) for P n+,n (n / x) and divide by P n,n (n / x), then the latter asymptotic result gives P n +,n lim (n / x) n n / P n,n (n / x) = x F(x) = Φ(x), (4.5) uniformly on K. In a similar way as before, the nearest neighbor recurrence relations for (n +, n) can be transformed to xp n+,n (x) = P n+,n+ (x) + d n+,n P n+,n (x) + (a n+,n + b n+,n )P n,n (x) From Theorem. we now use d n+,n = b 0 e πi/, c n+,n = b e πi/, + b n+,n (d n,n c n,n )P n,n (x). so that d n,n c n,n = (b 0 +b )e πi/ = e πi/ /(a ), where we used (.) with n =. We also have a n+,n = (n + )a e πi/, b n+,n = na e πi/, so that the recurrence relation becomes xp n+,n (x) = P n+,n+ (x) Γ(/) Γ(/) eπi/ P n+,n (x) + a e πi/ P n,n (x) + n P n,n (x). (4.6) Consider this for n / x and divide by P n+,n (n / x) then x = P n+,n+ (n / x) Γ(/) n / P n+,n (n / x) n / Γ(/) eπi/ + a P n,n(n / x) n /e πi/ P n+,n (n / x) 6 + n n / P n,n (n / x) P n+,n (n / x)
17 and by using (4.) and (4.5) we find uniformly on K. Now use the relation lim n n / P n +,n +(n / x) P n +,n (n / x) = x Φ(x), (4.7) P n+,n+ (n / x) = P n+,n+ (n / x) n / P n,n (n / x) n / P n+,n (n / x) P n+,n (n / x) n / P n,n (n / x) and let n through the subsequence (n ) N, then (4.9), (4.) and (4.7) show that Φ(x) = ( x ). (4.8) Φ(x) The cubic equation (4.8) has one solution Φ which behaves for x as Φ (x) = x + O(/x), x. There are two other solutions Φ, which behave as /(x) as x Φ (x) = x + 7x 5/ + O(/x4 ), Φ (x) = x 7x 5/ + O(/x4 ), x. Recall that our Φ satisfies Φ(x) = x F(x)/, where F(x) = O(/x), so that we need the solution Φ. The discriminant of (4.8) is (4x 9)/7 so that Φ has branch points at (9/4) /, (9/4) / e ±πi/, which are three points on Γ 0, Γ, Γ respectively, and since all the zeros of P n,n are on Γ 0 Γ Γ, we conclude that the scaled zeros x j,n /n / are dense on the three segments [0, (9/4) / ] [0, (9/4) / e πi/ ] [0, (9/4) / e πi/ ]. The cubic equation can be solved explicitly by using Cardano s formula: let y = x /(Φ) and z = /y, then the cubic equation (4.8) becomes z xz + = 0, and the solutions are z = ω j u / +ω j v / (j = 0,, ), where ω = e πi/, u +v = and uv = x, i.e., u = + 9 4x, v = 9 4x = x + 9 4x. The solution Φ corresponds to the solution with z(x) = /x + O(/x 4 ), and this is z(x) = ω u / + ω v /, and since Φ = y, we find Φ (x) = (ω u / + ω v / ) = ωu / + ω v / + x. 7
18 Corollary 4.. Let {x j,n, j =,,..., n} be the zeros of P n,n and µ n be the normalized counting measure of the scaled zeros x j,n /n /, µ n = n n δ xj,n /n /. Then the sequence (µ n ) n converges wealy to the probability measure µ for which f(x)dµ(x) = where ω = e πi/ and with (9/4) / 0 v(x) = v(x)f(x) dx + j= ω(9/4) / 0 v(x)f(x) dx + ω (9/4) / 0 v(x)f(x) dx, ( ) + x[a(x) + b(x)] [b(x) a(x)], (4.9) 4π ( )/ 9 4x a(x) =, b(x) = Proof. The Stieltjes transform of the measure µ n is ( + 9 4x x t dµ n(t) = P n,n(n / x) n / P n,n (n / x), )/. (4.0) hence Theorem 4. gives lim n x t dµ n(t) = F(x) = ( ) x Φ(x), uniformly on compact sets of C \ (Γ 0 Γ Γ ). The Grommer-Hamburger theorem [] then implies that µ n converges wealy to a measure µ for which x t dµ(t) = ( ) x Φ(x). The function Φ is analytic in C \ ([0, (9/4) / ] [0, ω(9/4) / ] [0, ω (9/4) / ]), hence the measure µ is supported on [0, (9/4) / ] [0, ω(9/4) / ] [0, ω (9/4) / ]. Furthermore it is absolutely continuous and we can find the density by using the Stieltjes inversion formula v(x) = ( ) π lim I (x + iɛ) Φ(x + iɛ). ɛ 0+ Due to the ω-symmetry, it is sufficient to determine v(x) for x [0, (9/4) / ]. Clearly v(x) = π lim IΦ(x + iɛ) = ɛ 0+ π I (ω a + ω b), with a and b given in (4.0). Then by some elementary (complex) calculus, using (a ab +b )(a + b) = a +b and ab = x, one finds the expression (4.9) for the density v. 8
19 x~ Figure 5: Histogram of the real zeros of P 600,600 and the density v(x) on [0, (9/4) / ] In Figure 5 we have given a histogram of the 400 real zeros of P 600,600 together with the density v, scaled so as to have total mass one for all the real zeros. There are 400 zeros on the interval [0, ω(9/4) / ] and 400 zeros on the interval [0, ω (9/4) / ] and these zeros are obtained by rotating the real zeros over an angle ±π/. The density v has a finite non-zero value at the origin v(0) = / /4π = and tends to zero as (9/4) / x when x (9/4) /. If l, the rotational symmetry of the zeros is broen. Suppose that l, then we see numerically that zeros of P,l lie on the line containing Γ (some zeros are in fact on Γ ), while the other zeros are distributed on a complex contour in the lower half plane; see Figure 6. Similarly, if l, then l zeros of P,l lie on the line containing Γ (again some zeros are on Γ ), and the other zeros are distributed on a complex contour in the upper half plane, as illustrated in Figure 7. Indeed, from Theorem., it is easily seen that the zeros of P,l are complex conjugates of the zeros of P l,. We expect that the asymptotic zero distribution of P,l will depend on the limit of the ratio /l. Figure 6: Zeros of P 7,0 (left) and P 4,5 (right) after scaling 9
20 Figure 7: Zeros of P 0,7 (left) and P 6,4 (right) after scaling 5 Conclusions and outloo In this paper, we have introduced the multiple orthogonal polynomials associated with an exponential cubic weight e x over two contours in the complex plane. The basic properties of these polynomials are studied, which include the Rodrigues formula and nearest-neighbor recurrence relations. These results then allow us to perform numerical studies of the recurrence coefficients and zero distributions of the multiple orthogonal polynomials. Moreover, the recurrence coefficients are related to a discrete Painlevé equation. One can also consider the more general exponential cubic weight e x +tx, where t R, and the associated multiple orthogonal polynomials have similar Rodrigues formulas and nearest-neighbor recurrence relations. Indeed, with e x replaced by e x +tx in Theorem., the difference-differential equations (.6) and (.) now read This then implies that P n,n+m (x) = (x t )P n,n+m (x) P n,n+m (x), P n+m,n (x) = (x t )P n+m,n (x) P n+m,n (x). δ n,n+m = δ 0,m, ε n,n+m = ε 0,m t, δ n+m,n = δ m,0, ε n+m,n = ε m,0 t, where δ,l and ε,l are defined in (.9) and now depend on t. Following the same strategy as in the proof of Theorem. and using the string equation (A.7) at the final stage, we obtain (.) and (.) with { β () 0 (t), m = 0, c n,n+m = β () m d n,n+m = β m () (t), a n,n+m = b n,n+m = (t), m > 0, { n β () 0 (t) β() 0 n m { (γ() n β () 0 (t) β() 0 (n+m) m (γ() (t), m = 0, m (t)) t, m > 0, 0 (t), m = 0, m (t)), m > 0,
21 where β n () (t), (γ n () ) (t) are the recurrence coefficients of the monic orthogonal polynomials with respect to the weight e x +tx on Γ 0 Γ, and β n () (t), (γ n () ) (t) are the recurrence coefficients of the monic orthogonal polynomials with respect to the weight e x +tx on Γ 0 Γ. The recurrence relations (.8) and (.9) for this general case can be obtained similarly, but we omit the results here. Clearly, in the general case, we lose the nice structure of the recurrence coefficients as stated in Theorem., and more importantly we lose the symmetry given in Proposition 4., which is why we focus on the weight e x in this paper. The challenging problem is to establish the asymptotic zero distribution of P,l (x) for the non-symmetric case. At present we are unable to find an analogue of Theorem 4. and Corollary 4. because of two reasons: first one needs the asymptotic behavior of the recurrence coefficients and at present we can only conjecture the behavior (see Proposition.). If we assume this to be correct, then the proof of Theorem 4. can be used to find the asymptotic behavior, away from the set where the zeros of the multiple orthogonal polynomials accumulate, in terms of an algebraic function Φ satisfying a cubic equation. But the second reason is that we don t now where the zeros of the multiple orthogonal polynomials accumulate. The discriminant of the cubic equation is a quartic polynomial in x and the four roots are branch points of the algebraic function Φ. The zeros will accumulate on two curves, each connecting two points in the complex plane, see Figures 6 and 7. One of the curves is a straight line, the other is a curved line connecting two of the four branch points. The straight line, however, does not connect the other two branch points but starts from one branch point and stops before the second branch point is reached. This suggests that a vector equilibrium problem is involved, for two measures living on curves connecting four branch points, with an external field x induced by the weight e x. In order to characterize the limiting zero distribution, one may need to extend the concept of S-property (cf. [4, 5]) for orthogonal polynomials and equilibrium measures to this setting for multiple orthogonal polynomials and vector equilibrium problems. Once that is obtained, one may be able to use the Riemann-Hilbert method to find the asymptotic behavior of the multiple orthogonal polynomials. Appendix A Derivation of the string equations In this appendix, we give an alternative proof of the string equations (.4) (.5) using ladder operators for orthogonal polynomials. Note that the ladder operators for multiple orthogonal polynomials and their compatibility conditions can be found in [0]. Following the general set-up (cf. [5]), if the weight function w vanishes at the endpoints of the orthogonality interval, the lowering and raising ladder operators for the associated monic polynomials p n are given by ( ) d dx + B n(x) p n (x) = γn A n(x)p n (x), (A.) ( ) d dx B n(x) v (x) p n (x) = A n (x)p n (x), (A.)
22 with and where v(x) := lnw(x), A n (x) := v (x) v (y) [p n (y)] w(y)dy, (A.) h n x y B n (x) := v (x) v (y) p n (y)p n (y)w(y)dy, (A.4) h n x y p m (x)p n (x)ω(x)dx = h n δ m,n, m, n = 0,,,.... (A.5) Note that A n and B n are not independent, but satisfy the following compatibility conditions [8, Lemma.. and Theorem..4]. Proposition A.. The functions A n and B n defined in (A.) and (A.4) satisfy B n+ (x) + B n (x) = (x β n )A n (x) v (x), (S ) + (x β n )[B n+ (x) B n (x)] = γ n+a n+ (x) γ na n (x). (S ) Now we consider a more general exponential cubic weight e x +tx, with parameter t R. Then v(x) = lnw(x) = x tx, and It then follows from (A.) (A.5) that v (x) v (y) x y = (x + y). A n (x) = (x + β n ), B n (x) = γ n. (A.6) Substituting (A.6) into (S ) and comparing the coefficients of the constant term, we have γ n + γ n+ + β n t = 0. (A.7) From (S ) we similarly get γ n(β n + β n ) = n. (A.8) Note that in this case, the recurrence coefficients β n and γ n all depend on t. By setting t = 0 in (A.7) and (A.8), we recover the string equations (.4) and (.5). The weight e x +tx is a modification of the weight e x with an exponential factor e tx, and as a consequence the recurrence coefficients satisfy the Toda equations [8,.8] d dt γ n = γ n (β n β n ), d dt β n = γn+ γn. If we differentiate (A.0) and then use (A.9), we find β n(t) = γ n+(β n+ β n ) γ n(β n β n ). (A.9) (A.0)
23 Then use (A.7) and (A.8) to find β n (t) = β n t β n + n +, which is the Painlevé II equation. The equations (A.7) and (A.8) give n + n + + βn β n + β n β n+ + β = t n which also follows from the Bäclund transformation of the second Painlevé equation (see [] and [7]), hence it is not surprising to find that β n satisfies the Painlevé II equation. If we write x n (t) = aβ n ( at), where a = (/) /, then x n(t) = x n + tx n + n +, which is the Painlevé II equation in standard form and with parameter α = n +. The second Painlevé equation is closely related to the Airy equation and has special solutions in terms of Airy functions for the parameter values α = n +, with n Z [6, 7.]. This relation with the Airy function was to be expected since one has the integral representation Ai(t) = e πi/ e z / tz dz, πi e πi/ (see Eq of the NIST Digital Library of Mathematical Functions ) which contains (a slight variation of) the weight e z +tz. The special solution of Painlevé II in terms of Airy functions is x n (t) = τ n(t) τ n (t) τ n+ (t) τ n+ (t), where τ n is the Hanel matrix ϕ(t) ϕ (t) ϕ (n ) (t) ϕ (t) ϕ (t) ϕ (n) (t) τ n (t) =.. ϕ (n ) (t) ϕ (n) (t)., ϕ (n ) (t) and ϕ is a solution of the Airy equation ϕ + tϕ = 0. This solution of P II coincides with the well nown solution β n = δ n δ n+, where δ n is the coefficient of x n for the monic orthogonal polynomial P n (x) = x n + δ n x n +. One has δ n = n / n, where n is the Hanel matrix containing the moments m 0 m m n m m m n n =... m n m n m n
24 and n is a similar determinant but with the last column replaced by m n, m n+,..., m n respectively. Even though this is an explicit solution, it is not convenient for finding the recurrence coefficients when n is large because of the high number of computations involved, whereas the relations (A.7) (A.8) have a low computational complexity. The explicit solution is also not convenient for obtaining the asymptotic behavior of β n and a n as n, which is easier to obtain from the string equations (see Proposition.). Acnowledgements The authors are grateful to the referees for their careful reading and constructive suggestions. WVA is supported by KU Leuven research grant OT//07, FWO research grant G.094. and the Belgian Interuniversity Attraction Poles Programme P7/8. GF is supported by the MNiSzW Iuventus Plus grant Nr 04/IP/0/7. LZ was a Postdoctoral Fellow of FWO, and is also partially supported by The Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning (No. SHH4007) and by Grant SGST DZ 7800 from Fudan University. References [] A.I. Aptearev, Multiple orthogonal polynomials, J. Comput. Appl. Math. 99 (998), [] A.I. Aptearev, A. Branquinho and W. Van Assche, Multiple orthogonal polynomials for classical weights, Trans. Amer. Math. Soc. 55 (00), [] Y. Ben Cheih and N. Ben Romdhane, On d-symmetric classical d-orthogonal polynomials, J. Comput. Appl. Math. 6 (0), [4] P. Bleher and A. Deaño, Topological expansion in the cubic random matrix model, Int. Math. Res. Not. 0, no., [5] Y. Chen and M. Ismail, Jacobi polynomials from compatibility conditions, Proc. Amer. Math. Soc. (005), [6] P.A. Clarson, Painlevé equations Nonlinear special functions, in Orthogonal Polynomials and Special Functions: Computation and Applications (F. Marcellán, W. Van Assche, Eds.), Lecture Notes in Mathematics 88, Springer-Verlag, Berlin, 006, pp. 4. [7] P. Clarson, E.L. Mansfield, H.N. Webster, On the relation between discrete and continuous Painlevé equations, Theoret. and Math. Phys. (000), no., 6. [8] A. Deaño, D. Huybrechs and A.B.J. Kuijlaars, Asymptotic zero distribution of complex orthogonal polynomials associated with Gaussian quadrature, J. Approx. Theory 6 (00), 0 4. [9] D. Dominici, Asymptotic analysis of generalized Hermite polynomials, Analysis 8 (008),
25 [0] G. Filipu, W. Van Assche and L. Zhang, Ladder operators and differential equations for multiple orthogonal polynomials, J. Phys. A: Math. Theor., 46 (0), 0504, 4 pp. [] A.S. Foas, B. Grammaticos and A. Ramani, From continuous to discrete Painlevé equations, J. Math. Anal. Appl. 80 (99), no., [] G. Freud, On the coefficients in the recursion formulae of orthogonal polynomials, Proc. Royal Irish Acad. Sect. A 76 (976), 6. [] J.S. Geronimo and T.P. Hill, Necessary and sufficient condition that the limit of Stieltjes transforms is a Stieltjes transform, J. Approx. Theory (00), no., [4] A.A. Gonchar and E.A. Rahmanov, Equilibrium distributions and the rate of rational approximation of analytic functions, Mat. Sb. (N.S.) 4(76) (987), 06 5, 447; translation in Math. USSR-Sb. 6 (989), no., [5] H.W. Gould and A.T. Hopper, Operational formulas connected with two generalizations of Hermite polynomials, Due Math. J. 9 (96), 5 6. [6] B. Grammaticos and A. Ramani, The hunting for the discrete Painlevé equations, Regul. Chaotic Dyn. 5 (000), [7] B. Grammaticos and A. Ramani, Discrete Painlevé equations: a review, Lect. Notes Phys. 644 (004), 45. [8] M.E.H. Ismail, Classical and Quantum Orthogonal Polynomials in One Variable, Encyclopedia of Mathematics and its Applications 98, Cambridge University Press, 005. [9] A.P. Magnus, Painlevé type differential equations for the recurrence coefficients of semi-classical orthogonal polynomials, J. Comput. Appl. Math. 57 (995), 5 7. [0] F. Nijhoff, J. Satsuma, K. Kajiwara, B. Grammaticos and A. Ramani, A study of the alternate discrete Painlevé II equation, Inverse Problems (996), [] E.M. Niishin and V.N. Soroin, Rational Approximations and Orthogonality, in: Translations of Mathematical Monographs 9, Amer. Math. Soc. Providence RI, 99. [] R.B. Paris, The asymptotics of the generalised Hermite-Bell polynomials, J. Comput. Appl. Math. (009), 6 6. [] G. Pólya, Über die Nullstellen suzessiver Derivierten, Math. Z. (9), no., [4] G. Pólya and G. Szegő, Problems and Theorems in Analysis II, Springer-Verlag, Berlin, 976 (revised and enlarged translation of Ausgaben und Lehrsätze aus der Analysis II, 4th edition, 97). 5
26 [5] H. Stahl, Orthogonal polynomials with complex-valued weight function. I, II, Constr. Approx. (986), 5 40, 4 5. [6] W. Van Assche, Discrete Painlevé equations for recurrence coefficients of orthogonal polynomials, in Difference Equations, Special Functions and Orthogonal Polynomials (S. Elaydi et al., Eds.), World Scientific, 007, pp [7] W. Van Assche, Nearest neighbor recurrence relations for multiple orthogonal polynomials, J. Approx. Theory 6 (0), [8] W. Van Assche and E. Coussement, Some classical multiple orthogonal polynomials, Numerical analysis 000, Vol. V, Quadrature and orthogonal polynomials. J. Comput. Appl. Math. 7 (00),
Multiple Orthogonal Polynomials
Summer school on OPSF, University of Kent 26 30 June, 2017 Introduction For this course I assume everybody is familiar with the basic theory of orthogonal polynomials: Introduction For this course I assume
More informationSimultaneous Gaussian quadrature for Angelesco systems
for Angelesco systems 1 KU Leuven, Belgium SANUM March 22, 2016 1 Joint work with Doron Lubinsky Introduced by C.F. Borges in 1994 Introduced by C.F. Borges in 1994 (goes back to Angelesco 1918). Introduced
More informationRecurrence Coef f icients of a New Generalization of the Meixner Polynomials
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 7 (011), 068, 11 pages Recurrence Coef f icients of a New Generalization of the Meixner Polynomials Galina FILIPUK and Walter VAN ASSCHE
More informationJacobi-Angelesco multiple orthogonal polynomials on an r-star
M. Leurs Jacobi-Angelesco m.o.p. 1/19 Jacobi-Angelesco multiple orthogonal polynomials on an r-star Marjolein Leurs, (joint work with Walter Van Assche) Conference on Orthogonal Polynomials and Holomorphic
More informationDifference Equations for Multiple Charlier and Meixner Polynomials 1
Difference Equations for Multiple Charlier and Meixner Polynomials 1 WALTER VAN ASSCHE Department of Mathematics Katholieke Universiteit Leuven B-3001 Leuven, Belgium E-mail: walter@wis.kuleuven.ac.be
More informationLAMÉ DIFFERENTIAL EQUATIONS AND ELECTROSTATICS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Xxxx XXXX, Pages 000 000 S 0002-9939(XX)0000-0 LAMÉ DIFFERENTIAL EQUATIONS AND ELECTROSTATICS DIMITAR K. DIMITROV AND WALTER VAN ASSCHE
More informationMultiple orthogonal polynomials. Bessel weights
for modified Bessel weights KU Leuven, Belgium Madison WI, December 7, 2013 Classical orthogonal polynomials The (very) classical orthogonal polynomials are those of Jacobi, Laguerre and Hermite. Classical
More informationOPSF, Random Matrices and Riemann-Hilbert problems
OPSF, Random Matrices and Riemann-Hilbert problems School on Orthogonal Polynomials in Approximation Theory and Mathematical Physics, ICMAT 23 27 October, 2017 Plan of the course lecture 1: Orthogonal
More informationPositivity of Turán determinants for orthogonal polynomials
Positivity of Turán determinants for orthogonal polynomials Ryszard Szwarc Abstract The orthogonal polynomials p n satisfy Turán s inequality if p 2 n (x) p n 1 (x)p n+1 (x) 0 for n 1 and for all x in
More informationUNIFORM BOUNDS FOR BESSEL FUNCTIONS
Journal of Applied Analysis Vol. 1, No. 1 (006), pp. 83 91 UNIFORM BOUNDS FOR BESSEL FUNCTIONS I. KRASIKOV Received October 8, 001 and, in revised form, July 6, 004 Abstract. For ν > 1/ and x real we shall
More informationMultiple Orthogonal Polynomials
Summer school on OPSF, University of Kent 26 30 June, 2017 Plan of the course Plan of the course lecture 1: Definitions + basic properties Plan of the course lecture 1: Definitions + basic properties lecture
More informationP -adic root separation for quadratic and cubic polynomials
P -adic root separation for quadratic and cubic polynomials Tomislav Pejković Abstract We study p-adic root separation for quadratic and cubic polynomials with integer coefficients. The quadratic and reducible
More informationMultiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv: v2 [math.ca] 8 Nov 2013
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv:1310.0982v2 [math.ca] 8 Nov 2013 F Ndayiragie 1 and W Van Assche 2 1 Département de Mathématiques, Université du Burundi, Campus
More informationEntropy of Hermite polynomials with application to the harmonic oscillator
Entropy of Hermite polynomials with application to the harmonic oscillator Walter Van Assche DedicatedtoJeanMeinguet Abstract We analyse the entropy of Hermite polynomials and orthogonal polynomials for
More informationPainlevé VI and Hankel determinants for the generalized Jacobi weight
arxiv:98.558v2 [math.ca] 3 Nov 29 Painlevé VI and Hankel determinants for the generalized Jacobi weight D. Dai Department of Mathematics, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationUltraspherical moments on a set of disjoint intervals
Ultraspherical moments on a set of disjoint intervals arxiv:90.049v [math.ca] 4 Jan 09 Hashem Alsabi Université des Sciences et Technologies, Lille, France hashem.alsabi@gmail.com James Griffin Department
More informationMath 324 Summer 2012 Elementary Number Theory Notes on Mathematical Induction
Math 4 Summer 01 Elementary Number Theory Notes on Mathematical Induction Principle of Mathematical Induction Recall the following axiom for the set of integers. Well-Ordering Axiom for the Integers If
More informationOPSF, Random Matrices and Riemann-Hilbert problems
OPSF, Random Matrices and Riemann-Hilbert problems School on Orthogonal Polynomials in Approximation Theory and Mathematical Physics, ICMAT 23 27 October, 207 Plan of the course lecture : Orthogonal Polynomials
More informationAsymptotics of Integrals of. Hermite Polynomials
Applied Mathematical Sciences, Vol. 4, 010, no. 61, 04-056 Asymptotics of Integrals of Hermite Polynomials R. B. Paris Division of Complex Systems University of Abertay Dundee Dundee DD1 1HG, UK R.Paris@abertay.ac.uk
More informationOn rational approximation of algebraic functions. Julius Borcea. Rikard Bøgvad & Boris Shapiro
On rational approximation of algebraic functions http://arxiv.org/abs/math.ca/0409353 Julius Borcea joint work with Rikard Bøgvad & Boris Shapiro 1. Padé approximation: short overview 2. A scheme of rational
More informationarxiv: v1 [cs.it] 12 Jun 2016
New Permutation Trinomials From Niho Exponents over Finite Fields with Even Characteristic arxiv:606.03768v [cs.it] 2 Jun 206 Nian Li and Tor Helleseth Abstract In this paper, a class of permutation trinomials
More informationConstrained Leja points and the numerical solution of the constrained energy problem
Journal of Computational and Applied Mathematics 131 (2001) 427 444 www.elsevier.nl/locate/cam Constrained Leja points and the numerical solution of the constrained energy problem Dan I. Coroian, Peter
More informationLegendre s Equation. PHYS Southern Illinois University. October 18, 2016
Legendre s Equation PHYS 500 - Southern Illinois University October 18, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 18, 2016 1 / 11 Legendre s Equation Recall We are trying
More informationMath 259: Introduction to Analytic Number Theory More about the Gamma function
Math 59: Introduction to Analytic Number Theory More about the Gamma function We collect some more facts about Γs as a function of a complex variable that will figure in our treatment of ζs and Ls, χ.
More informationViscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces
Viscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces YUAN-HENG WANG Zhejiang Normal University Department of Mathematics Yingbing Road 688, 321004 Jinhua
More information4. Ergodicity and mixing
4. Ergodicity and mixing 4. Introduction In the previous lecture we defined what is meant by an invariant measure. In this lecture, we define what is meant by an ergodic measure. The primary motivation
More informationOn Some Mean Value Results for the Zeta-Function and a Divisor Problem
Filomat 3:8 (26), 235 2327 DOI.2298/FIL6835I Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat On Some Mean Value Results for the
More informationSUMS OF ENTIRE FUNCTIONS HAVING ONLY REAL ZEROS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 SUMS OF ENTIRE FUNCTIONS HAVING ONLY REAL ZEROS STEVEN R. ADAMS AND DAVID A. CARDON (Communicated
More informationReferences 167. dx n x2 =2
References 1. G. Akemann, J. Baik, P. Di Francesco (editors), The Oxford Handbook of Random Matrix Theory, Oxford University Press, Oxford, 2011. 2. G. Anderson, A. Guionnet, O. Zeitouni, An Introduction
More information#A5 INTEGERS 18A (2018) EXPLICIT EXAMPLES OF p-adic NUMBERS WITH PRESCRIBED IRRATIONALITY EXPONENT
#A5 INTEGERS 8A (208) EXPLICIT EXAMPLES OF p-adic NUMBERS WITH PRESCRIBED IRRATIONALITY EXPONENT Yann Bugeaud IRMA, UMR 750, Université de Strasbourg et CNRS, Strasbourg, France bugeaud@math.unistra.fr
More informationA NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS. Masaru Nagisa. Received May 19, 2014 ; revised April 10, (Ax, x) 0 for all x C n.
Scientiae Mathematicae Japonicae Online, e-014, 145 15 145 A NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS Masaru Nagisa Received May 19, 014 ; revised April 10, 014 Abstract. Let f be oeprator monotone
More informationOn integral representations of q-gamma and q beta functions
On integral representations of -gamma and beta functions arxiv:math/3232v [math.qa] 4 Feb 23 Alberto De Sole, Victor G. Kac Department of Mathematics, MIT 77 Massachusetts Avenue, Cambridge, MA 239, USA
More informationSOLUTION TO RUBEL S QUESTION ABOUT DIFFERENTIALLY ALGEBRAIC DEPENDENCE ON INITIAL VALUES GUY KATRIEL PREPRINT NO /4
SOLUTION TO RUBEL S QUESTION ABOUT DIFFERENTIALLY ALGEBRAIC DEPENDENCE ON INITIAL VALUES GUY KATRIEL PREPRINT NO. 2 2003/4 1 SOLUTION TO RUBEL S QUESTION ABOUT DIFFERENTIALLY ALGEBRAIC DEPENDENCE ON INITIAL
More informationDaniel M. Oberlin Department of Mathematics, Florida State University. January 2005
PACKING SPHERES AND FRACTAL STRICHARTZ ESTIMATES IN R d FOR d 3 Daniel M. Oberlin Department of Mathematics, Florida State University January 005 Fix a dimension d and for x R d and r > 0, let Sx, r) stand
More informationP -ADIC ROOT SEPARATION FOR QUADRATIC AND CUBIC POLYNOMIALS. Tomislav Pejković
RAD HAZU. MATEMATIČKE ZNANOSTI Vol. 20 = 528 2016): 9-18 P -ADIC ROOT SEPARATION FOR QUADRATIC AND CUBIC POLYNOMIALS Tomislav Pejković Abstract. We study p-adic root separation for quadratic and cubic
More informationReal Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi
Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.
More informationGALOIS GROUPS OF CUBICS AND QUARTICS (NOT IN CHARACTERISTIC 2)
GALOIS GROUPS OF CUBICS AND QUARTICS (NOT IN CHARACTERISTIC 2) KEITH CONRAD We will describe a procedure for figuring out the Galois groups of separable irreducible polynomials in degrees 3 and 4 over
More informationMATH 324 Summer 2011 Elementary Number Theory. Notes on Mathematical Induction. Recall the following axiom for the set of integers.
MATH 4 Summer 011 Elementary Number Theory Notes on Mathematical Induction Principle of Mathematical Induction Recall the following axiom for the set of integers. Well-Ordering Axiom for the Integers If
More informationOn the Lebesgue constant of barycentric rational interpolation at equidistant nodes
On the Lebesgue constant of barycentric rational interpolation at equidistant nodes by Len Bos, Stefano De Marchi, Kai Hormann and Georges Klein Report No. 0- May 0 Université de Fribourg (Suisse Département
More informationSome tight polynomial-exponential lower bounds for an exponential function
Some tight polynomial-exponential lower bounds for an exponential function Christophe Chesneau To cite this version: Christophe Chesneau. Some tight polynomial-exponential lower bounds for an exponential
More informationMATH 117 LECTURE NOTES
MATH 117 LECTURE NOTES XIN ZHOU Abstract. This is the set of lecture notes for Math 117 during Fall quarter of 2017 at UC Santa Barbara. The lectures follow closely the textbook [1]. Contents 1. The set
More informationA Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)!
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 8 (2005), Article 05.2.3 A Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)! Ira M. Gessel 1 and Guoce Xin Department of Mathematics Brandeis
More informationULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR ORTHOGONAL POLYNOMIALS
ULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR ORTHOGONAL POLYNOMIALS arxiv:083666v [mathpr] 8 Jan 009 Abstract We characterize, up to a conjecture, probability distributions of finite all order moments
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematics ASYMPTOTIC EXPANSION OF THE EQUIPOISE CURVE OF A POLYNOMIAL INEQUALITY ROGER B. EGGLETON AND WILLIAM P. GALVIN Department of Mathematics Illinois
More informationF (z) =f(z). f(z) = a n (z z 0 ) n. F (z) = a n (z z 0 ) n
6 Chapter 2. CAUCHY S THEOREM AND ITS APPLICATIONS Theorem 5.6 (Schwarz reflection principle) Suppose that f is a holomorphic function in Ω + that extends continuously to I and such that f is real-valued
More informationRandom Bernstein-Markov factors
Random Bernstein-Markov factors Igor Pritsker and Koushik Ramachandran October 20, 208 Abstract For a polynomial P n of degree n, Bernstein s inequality states that P n n P n for all L p norms on the unit
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationbe the set of complex valued 2π-periodic functions f on R such that
. Fourier series. Definition.. Given a real number P, we say a complex valued function f on R is P -periodic if f(x + P ) f(x) for all x R. We let be the set of complex valued -periodic functions f on
More informationB n (x) zn n! n=0. E n (x) zn n! n=0
UDC 517.9 Q.-M. Luo Chongqing Normal Univ., China) q-apostol EULER POLYNOMIALS AND q-alternating SUMS* q-полiноми АПОСТОЛА ЕЙЛЕРА ТА q-знакозмiннi СУМИ We establish the basic properties generating functions
More informationRiemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,
Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of
More informationarxiv: v3 [math.ra] 10 Jun 2016
To appear in Linear and Multilinear Algebra Vol. 00, No. 00, Month 0XX, 1 10 The critical exponent for generalized doubly nonnegative matrices arxiv:1407.7059v3 [math.ra] 10 Jun 016 Xuchen Han a, Charles
More informationNonstationary Subdivision Schemes and Totally Positive Refinable Functions
Nonstationary Subdivision Schemes and Totally Positive Refinable Functions Laura Gori and Francesca Pitolli December, 2007 Abstract In this paper we construct a class of totally positive refinable functions,
More informationChapter 8. P-adic numbers. 8.1 Absolute values
Chapter 8 P-adic numbers Literature: N. Koblitz, p-adic Numbers, p-adic Analysis, and Zeta-Functions, 2nd edition, Graduate Texts in Mathematics 58, Springer Verlag 1984, corrected 2nd printing 1996, Chap.
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 9, 1999, pp. 39-52. Copyright 1999,. ISSN 1068-9613. ETNA QUADRATURE FORMULAS FOR RATIONAL FUNCTIONS F. CALA RODRIGUEZ, P. GONZALEZ VERA, AND M. JIMENEZ
More informationINTRODUCTION TO REAL ANALYTIC GEOMETRY
INTRODUCTION TO REAL ANALYTIC GEOMETRY KRZYSZTOF KURDYKA 1. Analytic functions in several variables 1.1. Summable families. Let (E, ) be a normed space over the field R or C, dim E
More informationThe Continuing Story of Zeta
The Continuing Story of Zeta Graham Everest, Christian Röttger and Tom Ward November 3, 2006. EULER S GHOST. We can only guess at the number of careers in mathematics which have been launched by the sheer
More informationAn Approach to Hensel s Lemma
Irish Math. Soc. Bulletin 47 (2001), 15 21 15 An Approach to Hensel s Lemma gary mcguire Abstract. Hensel s Lemma is an important tool in many ways. One application is in factoring polynomials over Z.
More informationTHE NEARLY ADDITIVE MAPS
Bull. Korean Math. Soc. 46 (009), No., pp. 199 07 DOI 10.4134/BKMS.009.46..199 THE NEARLY ADDITIVE MAPS Esmaeeil Ansari-Piri and Nasrin Eghbali Abstract. This note is a verification on the relations between
More informationWe are going to discuss what it means for a sequence to converge in three stages: First, we define what it means for a sequence to converge to zero
Chapter Limits of Sequences Calculus Student: lim s n = 0 means the s n are getting closer and closer to zero but never gets there. Instructor: ARGHHHHH! Exercise. Think of a better response for the instructor.
More informationON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION
Annales Univ. Sci. Budapest., Sect. Comp. 33 (2010) 273-284 ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION L. László (Budapest, Hungary) Dedicated to Professor Ferenc Schipp on his 70th
More information7 Asymptotics for Meromorphic Functions
Lecture G jacques@ucsd.edu 7 Asymptotics for Meromorphic Functions Hadamard s Theorem gives a broad description of the exponential growth of coefficients in power series, but the notion of exponential
More informationThe chain rule for F-differentiation
Irish Math. Soc. Bulletin Number 77, Summer 2016, 19 34 ISSN 0791-5578 The chain rule for -differentiation T. CHAOBANKOH, J.. EINSTEIN AND S. MORLEY Abstract. Let X be a perfect, compact subset of the
More informationMATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5
MATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5.. The Arzela-Ascoli Theorem.. The Riemann mapping theorem Let X be a metric space, and let F be a family of continuous complex-valued functions on X. We have
More informationEQUIVALENCE OF TOPOLOGIES AND BOREL FIELDS FOR COUNTABLY-HILBERT SPACES
EQUIVALENCE OF TOPOLOGIES AND BOREL FIELDS FOR COUNTABLY-HILBERT SPACES JEREMY J. BECNEL Abstract. We examine the main topologies wea, strong, and inductive placed on the dual of a countably-normed space
More informationFourier series of sums of products of Bernoulli functions and their applications
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 7, 798 85 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa Fourier series of sums of products
More informationComputing a Lower Bound for the Canonical Height on Elliptic Curves over Q
Computing a Lower Bound for the Canonical Height on Elliptic Curves over Q John Cremona 1 and Samir Siksek 2 1 School of Mathematical Sciences, University of Nottingham, University Park, Nottingham NG7
More informationZeros of Polynomials: Beware of Predictions from Plots
[ 1 / 27 ] University of Cyprus Zeros of Polynomials: Beware of Predictions from Plots Nikos Stylianopoulos a report of joint work with Ed Saff Vanderbilt University May 2006 Five Plots Fundamental Results
More informationPADÉ INTERPOLATION TABLE AND BIORTHOGONAL RATIONAL FUNCTIONS
ELLIPTIC INTEGRABLE SYSTEMS PADÉ INTERPOLATION TABLE AND BIORTHOGONAL RATIONAL FUNCTIONS A.S. ZHEDANOV Abstract. We study recurrence relations and biorthogonality properties for polynomials and rational
More informationPolyexponentials. Khristo N. Boyadzhiev Ohio Northern University Departnment of Mathematics Ada, OH
Polyexponentials Khristo N. Boyadzhiev Ohio Northern University Departnment of Mathematics Ada, OH 45810 k-boyadzhiev@onu.edu 1. Introduction. The polylogarithmic function [15] (1.1) and the more general
More informationSYMMETRY AND SPECIALIZABILITY IN THE CONTINUED FRACTION EXPANSIONS OF SOME INFINITE PRODUCTS
SYMMETRY AND SPECIALIZABILITY IN THE CONTINUED FRACTION EXPANSIONS OF SOME INFINITE PRODUCTS J MC LAUGHLIN Abstract Let fx Z[x] Set f 0x = x and for n 1 define f nx = ff n 1x We describe several infinite
More informationSEVERAL PRODUCTS OF DISTRIBUTIONS ON MANIFOLDS
Novi Sad J. Math. Vol. 39, No., 2009, 3-46 SEVERAL PRODUCTS OF DISTRIBUTIONS ON MANIFOLDS C. K. Li Abstract. The problem of defining products of distributions on manifolds, particularly un the ones of
More informationOpen problems. Christian Berg a a Department of Mathematical Sciences, University of. Copenhagen, Copenhagen, Denmark Published online: 07 Nov 2014.
This article was downloaded by: [Copenhagen University Library] On: 4 November 24, At: :7 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 72954 Registered office:
More informationInequalities for logarithmic and exponential functions
General Mathematics Vol. 12, No. 2 (24), 47 52 Inequalities for logarithmic and exponential functions Eugen Constantinescu Abstract In [1] it was announced the precise inequality : for x 1 4 ( ) ( ) 12x
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Title Composition operators on Hilbert spaces of entire functions Author(s) Doan, Minh Luan; Khoi, Le Hai Citation
More informationSome Explicit Biorthogonal Polynomials
Some Explicit Biorthogonal Polynomials D. S. Lubinsky and H. Stahl Abstract. Let > and ψ x) x. Let S n, be a polynomial of degree n determined by the biorthogonality conditions Z S n,ψ,,,..., n. We explicitly
More informationBounds on Turán determinants
Bounds on Turán determinants Christian Berg Ryszard Szwarc August 6, 008 Abstract Let µ denote a symmetric probability measure on [ 1, 1] and let (p n ) be the corresponding orthogonal polynomials normalized
More informationChange of variable formulæ for regularizing slowly decaying and oscillatory Cauchy and Hilbert transforms
Change of variable formulæ for regularizing slowly decaying and oscillatory Cauchy and Hilbert transforms Sheehan Olver November 11, 2013 Abstract Formulæ are derived for expressing Cauchy and Hilbert
More informationNOTES Edited by William Adkins. On Goldbach s Conjecture for Integer Polynomials
NOTES Edited by William Adkins On Goldbach s Conjecture for Integer Polynomials Filip Saidak 1. INTRODUCTION. We give a short proof of the fact that every monic polynomial f (x) in Z[x] can be written
More informationTHE MULTIPLICATIVE ERGODIC THEOREM OF OSELEDETS
THE MULTIPLICATIVE ERGODIC THEOREM OF OSELEDETS. STATEMENT Let (X, µ, A) be a probability space, and let T : X X be an ergodic measure-preserving transformation. Given a measurable map A : X GL(d, R),
More informationMADHAVA MATHEMATICS COMPETITION, December 2015 Solutions and Scheme of Marking
MADHAVA MATHEMATICS COMPETITION, December 05 Solutions and Scheme of Marking NB: Part I carries 0 marks, Part II carries 30 marks and Part III carries 50 marks Part I NB Each question in Part I carries
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationTHE DENOMINATORS OF POWER SUMS OF ARITHMETIC PROGRESSIONS. Bernd C. Kellner Göppert Weg 5, Göttingen, Germany
#A95 INTEGERS 18 (2018) THE DENOMINATORS OF POWER SUMS OF ARITHMETIC PROGRESSIONS Bernd C. Kellner Göppert Weg 5, 37077 Göttingen, Germany b@bernoulli.org Jonathan Sondow 209 West 97th Street, New Yor,
More information1 Boas, problem p.564,
Physics 6C Solutions to Homewor Set # Fall 0 Boas, problem p.564,.- Solve the following differential equations by series and by another elementary method and chec that the results agree: xy = xy +y ()
More informationA LINEAR BINOMIAL RECURRENCE AND THE BELL NUMBERS AND POLYNOMIALS
Applicable Analysis and Discrete Mathematics, 1 (27, 371 385. Available electronically at http://pefmath.etf.bg.ac.yu A LINEAR BINOMIAL RECURRENCE AND THE BELL NUMBERS AND POLYNOMIALS H. W. Gould, Jocelyn
More information(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.
1. Complex numbers A complex number z is defined as an ordered pair z = (x, y), where x and y are a pair of real numbers. In usual notation, we write z = x + iy, where i is a symbol. The operations of
More informationAn Application of Catalan Numbers on Cayley Tree of Order 2: Single Polygon Counting
BULLETIN of the Malaysian Mathematical Sciences Society http://math.usm.my/bulletin Bull. Malays. Math. Sci. Soc. (2) 31(2) (2008), 175 183 An Application of Catalan Numbers on Cayley Tree of Order 2:
More informationGeneralized Akiyama-Tanigawa Algorithm for Hypersums of Powers of Integers
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 16 (2013, Article 1332 Generalized Aiyama-Tanigawa Algorithm for Hypersums of Powers of Integers José Luis Cereceda Distrito Telefónica, Edificio Este
More informationSums of Products of Bernoulli Numbers*
journal of number theory 60, 234 (996) article no. 00 Sums of Products of Bernoulli Numbers* Karl Dilcher Department of Mathematics, Statistics, and Computing Science, Dalhousie University, Halifax, Nova
More informationHERMITE-PADÉ APPROXIMANTS FOR A PAIR OF CAUCHY TRANSFORMS WITH OVERLAPPING SYMMETRIC SUPPORTS
This is the author's manuscript of the article published in final edited form as: Aptekarev, A. I., Van Assche, W. and Yattselev, M. L. 07, Hermite-Padé Approximants for a Pair of Cauchy Transforms with
More informationPainlevé equations and orthogonal polynomials
KU Leuven, Belgium Kapaev workshop, Ann Arbor MI, 28 August 2017 Contents Painlevé equations (discrete and continuous) appear at various places in the theory of orthogonal polynomials: Discrete Painlevé
More informationP AC COMMUTATORS AND THE R TRANSFORM
Communications on Stochastic Analysis Vol. 3, No. 1 (2009) 15-31 Serials Publications www.serialspublications.com P AC COMMUTATORS AND THE R TRANSFORM AUREL I. STAN Abstract. We develop an algorithmic
More informationSOME VARIANTS OF LAGRANGE S FOUR SQUARES THEOREM
Acta Arith. 183(018), no. 4, 339 36. SOME VARIANTS OF LAGRANGE S FOUR SQUARES THEOREM YU-CHEN SUN AND ZHI-WEI SUN Abstract. Lagrange s four squares theorem is a classical theorem in number theory. Recently,
More informationq-series Michael Gri for the partition function, he developed the basic idea of the q-exponential. From
q-series Michael Gri th History and q-integers The idea of q-series has existed since at least Euler. In constructing the generating function for the partition function, he developed the basic idea of
More informationON LINEAR COMBINATIONS OF
Física Teórica, Julio Abad, 1 7 (2008) ON LINEAR COMBINATIONS OF ORTHOGONAL POLYNOMIALS Manuel Alfaro, Ana Peña and María Luisa Rezola Departamento de Matemáticas and IUMA, Facultad de Ciencias, Universidad
More informationLECTURE 5, FRIDAY
LECTURE 5, FRIDAY 20.02.04 FRANZ LEMMERMEYER Before we start with the arithmetic of elliptic curves, let us talk a little bit about multiplicities, tangents, and singular points. 1. Tangents How do we
More informationTheorem 5.3. Let E/F, E = F (u), be a simple field extension. Then u is algebraic if and only if E/F is finite. In this case, [E : F ] = deg f u.
5. Fields 5.1. Field extensions. Let F E be a subfield of the field E. We also describe this situation by saying that E is an extension field of F, and we write E/F to express this fact. If E/F is a field
More informationDecomposition of a recursive family of polynomials
Decomposition of a recursive family of polynomials Andrej Dujella and Ivica Gusić Abstract We describe decomposition of polynomials f n := f n,b,a defined by f 0 := B, f 1 (x := x, f n+1 (x = xf n (x af
More informationON `-TH ORDER GAP BALANCING NUMBERS
#A56 INTEGERS 18 (018) ON `-TH ORDER GAP BALANCING NUMBERS S. S. Rout Institute of Mathematics and Applications, Bhubaneswar, Odisha, India lbs.sudhansu@gmail.com; sudhansu@iomaorissa.ac.in R. Thangadurai
More informationAbstract. 2. We construct several transcendental numbers.
Abstract. We prove Liouville s Theorem for the order of approximation by rationals of real algebraic numbers. 2. We construct several transcendental numbers. 3. We define Poissonian Behaviour, and study
More information