Painlevé equations and orthogonal polynomials

Size: px
Start display at page:

Download "Painlevé equations and orthogonal polynomials"

Transcription

1 KU Leuven, Belgium Kapaev workshop, Ann Arbor MI, 28 August 2017

2 Contents Painlevé equations (discrete and continuous) appear at various places in the theory of orthogonal polynomials: Discrete Painlevé equations for the recurrence coefficients of orthogonal polynomials Painlevé differential equations for the recurrence coefficients (Toda flows) Rational solutions of Painlevé equations [Special function solutions of Painlevé equations] [Local asymptotics for orthogonal polynomials at critical points] In preparation: Orthogonal Polynomials and Painlevé Equations Australian Mathematical Society Lecture Series

3 Orthogonal polynomials p n (x)p m (x)w(x) dx = δ m,n, p n (x) = γ n x n + p n (x k )p m (x k )w k = δ m,n. k=0 Three term recurrence relation xp n (x) = a n+1 p n+1 (x) + b n p n (x) + a n p n 1 (x), orthonormal xp n (x) = P n+1 (x)+b n P n (x)+anp 2 n 1 (x), monic polynomials Orthogonal polynomials on the unit circle (OPUC) 1 2π 2π 0 ϕ n (z)ϕ m (z)v(θ) dθ = δ m,n, z = e iθ, ϕ n (z) = κ n z n + zφ n (z) = Φ n+1 (z) + α n Φ n(z), Φ n(z) = z n Φ n (1/z).

4 discrete Painlevé I Let us consider w(x) = e x4 +tx 2 on (, ). The symmetry w( x) = w(x) implies b n = 0, hence The structure relation is xp n (x) = a n+1 p n+1 (x) + a n p n 1 (x). p n(x) = A n p n 1 (x) + C n p n 3 (x). Compatibility between these two relations gives ( 4an 2 an an 2 + an 1 2 t ) = n 2 This is known as discrete Painlevé I (d-p I ).

5 Some history Already in old work of Laguerre (1885), Also in work of Shohat (1939), Independently found by Freud (1976), Positive solution analyzed by Nevai (1983). Asymptotic expansion by Máté-Nevai-Zaslavsky (1985) Recognized as d-p I by Fokas-Its-Kitaev (1991). Relation with continuous Painlevé IV (Magnus, 1995)

6 Asymptotic behavior of a 2 n Theorem (Freud) The recurrence coefficients for the weight w(x) = e x4 +tx 2 satisfy lim n a n n 1/4 =

7 Special solution Put x n = a 2 n, then (for t = 0) x n (x n+1 + x n + x n 1 ) = an, a = 1/4. (1) For orthogonal polynomials we want a solution with x 0 = 0 (because a 2 0 = 0) and x n > 0 for n 1. Theorem (Lew and Quarles, Nevai) There is a unique solution of (1) for which x 0 = 0 and x n > 0 for all n 1.

8 Langmuir lattice or Kac-van Moerbeke equations General orthogonal polynomials P n (x)p m (x) dµ(x) = 0, m n. R Theorem Let µ be a symmetric positive measure on R for which all the moments exist and let µ t be the measure for which dµ t (x) = e tx2 dµ(x), where t R is such that all the moments of µ t exist. Then the recurrence coefficients of the orthogonal polynomials for µ t satisfy the differential-difference equations d dt a2 n = a 2 n(a 2 n+1 a 2 n 1), n 1.

9 Painlevé IV Put a 2 n = x n, then n = 4x n (x n+1 + x n + x n 1 t/2), (2) x n = x n (x n+1 x n 1 ) (3) Eliminate x n+1 and x n 1 using (2) (3) x n = (x n) 2 + 3x 3 ( ) n n 2x n 2 tx n 2 + x n 2 + t2 n2 8 32x n This is Painlevé IV.

10 Orthogonal polynomials on the imaginary line The weight e x4 +tx 2 behaves well on R, but it also tends to 0 along the imaginary axis ir. Orthogonal polynomials (Q n ) n on the imaginary axis +i i Q n (x)q m (x)e x4 +tx 2 dx = 0, n m, xq n (x) = Q n+1 (x) b n Q n 1 (x). One has Q n (x; t) = ( i) n P n (ix; t). Hence b n (t) = a 2 n( t). The (b n ) n is the unique negative solution of d-p I with b 0 = 0.

11 Orthogonal polynomials on the cross We can combine R and ir and look for orthogonal polynomials (R n ) n for which (n m) α R n (x)r m (x)e x4 +tx 2 dx+β +i i with α, β > 0. xr n (x) = R n+1 (x) c n R n 1 (x) R n (x)r m (x)e x4 +tx 2 dx = 0, The (c n ) n still satisfy d-p I, but with initial conditions where c 0 = 0, m 2k (t) = c 1 = αm 2(t) βm 2 ( t) αm 0 (t) + βm 0 ( t), x 2k e x4 +tx 2 dx are in terms of parabolic cylinder functions.

12 Orthogonal polynomials on the cross The solution of d-p I behaves very different. (t = 0, β/α = 1/2)

13 Orthogonal polynomials on the cross ( 4x n x n+1 + x n + x n 1 t ) = n 2 The case α = β and t = 0 is special: the initial values are x 0 = x 1 = 0. 0 is a singularity of d-p I and gives x 2 =, hence R 3 does not exist. The singularity is confined to a finite number of terms. Property For t = 0 and α = β one has that R 4n 1 does not exist for n 1. Furthermore R 4n (x) = r n (x 4 ), R 4n+1 (x) = xs n (x 4 ), R 4n+2 (x) = x 2 s n (x 4 ).

14 discrete Painlevé II Let v(θ) = e t cos θ on the unit circle Trigonometric moments: modified Bessel functions 1 2π 2π 0 e inθ v(θ) dθ = I n (t), v( θ) = v(θ) implies that α n (t) are real valued. v(θ) = ˆv(z), z = e iθ, ˆv(z) = exp ( t z + 1 z 2 ). Pearson equation ˆv (z) = t 2 ( 1 1 ) z 2 ˆv(z).

15 Structure relation Property The monic orthogonal polynomials for v(θ) = e t cos θ satisfy Φ n(z) = nφ n 1 (z) + B n Φ n 2 (z), for a sequence (B n ) n. In fact, one has B n = t 2 κ 2 n 2 κ 2. n

16 Compatibility zφ n (z) = Φ n+1 (z) + α n Φ n(x) Φ n(z) = nφ n 1 (z) + B n Φ n 2 (z) Theorem (Periwal and Shevitz) The Verblunsky coefficients for the weight v(θ) = e t cos θ satisfy with initial values t 2 (α n+1 + α n 1 ) = (n + 1)α n 1 αn 2, α 1 = 1, α 0 = I 1(t) I 0 (t). This is discrete Painlevé II (d-p II )

17 Discrete Painlevé II Let x n = α n 1, then x n+1 + x n 1 = αnx n 1 x 2 n (4) We need a solution with x 0 = 1 and x n < 1 for n 1. Such a solution is unique. Theorem Suppose α > 0. Then there is a unique solution of (4) for which x 0 = 1 and 1 < x n < 1. The solution corresponds to x 1 = I 1 (2α)/I 0 (2/α) and is positive for every n 0.

18 Asymptotic behavior This solution converges to zero (fast) x n 0. Property The solution of d-p II with x 0 = 1 and 0 < x n < 1 for n 1 satisfies 1 α n n! x n n 4n n! k=0 α k α n (2n)! 1 ( e ) n. k! 2 αn

19 The Ablowitz-Ladik lattice The lattice equations corresponding to orthogonal polynomials on the unit circle are Ablowitz-Ladik lattice equations (or Schur flow). Theorem Let ν be a positive measure on the unit circle which is symmetric (the Verblunsky coefficients are real). Let ν t be the modified measure dν t (θ) = e t cos θ dν(θ), with t R. The Verblunsky coefficients (α n (t)) n for the measure ν t then satisfy 2α n = (1 α 2 n)(α n+1 α n 1 ), n 0.

20 Painlevé V and III d-p II gives and Ablowitz-Ladik gives α n+1 + α n 1 = 2nα n t(1 α 2 n) α n+1 α n 1 = 2α n 1 αn 2. Eliminate α n+1 and α n 1 to find α α n n = 1 αn 2 (α n) 2 α n t α n (1 α 2 n) + (n + 1)2 t 2 α n 1 αn 2. If we put α n = 1+y 1 y, then y satisfies Painlevé V with γ = 0. The ratio w n = α n /α n 1 satisfies Painlevé III [Hisakado, Tracy-Widom]

21 Generalized Charlier polynomials Discrete orthogonal polynomials c k P n (k)p m (k) (β) k k! k=0 = 0, n m. b n + b n 1 n + β = cn an 2, (an+1 2 c)(an 2 c) = c(b n n)(b n n + β 1). This corresponds to a limiting case of discrete Painlevé with surface/symmetry D (1) 4 in Sakai s classification.

22 Toda Lattice Put c = c 0 e t, then c k (β) k k! = etk c k 0 (β) k k! Theorem Suppose µ is a positive measure on the real line and dµ t = e tx dµ(x), where t is such that all the moments of µ t exist. Then the recurrence coefficients b n (t) and a 2 n(t) of the orthogonal polynomials for the measure µ t satisfy d dt a2 n = an(b 2 n b n 1 ), n 1, d dt b n = an+1 2 an, 2 n 0.

23 generalized Charlier polynomials Put x n = a 2 n and y n = b n and x n = dx n /da, y n = dy n /da, then (x n a)(x n+1 a) = a(y n n)(y n n + β 1), y n + y n 1 n + β = an x n and the Toda lattice equations are ax n = x n (y n y n 1 ), ay n = x n+1 x n. Eliminate y n 1 and x n+1, and put x n = y = 1 2 ( 1 2y + 1 ) (y ) 2 y ( y)2 n 2 y +(1 y 1 a a 2 2 This is Painlevé V with δ = 0. This can be transformed to Painlevé III. a 1 y, then y(a) satisfies ) (β 1)2 2y 2y a

24 generalized Meixner polynomials Discrete orthogonal polynomials P n (k)p m (k) (γ) ka k = 0, n m. (β) k k! k=0 Put an 2 = na (γ 1)u n, and b n = n + γ β + a γ 1 a v n, then (u n + v n )(u n+1 + v n ) = γ 1 ( a 2 v n (v n a) v n a γ β ), γ 1 ( u n (u n + v n )(u n + v n 1 ) = u n an (u n + a) u n + a γ β ). γ 1 γ 1 Initial values a 0 = 0, b 0 = γa M(γ + 1, β + 1, a) β M(γ, β, a) This is asymmetric discrete Painlevé IV or d-p(e (1) 6 /A(1) 2 ).

25 generalized Meixner If we put v n (a) = ( ) a ay (1 + β 2γ)y 2 + (n + 1 a + β 2γ)y n, 2(γ 1)(y 1)y then y = with ( 1 2y + 1 ) (y ) 2 y y 1 a A = +(y 1)2 a 2 ( Ay+ B y ) Cy + 1) + +Dy(y a y 1 (β 1)2, B = n2 2 2, C = n β + 2γ, D = 1 2, which is Painlevé V.

26 Other examples Chen and Its (2010): w(x) = x α e x e t/x on [0, ) Put b n = 2n + α c n, an 2 = n(n + α) + y n + n 1 j=0 c j, and c n = 1/x n d-p((2a 1 ) (1) /D (1) 6 ) x n + x n 1 = nt (2n + α)y n y n (y n t) y n + y n+1 = t 2n + α + 1 x n 1 x 2 n c n = (c n) 2 c n c n t + (2n + α + 1)c2 n t 2 + c3 n t 2 + α t 1 c n which is Painlevé III.

27 Other examples Basor, Chen, Ehrhardt (2010): w(x) = (1 x) α (1 + x) β e tx tb n = 2n α + β t 2R n, t(t + R n )a 2 n = n(n + β) (2n + α + β)r n tr n(r n + α) R n 2t(r n+1 + r n ) = 4Rn 2 2R n (2n α + β t) 2αt, ( t 2 n(n + β) (2n + α + β)r n = r n (r n + α) + t + t ) R n R n 1 R n and for y = 1 + t/r n y = 3y 1 2y(y 1) (y ) 2 y t +2(2n+1+α+β)y t which is Painlevé V. R n 1 ( + 1) (y 1)2 α 2 ) y 2y(y + y 1 t 2 2 β2 2y

28 Other examples q-orthogonal polynomials: w(x) = x α ( x 2 ; q 2 ) ( q 2 /x 2 ; q 2 ), x [0, ) gives rise to q-discrete Painlevé III x n 1 x n+1 = (x n + q α ) 2 (q n+α x n + 1) 2. x α ( p/x 2 ; q 2 ) w(x) = ( x 2 ; q 2 ) ( q 2 /x 2 ; q 2, x [0, ) ) gives rise to q-discrete Painlevé V (z n z n 1 1)(z n z n+1 1) = (z n + q 2 α /p) 2 (z n pq α 2 ) 2 (q n+α/2 1 pz n + 1) 2. w(x) = x α (q 2 x 2 ; q 2 ), x {q k, k = 0, 1, 2, 3,..., } gives again q-discrete Painlevé V

29 Other examples Bi-orthogonal polynomials on the unit circle (Forrester and Witte, 2006) { w(z) = z µ ω (1+z) 2ω 1 (1+tz) 2µ 1, θ / (π φ, π), 1 ξ, θ (π φ, π)., t = eiφ, gives rise to g n+1 g n = t (f n + n)(f n + n + 2µ), f n (f n 2ω 1 ) f n + f n 1 = 2ω 1 + n 1 + µ + ω + g n 1 discrete Painlevé V (Sakai s surface D (1) 4 ). (n + µ + ω)t g n t

30 Other examples Biane (2014) worked out a q-version w(e iθ ) = (ae iθ ; q) (be iθ ; q) and found discrete Painlevé equations corresponding to A (1) 3 r n+2 (a bq n+2 )+qr n (a bq n ) = 2r ( ) n+1 1 rn+1 2 (1 q)(a+bq n+1 )β n+1 +(1 q n+1 )(ab+q), 2, with n β n = r r r k 1, r n = Φ n (0). k=1

31 Other examples Witte (2015) worked out the most general case (deformed Askey-Wilson polynomials): Pearson-type equation for the weight D x w(x) = 2V (x) W (x) M xw(x), where D x is a divided difference operator on a special nonuniform lattice obtained from a quadratic equation Ay 2 + 2Bxy + Cx 2 + 2Dy + 2Ex + F = 0, and M x is an average operator for the forward and backward moves on the lattice, and V, W are polynomials. W (z) ± yv (z) = z 3 6 (1 a j q 1/2 z ±1 ), gives the q-discrete Painlevé cases for E (1) 7 in Sakai s scheme. j=1

32 Rational solutions of Painlevé equations

33 Rational solutions of Painlevé II y = 2y 3 + xy + α has rational solutions if and only if α = n Z. The solutions are y n = d dz log Q n 1(z) Q n (z) where Q n are the Yablonskii-Vorobiev polynomials p 1 p 3 p 5 p 2n 1 p 1 p 3 p 5 p 2n 1 Q n (z) = c n det.... p (n 1) 1 p (n 1) 3 p (n 1) 5 p (n 1) 2n 1 with polynomial p k given by p k (z)λ k = exp(zλ 4 3 λ3 ). k=0

34 Rational solutions of Painlevé II The polynomials p k are multiple orthogonal polynomials on the 3-star 2 Γ = {z C : arg(z) = 2πk 3 }, k=0 p n (z)ai(2 2/3 z )z k dz = 0, 0 k n 2 1, Γ Γ p n (z) z z Ai (2 2/3 z )z k dz = 0, 0 k n 2 1.

35 Rational solutions of Painlevé II Bertola and Bothner (2015) gave another determinant representation µ 0 (z) µ 1 (z) µ 2 (z) µ n 1 (z) Qn 1(z) 2 µ 1 (z) µ 2 (z) µ 3 (z) µ n (z) = d n det.... µ n 1 (z) µ n (z) µ n+1 (z)... µ 2n 2 (z) with µ n (z) = p n (2 2/3 z)2 2n/3. Advantage: this is a Hankel matrix, and hence allows to use the moment problem and orthogonal polynomials µ 0 (z) µ 1 (z) µ 2 (z) µ n (z) d µ 1 (z) µ 2 (z) µ 3 (z) µ n+1 (z) n P n (x) = Qn 1 2 det (z).... µ n 1 (z) µ n (z) µ n+1 (z)... µ 2n 1 (z) 1 x x 2 x n

36 Rational solutions of Painlevé II z is a zero of Q n 1 (and a pole of the rational solution) if and only if the monic orthogonal polynomial P n (x) does not exist. This allows to locate the zeros of Q n 1 by investigating the Riemann-Hilbert problem for the orthogonal polynomials (P k ) k N. Theorem (Bertola-Bothner, Buckingham-Miller) The zeros of Q n (n 2/3 z) accumulate on a triangular shaped region with corners x k = 3/ 3 2e 2πi 3 k (k = 0, 1, 2). The sides of the region are given by ( R 2 log a 3 ia + ) 1 + 2a 3 4a3 1 2a 3a 3 = 0, where a = a(x) is a solution of the cubic equation 1 + 2xa 2 4a 3 = 0.

37

38 Rational solutions of Painlevé III y = (y ) 2 y y z + αy 2 + β + y 3 1 z y Has rational solutions if and only if α + β = 4n or α β = 4n for n Z. For α = 2n + 2µ 1 and β = 2n 2µ 1 y n = 1 + d dz S n 1 (z; µ 1). S n (z; µ) and for α = 2n + 2µ 1 and β = 2n 2µ + 1 y n = 1 d dz S n 1 (z; µ) S n (z; µ 1).

39 Rational solutions of Painlevé III S n (z; µ) are Umemura polynomials S n (z; µ) = c n det L µ n n ( z) L µ n 1 n+1 ( z) L µ 2n+1 2n 1 n 2 ( z) L µ n+1 n 1 ( z) L µ 2n+3 L µ n+2. L µ+n 2 n+2 ( z) with L α n (x) the Laguerre polynomials. Lµ+n 3 n+3 2n 3 ( z) ( z).. ( z) Lµ 1 1 ( z)

40 Rational solutions of Painlevé III Wronskian formula: L µ 1 1 ( z) L µ 3 3 ( z) L µ 2n+1 2n 1 ( z) [L µ 1 1 ( z)] [L µ 3 3 ( z)] [L µ 2n+1 2n 1 ( z)] S n (z; µ) = c n det... [L µ 1 1 ( z)] (n 1) [L µ 3 3 ( z)] (n 1) [L µ 2n+1 2n 1 with L α n (x) the Laguerre polynomials. ( z)] (n 1)

41 Zeros of Umemura polynomial S 10 (µ), 9 µ 9

42 Rational solutions of Painlevé IV y = (y ) 2 2y Has rational solutions if and only if d dz log H m+1,n(z) H m,n, y 3 + 4zy 2 + 2(z 2 α)y + β y α = m, β = 2(2n m + 1) 2 α = m, β = 2(2n m )2. d dz log H m,n(z) H m,n+1 (z), with H m,n the generalized Hermite polynomials 2z + d d log H m,n+1(z) H m+1,n (z) 2 3 z+ d dz log Q m+1,n(z) Q m,n (z), 2 3 z+ d dz log Q m,n(z) Q m,n+1 (z), 2 3 z+ d dz log Q m,n+1(z) Q m+1,n (z) with Q m,n the generalized Okamoto polynomials

43 Rational solutions of Painlevé IV The generalized Hermite polynomials are H m (z) H m+1 (z) H m+n 1 (z) H m+1 (z) H m+2 (z) H m+n (z) H m,n (z) = det... H m+n 1 (z) H m+n (z) H m+2n 2 (z) and H m (z) are Hermite polynomials.

44 Zeros of generalized Hermite polynomials

45 Zeros of generalized Hermite polynomials Theorem (Buckingham) The zeros of the generalized Hermite polynomial H m,n (m 1/2 z) for the rational solution of Painlevé IV accumulate in a square shaped region with corners given by four solutions of the equation r 4 x 8 24r 2 (r 3 +r+1)x 4 +32r(2r 3 +3r 2 3r 2)x 2 48(r 2 +r+1) 2 = 0 (not real or imaginary), and sides given by the curves R ( (1 + r)r 1/2 ( ) x 4 R (1 + r) log 2R 2 (1 + r)q S + (r 1) log ( (1 + r)q 3 + (1 + r)q 2 R + s ) + log(s 2 4Q 2 ) ) = 0 where r = m/n 1.

46 Zeros of generalized Hermite polynomials Theorem (continued) Q = Q(x, r) is the unique solution of 3(1 + r) 2 Q 4 + 8(1 + r)r 1/2 xq 3 + 4(r 1 + rx 2 )Q 2 4 = 0 for which Q(x, r) = x + O(x 2 ) as x, S = (1 + r)q 3 + 2r 1/2 Q 2 and ( (1 + r) 2 Q 4 + 2(1 + r)qs + 4 ) 1/2 R =. (1 + r)q

47 Zeros of generalized Hermite polynomials

48 Rational solutions of Painlevé IV The Okamoto polynomials are Q n (z) = Q n+1,0 (z) = Wr(H 3n 1, H 3n 4,..., H 2 ), R n (z) = Q n,1 (z) = Wr(H 3n 2, H 3n 5,..., H 1 ) The generalized Okamoto polynomials are Q m, n (z) = H λ (z), λ = (m + 2n, m + 2n 2,..., m + 2, m, m, m 1, m 1,..., 1, 1), Q n,m+1 (z) = H λ (z), λ = (m + 2n 1, m + 2n 3,..., m + 1, m, m, m 1, m 1,..., 1, 1). H λ1 H λ1+1 H λ1+2 H λ1+k 1 H λ2 1 H λ2 H λ2+1 H λ2+k 2 H λ (z) = det H λ3 2 H λ3 1 H λ3 2 H λ3+k H λk k+1 H λk k+2 H λk k+3 H λk

49 Zeros of generalized Okamoto polynomials

50 Rational solutions of Painlevé V y = ( 1 2y + 1 ) (y ) 2 y ( 1)2 +(y y 1 z z 2 αy + β ) + γy + 1) y(y y z 2(y 1) Has rational solutions if 1 α = 1 2 (m ± 1)2 and β = n2 2, with n > 0, m + n odd, α 0 when m < n; 2 α = n2 2, β = 1 2 (m ± 1)2, with n > 0, m + n odd, β 0 when m < n; 3 α = a 2, β = 1 2 (a + n)2, and γ = m, with m + n even; 4 α = (b + n) 2, β = b2, and γ = m, with m + n even; 5 α = 1 8 (2m + 1)2, β = 1 8 (2n + 1)2. 2

51 Rational solutions of Painlevé V In terms of Umemura polynomials L r n(x) L r n+1 (x) Lr 2n 1 (x) L r n 2 U 0,n (x; µ) = C n det (x) Lr n 1 (x) Lr 2n 3 (x)..., L r 2 n (x) Lr 3 n (x) Lr 1 (x) L r m( x) L r m+1 ( x) Lr 2m 1 ( x) L r m 2 U m,0 (x; µ) = D m det ( x) Lr m 1 ( x) Lr 2m 3 ( x)..., L r 2 m ( x) Lr 3 m ( x) Lr 1 ( x) with L α n (x) the Laguerre polynomials. r = µ+n 1 r = µ m 1

52 Rational solutions of Painlevé V and generalized Umemura polynomials U m,n (x; µ) L r 1 ( x) Lr 0 ( x) Lr m+2 ( x) Lr m+1 ( x) Lr m n+2 ( x) L r 3 ( x) Lr 2 ( x) Lr m+4 ( x) Lr m+3 ( x) Lr m n+4 ( x)..... L r det 2m 1 ( x) Lr 2m 2 ( x) Lr m( x) L r m 1 ( x) Lr m n ( x) L r n m (x) Lr n m+1 (x) Lr n 1 (x) Lr n(x) L r 2n 1 (x) L r n m 2 (x) Lr n m 1 (x) Lr n 3 (x) Lr n 2 (x) Lr 2n 3 (x)..... L r n m+2 (x) Lr n m+3 (x) Lr n+1 (x) Lr n+2 (x) Lr 1 (x) where r = µ m + n 1.

53 Zeros of U 8,10 (µ = 2k + 1), 11 k 9

54 Rational solutions of Painlevé VI Special function solution in terms of terminating hypergeometric functions These are determinants with Jacobi polynomials Little research has been done on these rational solutions.

55 Special function solutions of Painlevé equations

56 Special functions Painlevé I Painlevé II Airy functions Painlevé III Bessel functions Painlevé IV Parabolic cylinder functions Painlevé V confluent hypergeometric functions Painlevé VI hypergeometric functions seed function: comes from a Riccati equation Bäcklund transformations give all the special function solutions in terms of the seed function.

57 Connection with orthogonal polynomials Hankel determinant with moments of a measure µ m 0 m 1 m 2 m n 1 m 1 m 2 m 3 m n n = det m 2 m 3 m 4 m n+1, m k =.... m n 1 m n m n+1 m 2n 2 m 0 = and n is a Wronskian w(x)e xt dx, m k = d k m 0 dt k x k dµ(x). n = det d i+j m 0 dt i+j.

58 Connection with orthogonal polynomials Recurrence coefficients an 2 = n+1 n 1 2, b n = n+1 n, n n+1 n where m 0 m 1 m 2 m n 2 m n m 1 m 2 m 3 m n 1 m n+1 n = det m 2 m 3 m 4 m n m n+2 = d dt n m n 1 m n m n+1 m 2n 3 m 2n 1

59 Examples w(x) = e x4 +tx 2 m 0 = 2 1/4 πe t2 /8 U(0, t/ 2) P IV w(x) = e x3 /3+tx m 0 = 2πiAi(t) P II w k = ak (β) k k! on N m 0 = a (β 1)/2 Γ(β)I β 1 (2 a) P III w(x) = x α e x e s/x m k K k 1 (2 t) P III w(x) = x a 1 (1 x) b a 1 e xt m 0 = Γ(a)Γ(b a)m(a, b, t) P V w k = (γ) ka k (β) k k! on N m 0 = M(γ, β, a) P V w(x) = x b 1 (1 x) c b 1 (1 xt) a m 0 = Γ(b)Γ(b c) 2 F 1 (a, b; c; t) P VI

60 References P.A. Clarkson, Painlevé equations nonlinear special functions, Lecture Notes in Mathematics 1883, Springer, Berlin, 2006, pp K. Kajiwara, M. Noumi, Y. Yamada, Geometric aspects of Painlevé equations, J. Phys. A: Math. Theor. 50 (2017), no. 7, (164. pp.) W. Van Assche, Discrete Painlevé equations for recurrence coefficients of orthogonal polynomials, in Difference Equations, Special Functions and Orthogonal Polynomials, World Scientific, 2007, pp W. Van Assche, Orthogonal Polynomials and Painlevé Equations, Australian Mathematical Society Lecture Series, Cambridge University Press (to appear).

Multiple orthogonal polynomials. Bessel weights

Multiple 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 information

Painlevé VI and Hankel determinants for the generalized Jacobi weight

Painlevé 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 information

Recurrence Coef f icients of a New Generalization of the Meixner Polynomials

Recurrence 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 information

Special solutions of the third and fifth Painlevé equations and vortex solutions of the complex Sine-Gordon equations

Special solutions of the third and fifth Painlevé equations and vortex solutions of the complex Sine-Gordon equations Special solutions of the third and fifth Painlevé equations and vortex solutions of the complex Sine-Gordon equations Peter A Clarkson School of Mathematics, Statistics and Actuarial Science University

More information

Markov operators, classical orthogonal polynomial ensembles, and random matrices

Markov operators, classical orthogonal polynomial ensembles, and random matrices Markov operators, classical orthogonal polynomial ensembles, and random matrices M. Ledoux, Institut de Mathématiques de Toulouse, France 5ecm Amsterdam, July 2008 recent study of random matrix and random

More information

Fredholm determinant with the confluent hypergeometric kernel

Fredholm determinant with the confluent hypergeometric kernel Fredholm determinant with the confluent hypergeometric kernel J. Vasylevska joint work with I. Krasovsky Brunel University Dec 19, 2008 Problem Statement Let K s be the integral operator acting on L 2

More information

The Density Matrix for the Ground State of 1-d Impenetrable Bosons in a Harmonic Trap

The Density Matrix for the Ground State of 1-d Impenetrable Bosons in a Harmonic Trap The Density Matrix for the Ground State of 1-d Impenetrable Bosons in a Harmonic Trap Institute of Fundamental Sciences Massey University New Zealand 29 August 2017 A. A. Kapaev Memorial Workshop Michigan

More information

arxiv:nlin/ v1 [nlin.si] 29 May 2002

arxiv:nlin/ v1 [nlin.si] 29 May 2002 arxiv:nlin/0205063v1 [nlinsi] 29 May 2002 On a q-difference Painlevé III Equation: II Rational Solutions Kenji KAJIWARA Graduate School of Mathematics, Kyushu University 6-10-1 Hakozaki, Higashi-ku, Fukuoka

More information

This ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0

This ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0 Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x

More information

OPSF, Random Matrices and Riemann-Hilbert problems

OPSF, 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 information

Recurrence coefficients of generalized Meixner polynomials and Painlevé equations

Recurrence coefficients of generalized Meixner polynomials and Painlevé equations Recurrence coefficients of generalized Meixner polynomials and Painlevé equations Lies Boelen, Galina Filipuk, Walter Van Assche October 0, 010 Abstract We consider a semi-classical version of the Meixner

More information

Additional material: Linear Differential Equations

Additional material: Linear Differential Equations Chapter 5 Additional material: Linear Differential Equations 5.1 Introduction The material in this chapter is not formally part of the LTCC course. It is included for completeness as it contains proofs

More information

Nonlinear Integral Equation Formulation of Orthogonal Polynomials

Nonlinear Integral Equation Formulation of Orthogonal Polynomials Nonlinear Integral Equation Formulation of Orthogonal Polynomials Eli Ben-Naim Theory Division, Los Alamos National Laboratory with: Carl Bender (Washington University, St. Louis) C.M. Bender and E. Ben-Naim,

More information

Difference Equations for Multiple Charlier and Meixner Polynomials 1

Difference 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 information

Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium

Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium SEA 06@MIT, Workshop on Stochastic Eigen-Analysis and its Applications, MIT, Cambridge,

More information

Painlevé Transcendents and the Hankel Determinants Generated by a Discontinuous Gaussian Weight

Painlevé Transcendents and the Hankel Determinants Generated by a Discontinuous Gaussian Weight Painlevé Transcendents and the Hankel Determinants Generated by a Discontinuous Gaussian Weight arxiv:1803.10085v [math-ph] May 018 Chao Min and Yang Chen November 9, 018 Abstract This paper studies the

More information

QUADRATIC EQUATIONS AND MONODROMY EVOLVING DEFORMATIONS

QUADRATIC EQUATIONS AND MONODROMY EVOLVING DEFORMATIONS QUADRATIC EQUATIONS AND MONODROMY EVOLVING DEFORMATIONS YOUSUKE OHYAMA 1. Introduction In this paper we study a special class of monodromy evolving deformations (MED). Chakravarty and Ablowitz [4] showed

More information

Spectral Theory of Orthogonal Polynomials

Spectral Theory of Orthogonal Polynomials Spectral Theory of Orthogonal Polynomials Barry Simon IBM Professor of Mathematics and Theoretical Physics California Institute of Technology Pasadena, CA, U.S.A. Lecture 1: Introduction and Overview Spectral

More information

Jacobi-Angelesco multiple orthogonal polynomials on an r-star

Jacobi-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 information

arxiv:solv-int/ v1 20 Oct 1993

arxiv:solv-int/ v1 20 Oct 1993 Casorati Determinant Solutions for the Discrete Painlevé-II Equation Kenji Kajiwara, Yasuhiro Ohta, Junkichi Satsuma, Basil Grammaticos and Alfred Ramani arxiv:solv-int/9310002v1 20 Oct 1993 Department

More information

PHYS 404 Lecture 1: Legendre Functions

PHYS 404 Lecture 1: Legendre Functions PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions

More information

Series solutions of second order linear differential equations

Series solutions of second order linear differential equations Series solutions of second order linear differential equations We start with Definition 1. A function f of a complex variable z is called analytic at z = z 0 if there exists a convergent Taylor series

More information

Approximation theory

Approximation theory Approximation theory Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 1 1.3 6 8.8 2 3.5 7 10.1 Least 3squares 4.2

More information

SYMMETRIES IN THE FOURTH PAINLEVÉ EQUATION AND OKAMOTO POLYNOMIALS

SYMMETRIES IN THE FOURTH PAINLEVÉ EQUATION AND OKAMOTO POLYNOMIALS M. Noumi and Y. Yamada Nagoya Math. J. Vol. 153 (1999), 53 86 SYMMETRIES IN THE FOURTH PAINLEVÉ EQUATION AND OKAMOTO POLYNOMIALS MASATOSHI NOUMI and YASUHIKO YAMADA Abstract. The fourth Painlevé equation

More information

Bessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics

Bessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy

More information

Zeros of Polynomials: Beware of Predictions from Plots

Zeros 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 information

Free Meixner distributions and random matrices

Free Meixner distributions and random matrices Free Meixner distributions and random matrices Michael Anshelevich July 13, 2006 Some common distributions first... 1 Gaussian Negative binomial Gamma Pascal chi-square geometric exponential 1 2πt e x2

More information

Simultaneous Gaussian quadrature for Angelesco systems

Simultaneous 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 information

ORTHOGONAL POLYNOMIALS WITH EXPONENTIALLY DECAYING RECURSION COEFFICIENTS

ORTHOGONAL POLYNOMIALS WITH EXPONENTIALLY DECAYING RECURSION COEFFICIENTS ORTHOGONAL POLYNOMIALS WITH EXPONENTIALLY DECAYING RECURSION COEFFICIENTS BARRY SIMON* Dedicated to S. Molchanov on his 65th birthday Abstract. We review recent results on necessary and sufficient conditions

More information

Exponential tail inequalities for eigenvalues of random matrices

Exponential tail inequalities for eigenvalues of random matrices Exponential tail inequalities for eigenvalues of random matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify

More information

Central Limit Theorems for linear statistics for Biorthogonal Ensembles

Central Limit Theorems for linear statistics for Biorthogonal Ensembles Central Limit Theorems for linear statistics for Biorthogonal Ensembles Maurice Duits, Stockholm University Based on joint work with Jonathan Breuer (HUJI) Princeton, April 2, 2014 M. Duits (SU) CLT s

More information

Positivity of Turán determinants for orthogonal polynomials

Positivity 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 information

Orthogonal Polynomials, Perturbed Hankel Determinants. and. Random Matrix Models

Orthogonal Polynomials, Perturbed Hankel Determinants. and. Random Matrix Models Orthogonal Polynomials, Perturbed Hankel Determinants and Random Matrix Models A thesis presented for the degree of Doctor of Philosophy of Imperial College London and the Diploma of Imperial College by

More information

Triangular matrices and biorthogonal ensembles

Triangular matrices and biorthogonal ensembles /26 Triangular matrices and biorthogonal ensembles Dimitris Cheliotis Department of Mathematics University of Athens UK Easter Probability Meeting April 8, 206 2/26 Special densities on R n Example. n

More information

ODE Homework Series Solutions Near an Ordinary Point, Part I 1. Seek power series solution of the equation. n(n 1)a n x n 2 = n=0

ODE Homework Series Solutions Near an Ordinary Point, Part I 1. Seek power series solution of the equation. n(n 1)a n x n 2 = n=0 ODE Homework 6 5.2. Series Solutions Near an Ordinary Point, Part I 1. Seek power series solution of the equation y + k 2 x 2 y = 0, k a constant about the the point x 0 = 0. Find the recurrence relation;

More information

Multiple Orthogonal Polynomials

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 information

Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv: v2 [math.ca] 8 Nov 2013

Multiple 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 information

Alphonse Magnus, Université Catholique de Louvain.

Alphonse Magnus, Université Catholique de Louvain. Rational interpolation to solutions of Riccati difference equations on elliptic lattices. Alphonse Magnus, Université Catholique de Louvain. http://www.math.ucl.ac.be/membres/magnus/ Elliptic Riccati,

More information

Plancherel Rotach Asymptotic Expansion for Some Polynomials from Indeterminate Moment Problems

Plancherel Rotach Asymptotic Expansion for Some Polynomials from Indeterminate Moment Problems Constr Approx : DOI.7/s5--95- Plancherel Rotach Asymptotic Expansion for Some Polynomials from Indeterminate Moment Problems Dan Dai Mourad E.H. Ismail Xiang-Sheng Wang Received: 7 February / Accepted:

More information

Entropy of Hermite polynomials with application to the harmonic oscillator

Entropy 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 information

ORTHOGONAL POLYNOMIALS FOR THE OSCILLATORY-GEGENBAUER WEIGHT. Gradimir V. Milovanović, Aleksandar S. Cvetković, and Zvezdan M.

ORTHOGONAL POLYNOMIALS FOR THE OSCILLATORY-GEGENBAUER WEIGHT. Gradimir V. Milovanović, Aleksandar S. Cvetković, and Zvezdan M. PUBLICATIONS DE L INSTITUT MATHÉMATIQUE Nouvelle série, tome 8498 2008, 49 60 DOI: 02298/PIM0898049M ORTHOGONAL POLYNOMIALS FOR THE OSCILLATORY-GEGENBAUER WEIGHT Gradimir V Milovanović, Aleksandar S Cvetković,

More information

A trigonometric orthogonality with respect to a nonnegative Borel measure

A trigonometric orthogonality with respect to a nonnegative Borel measure Filomat 6:4 01), 689 696 DOI 10.98/FIL104689M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat A trigonometric orthogonality with

More information

Multiple Orthogonal Polynomials

Multiple 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 information

Prime Numbers and Irrational Numbers

Prime Numbers and Irrational Numbers Chapter 4 Prime Numbers and Irrational Numbers Abstract The question of the existence of prime numbers in intervals is treated using the approximation of cardinal of the primes π(x) given by Lagrange.

More information

Introduction to orthogonal polynomials. Michael Anshelevich

Introduction to orthogonal polynomials. Michael Anshelevich Introduction to orthogonal polynomials Michael Anshelevich November 6, 2003 µ = probability measure on R with finite moments m n (µ) = R xn dµ(x)

More information

MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS

MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS ABDUL HASSEN AND HIEU D. NGUYEN Abstract. This paper investigates a generalization the classical Hurwitz zeta function. It is shown that many of the properties

More information

arxiv: v2 [nlin.si] 3 Feb 2016

arxiv: v2 [nlin.si] 3 Feb 2016 On Airy Solutions of the Second Painlevé Equation Peter A. Clarkson School of Mathematics, Statistics and Actuarial Science, University of Kent, Canterbury, CT2 7NF, UK P.A.Clarkson@kent.ac.uk October

More information

Journal of Computational and Applied Mathematics

Journal of Computational and Applied Mathematics Journal of Computational and Applied Mathematics 233 2010) 1554 1561 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: wwwelseviercom/locate/cam

More information

Double contour integral formulas for the sum of GUE and one matrix model

Double contour integral formulas for the sum of GUE and one matrix model Double contour integral formulas for the sum of GUE and one matrix model Based on arxiv:1608.05870 with Tom Claeys, Arno Kuijlaars, and Karl Liechty Dong Wang National University of Singapore Workshop

More information

Numerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018

Numerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018 Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start

More information

Painlevé IV and degenerate Gaussian Unitary Ensembles

Painlevé IV and degenerate Gaussian Unitary Ensembles Painlevé IV and degenerate Gaussian Unitary Ensembles arxiv:math-ph/0606064v1 8 Jun 006 Yang Chen Department of Mathematics Imperial College London 180 Queen s Gates London SW7 BZ UK M. V. Feigin Department

More information

On the singularities of non-linear ODEs

On the singularities of non-linear ODEs On the singularities of non-linear ODEs Galina Filipuk Institute of Mathematics University of Warsaw G.Filipuk@mimuw.edu.pl Collaborators: R. Halburd (London), R. Vidunas (Tokyo), R. Kycia (Kraków) 1 Plan

More information

CALCULUS JIA-MING (FRANK) LIOU

CALCULUS JIA-MING (FRANK) LIOU CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion

More information

Lecture 4.6: Some special orthogonal functions

Lecture 4.6: Some special orthogonal functions Lecture 4.6: Some special orthogonal functions Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4340, Advanced Engineering Mathematics

More information

A Numerical Approach To The Steepest Descent Method

A Numerical Approach To The Steepest Descent Method A Numerical Approach To The Steepest Descent Method K.U. Leuven, Division of Numerical and Applied Mathematics Isaac Newton Institute, Cambridge, Feb 15, 27 Outline 1 One-dimensional oscillatory integrals

More information

Orthogonal polynomials with respect to generalized Jacobi measures. Tivadar Danka

Orthogonal polynomials with respect to generalized Jacobi measures. Tivadar Danka Orthogonal polynomials with respect to generalized Jacobi measures Tivadar Danka A thesis submitted for the degree of Doctor of Philosophy Supervisor: Vilmos Totik Doctoral School in Mathematics and Computer

More information

Elliptic functions. MATH 681: Lecture Notes. Part I. 1 Evaluating ds cos s. 18 June 2014

Elliptic functions. MATH 681: Lecture Notes. Part I. 1 Evaluating ds cos s. 18 June 2014 MATH 68: Lecture Notes 8 June 4 Part I Elliptic functions Evaluating cos s Consider cos s First we can substitute u cos s and arrive at u u p3 u, where p 3 is a third-degree polynomial We can also write

More information

Local and global finite branching of solutions of ODEs in the complex plane

Local and global finite branching of solutions of ODEs in the complex plane Local and global finite branching of solutions of ODEs in the complex plane Workshop on Singularities and Nonlinear ODEs Thomas Kecker University College London / University of Portsmouth Warszawa, 7 9.11.2014

More information

Problem Set 8 Mar 5, 2004 Due Mar 10, 2004 ACM 95b/100b 3pm at Firestone 303 E. Sterl Phinney (2 pts) Include grading section number

Problem Set 8 Mar 5, 2004 Due Mar 10, 2004 ACM 95b/100b 3pm at Firestone 303 E. Sterl Phinney (2 pts) Include grading section number Problem Set 8 Mar 5, 24 Due Mar 1, 24 ACM 95b/1b 3pm at Firestone 33 E. Sterl Phinney (2 pts) Include grading section number Useful Readings: For Green s functions, see class notes and refs on PS7 (esp

More information

Power Series Solutions to the Legendre Equation

Power Series Solutions to the Legendre Equation Department of Mathematics IIT Guwahati The Legendre equation The equation (1 x 2 )y 2xy + α(α + 1)y = 0, (1) where α is any real constant, is called Legendre s equation. When α Z +, the equation has polynomial

More information

Computation of Rational Szegő-Lobatto Quadrature Formulas

Computation of Rational Szegő-Lobatto Quadrature Formulas Computation of Rational Szegő-Lobatto Quadrature Formulas Joint work with: A. Bultheel (Belgium) E. Hendriksen (The Netherlands) O. Njastad (Norway) P. GONZÁLEZ-VERA Departament of Mathematical Analysis.

More information

Differential Equations and Associators for Periods

Differential Equations and Associators for Periods Differential Equations and Associators for Periods Stephan Stieberger, MPP München Workshop on Geometry and Physics in memoriam of Ioannis Bakas November 2-25, 26 Schloß Ringberg, Tegernsee based on: St.St.,

More information

Lecture No 1 Introduction to Diffusion equations The heat equat

Lecture No 1 Introduction to Diffusion equations The heat equat Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and

More information

Measures, orthogonal polynomials, and continued fractions. Michael Anshelevich

Measures, orthogonal polynomials, and continued fractions. Michael Anshelevich Measures, orthogonal polynomials, and continued fractions Michael Anshelevich November 7, 2008 MEASURES AND ORTHOGONAL POLYNOMIALS. µ a positive measure on R. A linear functional µ[p ] = P (x) dµ(x), µ

More information

Measures, orthogonal polynomials, and continued fractions. Michael Anshelevich

Measures, orthogonal polynomials, and continued fractions. Michael Anshelevich Measures, orthogonal polynomials, and continued fractions Michael Anshelevich November 7, 2008 MEASURES AND ORTHOGONAL POLYNOMIALS. MEASURES AND ORTHOGONAL POLYNOMIALS. µ a positive measure on R. 2 MEASURES

More information

MATH 312 Section 6.2: Series Solutions about Singular Points

MATH 312 Section 6.2: Series Solutions about Singular Points MATH 312 Section 6.2: Series Solutions about Singular Points Prof. Jonathan Duncan Walla Walla University Spring Quarter, 2008 Outline 1 Classifying Singular Points 2 The Method of Frobenius 3 Conclusions

More information

Chapter 2 Formulas and Definitions:

Chapter 2 Formulas and Definitions: Chapter 2 Formulas and Definitions: (from 2.1) Definition of Polynomial Function: Let n be a nonnegative integer and let a n,a n 1,...,a 2,a 1,a 0 be real numbers with a n 0. The function given by f (x)

More information

Material for review. By Lei. May, 2011

Material for review. By Lei. May, 2011 Material for review. By Lei. May, 20 You shouldn t only use this to do the review. Read your book and do the example problems. Do the problems in Midterms and homework once again to have a review. Some

More information

An old new class of meromorphic functions. Norbert Steinmetz. June 13, 2015

An old new class of meromorphic functions. Norbert Steinmetz. June 13, 2015 Norbert Steinmetz June 13, 2015 NEW & OLD NEW: Hamiltonian system p = H q = q 2 zp α q = H p = p 2 + zq + β with Hamiltonian H(z, p, q) = 1 3 (p3 + q 3 ) + zpq + αq + βp OLD: Almost Painlevé IV 2ww = w

More information

Elementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test Find the radius of convergence of the power series

Elementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test Find the radius of convergence of the power series Elementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test 2 SOLUTIONS 1. Find the radius of convergence of the power series Show your work. x + x2 2 + x3 3 + x4 4 + + xn

More information

Painlevé VI and the Unitary Jacobi ensembles

Painlevé VI and the Unitary Jacobi ensembles Painlevé VI and the Unitary Jacobi ensembles arxiv:0911.5636v3 [math.ca] 23 Dec 2009 Yang Chen Department of Mathematics, Imperial College London, 180 Queen s Gates, London SW7 2BZ, UK ychen@imperial.ac.uk

More information

Introduction to Orthogonal Polynomials: Definition and basic properties

Introduction to Orthogonal Polynomials: Definition and basic properties Introduction to Orthogonal Polynomials: Definition and basic properties Prof. Dr. Mama Foupouagnigni African Institute for Mathematical Sciences, Limbe, Cameroon and Department of Mathematics, Higher Teachers

More information

RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS

RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS David García-García May 13, 2016 Faculdade de Ciências da Universidade de Lisboa OVERVIEW Random Matrix Theory Introduction Matrix ensembles A sample computation:

More information

Rational Solutions and Associated Special Polynomials for the Painlevé Equations

Rational Solutions and Associated Special Polynomials for the Painlevé Equations Rational Solutions and Associated Special Polynomials for the Painlevé Equations Peter A Clarkson Institute of Mathematics, Statistics and Actuarial Science University of Kent, Canterbury, CT2 7NF, UK

More information

Math 234. What you should know on day one. August 28, You should be able to use general principles like. x = cos t, y = sin t, 0 t π.

Math 234. What you should know on day one. August 28, You should be able to use general principles like. x = cos t, y = sin t, 0 t π. Math 234 What you should know on day one August 28, 2001 1 You should be able to use general principles like Length = ds, Area = da, Volume = dv For example the length of the semi circle x = cos t, y =

More information

New Exact Solutions for the (2+1)-dimensional Boiti-Leon-Manna-Pempinelli Equation

New Exact Solutions for the (2+1)-dimensional Boiti-Leon-Manna-Pempinelli Equation Applied Mathematical Sciences, Vol. 6, 2012, no. 12, 579-587 New Exact Solutions for the (2+1)-dimensional Boiti-Leon-Manna-Pempinelli Equation Ying Li and Desheng Li School of Mathematics and System Science

More information

Section 5.2 Series Solution Near Ordinary Point

Section 5.2 Series Solution Near Ordinary Point DE Section 5.2 Series Solution Near Ordinary Point Page 1 of 5 Section 5.2 Series Solution Near Ordinary Point We are interested in second order homogeneous linear differential equations with variable

More information

Roots and Coefficients Polynomials Preliminary Maths Extension 1

Roots and Coefficients Polynomials Preliminary Maths Extension 1 Preliminary Maths Extension Question If, and are the roots of x 5x x 0, find the following. (d) (e) Question If p, q and r are the roots of x x x 4 0, evaluate the following. pq r pq qr rp p q q r r p

More information

New explicit solitary wave solutions for (2 + 1)-dimensional Boussinesq equation and (3 + 1)-dimensional KP equation

New explicit solitary wave solutions for (2 + 1)-dimensional Boussinesq equation and (3 + 1)-dimensional KP equation Physics Letters A 07 (00) 107 11 www.elsevier.com/locate/pla New explicit solitary wave solutions for ( + 1)-dimensional Boussinesq equation and ( + 1)-dimensional KP equation Yong Chen, Zhenya Yan, Honging

More information

References 167. dx n x2 =2

References 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

Automatic Proofs of Identities: Beyond A=B

Automatic Proofs of Identities: Beyond A=B Bruno.Salvy@inria.fr FPSAC, Linz, July 20, 2009 Joint work with F. Chyzak and M. Kauers 1 / 28 I Introduction 2 / 28 Examples of Identities: Definite Sums, q-sums, Integrals n k=0 + 0 1 2πi n=0 ( ) H n

More information

Series Solutions of Linear ODEs

Series Solutions of Linear ODEs Chapter 2 Series Solutions of Linear ODEs This Chapter is concerned with solutions of linear Ordinary Differential Equations (ODE). We will start by reviewing some basic concepts and solution methods for

More information

Normalization integrals of orthogonal Heun functions

Normalization integrals of orthogonal Heun functions Normalization integrals of orthogonal Heun functions Peter A. Becker a) Center for Earth Observing and Space Research, Institute for Computational Sciences and Informatics, and Department of Physics and

More information

Asymptotics of Integrals of. Hermite Polynomials

Asymptotics 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 information

18.S34 (FALL 2007) PROBLEMS ON ROOTS OF POLYNOMIALS

18.S34 (FALL 2007) PROBLEMS ON ROOTS OF POLYNOMIALS 18.S34 (FALL 2007) PROBLEMS ON ROOTS OF POLYNOMIALS Note. The terms root and zero of a polynomial are synonyms. Those problems which appeared on the Putnam Exam are stated as they appeared verbatim (except

More information

ON THE ASYMPTOTIC REPRESENTATION OF THE SOLUTIONS TO THE FOURTH GENERAL PAINLEVÉ EQUATION

ON THE ASYMPTOTIC REPRESENTATION OF THE SOLUTIONS TO THE FOURTH GENERAL PAINLEVÉ EQUATION IJMMS 003:13, 45 51 PII. S016117103060 http://ijmms.hindawi.com Hindawi Publishing Corp. ON THE ASYMPTOTIC REPRESENTATION OF THE SOLUTIONS TO THE FOURTH GENERAL PAINLEVÉ EQUATION YOUMIN LU Received 5 June

More information

On rational approximation of algebraic functions. Julius Borcea. Rikard Bøgvad & Boris Shapiro

On 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 information

Evaluación numérica de ceros de funciones especiales

Evaluación numérica de ceros de funciones especiales Evaluación numérica de ceros de funciones especiales Javier Segura Universidad de Cantabria RSME 2013, Santiago de Compostela J. Segura (Universidad de Cantabria) Ceros de funciones especiales RSME 2013

More information

ELLIPSES AND ELLIPTIC CURVES. M. Ram Murty Queen s University

ELLIPSES AND ELLIPTIC CURVES. M. Ram Murty Queen s University ELLIPSES AND ELLIPTIC CURVES M. Ram Murty Queen s University Planetary orbits are elliptical What is an ellipse? x 2 a 2 + y2 b 2 = 1 An ellipse has two foci From: gomath.com/geometry/ellipse.php Metric

More information

arxiv:math/ v1 [math.ca] 21 Mar 2006

arxiv:math/ v1 [math.ca] 21 Mar 2006 arxiv:math/0603516v1 [math.ca] 1 Mar 006 THE FOURTH-ORDER TYPE LINEAR ORDINARY DIFFERENTIAL EQUATIONS W.N. EVERITT, D. J. SMITH, AND M. VAN HOEIJ Abstract. This note reports on the recent advancements

More information

ENCYCLOPEDIA OF MATHEMATICS AND ITS APPLICATIONS. Special Functions GEORGE E. ANDREWS RICHARD ASKEY RANJAN ROY CAMBRIDGE UNIVERSITY PRESS

ENCYCLOPEDIA OF MATHEMATICS AND ITS APPLICATIONS. Special Functions GEORGE E. ANDREWS RICHARD ASKEY RANJAN ROY CAMBRIDGE UNIVERSITY PRESS ENCYCLOPEDIA OF MATHEMATICS AND ITS APPLICATIONS Special Functions GEORGE E. ANDREWS RICHARD ASKEY RANJAN ROY CAMBRIDGE UNIVERSITY PRESS Preface page xiii 1 The Gamma and Beta Functions 1 1.1 The Gamma

More information

Lecture 5: Moment generating functions

Lecture 5: Moment generating functions Lecture 5: Moment generating functions Definition 2.3.6. The moment generating function (mgf) of a random variable X is { x e tx f M X (t) = E(e tx X (x) if X has a pmf ) = etx f X (x)dx if X has a pdf

More information

Qualification Exam: Mathematical Methods

Qualification Exam: Mathematical Methods Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin

More information

Haydock s recursive solution of self-adjoint problems. Discrete spectrum

Haydock s recursive solution of self-adjoint problems. Discrete spectrum Haydock s recursive solution of self-adjoint problems. Discrete spectrum Alexander Moroz Wave-scattering.com wavescattering@yahoo.com January 3, 2015 Alexander Moroz (WS) Recursive solution January 3,

More information

ENGI 9420 Lecture Notes 1 - ODEs Page 1.01

ENGI 9420 Lecture Notes 1 - ODEs Page 1.01 ENGI 940 Lecture Notes - ODEs Page.0. Ordinary Differential Equations An equation involving a function of one independent variable and the derivative(s) of that function is an ordinary differential equation

More information

Spectral Theory of Orthogonal Polynomials

Spectral Theory of Orthogonal Polynomials Spectral Theory of Orthogonal Barry Simon IBM Professor of Mathematics and Theoretical Physics California Institute of Technology Pasadena, CA, U.S.A. Lecture 2: Szegö Theorem for OPUC for S Spectral Theory

More information

A. MT-03, P Solve: (i) = 0. (ii) A. MT-03, P. 17, Solve : (i) + 4 = 0. (ii) A. MT-03, P. 16, Solve : (i)

A. MT-03, P Solve: (i) = 0. (ii) A. MT-03, P. 17, Solve : (i) + 4 = 0. (ii) A. MT-03, P. 16, Solve : (i) Program : M.A./M.Sc. (Mathematics) M.A./M.Sc. (Previous) Paper Code:MT-03 Differential Equations, Calculus of Variations & Special Functions Section C (Long Answers Questions) 1. Solve 2x cos y 2x sin

More information

1 Math 241A-B Homework Problem List for F2015 and W2016

1 Math 241A-B Homework Problem List for F2015 and W2016 1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let

More information

March Algebra 2 Question 1. March Algebra 2 Question 1

March Algebra 2 Question 1. March Algebra 2 Question 1 March Algebra 2 Question 1 If the statement is always true for the domain, assign that part a 3. If it is sometimes true, assign it a 2. If it is never true, assign it a 1. Your answer for this question

More information

Geometric RSK, Whittaker functions and random polymers

Geometric RSK, Whittaker functions and random polymers Geometric RSK, Whittaker functions and random polymers Neil O Connell University of Warwick Advances in Probability: Integrability, Universality and Beyond Oxford, September 29, 2014 Collaborators: I.

More information