Padé Approximations of and preservation of quadratic Lyapunov functions

We investigate whether or not quadratic Lyapunov functions are preserved under Padé approximations. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Automatic Deformation of Riemann–Hilbert Problems with Applications to the Painlevé II Transcendents

The stability and convergence rate of Olver’s collocation method for the numerical solution of Riemann–Hilbert problems (RHPs) are known to depend very sensitively on the particular choice of contours used as data of the RHP. By manually performing contour deformations that proved to be successful in the asymptotic analysis of RHPs, such as the method of nonlinear steepest descent, the numerical method can basically be preconditioned, making it asymptotically stable. In this paper, however, we will show that most of these preconditioning deformations, including lensing, can be addressed in an automatic, completely algorithmic fashion that would turn the numerical method into a black-box solver. To this end, the preconditioning of RHPs is recast as a discrete, graph-based optimization problem: the deformed contours are obtained as a system of shortest paths within a planar graph weighted by the relative strength of the jump matrices. The algorithm is illustrated for the RHP representing the Painlevé II transcendents.     

Location of Poles for the Hastings–McLeod Solution to the Second Painlevé Equation

We show that the well-known Hastings–McLeod solution to the second Painlevé equation is pole-free in the region \(\arg x \in [-\frac{\pi }{3},\frac{\pi }{3}]\cup [\frac{2\pi }{3},\frac{ 4 \pi }{3}]\), which proves an important special case of a general conjecture concerning pole distributions of Painlevé transcedents proposed by Novokshenov. Our strategy is to construct explicit quasi-solutions approximating the Hastings–McLeod solution in different regions of the complex plane and estimate the errors rigorously. The main idea is very similar to the one used to prove Dubrovin’s conjecture for the first Painlevé equation, but there are various technical improvements.     

Asymptotic Behavior of the Fourth Painlevé Transcendents in the Space of Initial Values

We study the asymptotic behavior of solutions of the fourth Painlevé equation as the independent variable goes to infinity in its space of (complex) initial values, which is a generalization of phase space described by Okamoto. We show that the limit set of each solution is compact and connected and, moreover, that any solution that is not rational has an infinite number of poles and infinite number of zeros.     

Poles of tritronquée solution to the Painlevé I equation and cubic anharmonic oscillator

The distribution of poles of zero-parameter solution to Painlevé I, specified by P. Boutroux as intégrale tritronquée, is studied in the complex plane. This solution has regular asymptotics ? $ \sqrt {{z \mathord{\left/ {\vphantom {z 6}} \right. \kern-0em} 6}} $ + O(1) as z → ∞, | arg z| < 4π/5. At the sector | arg z| > 4π/5 it is a meromorphic function with regular asymptotic distribution of poles at infinity. This fact together with numeric simulations for |z| < const allowed B. Dubrovin to make a conjecture that all poles of the intégrale tritronquée belong to this sector. As a step to prove this conjecture, we study the Riemann-Hilbert (RH) problem related to the specified solution of the Painlevé I equation. It is “undressed” to a similar RH problem for the Schrödinger equation with cubic potential. The latter determines all coordinates of poles for the intégrale tritronquée via a Bohr-Sommerfeld quantization conditions.     

An Algorithm for the Computation of Hermite–Padé Approximations to the Exponential Function: Divided Differences and Hermite–Padé Forms

We present the first of two different algorithms for the explicit computation of Hermite–Padé forms (HPF) associated with the exponential function. Some roots of the algebraic equation associated with a given HPF are good approximants to the exponential in some subsets of the complex plane: they are called Hermite–Padé approximants (HPA) to this function. Our algorithm is recursive and based upon the expression of HPF as divided differences of the function texp(xt) at multiple integer nodes. Using this algorithm, we find again the results obtained by Borwein and Driver for quadratic HPF. As an example, we give an interesting family of quadratic HPA to the exponential.     

Convergence of simultaneous Padé approximants to two dilogarithms

. .     

The convergence of diagonal Padé approximants and the Padé Conjecture

Questions related to the convergence problem of diagonal Padé approximants are discussed. A central place is taken by the Padé Conjecture (also known as the Baker-Gammel-Wills Conjecture). Partial results concerning this conjecture are reviewed and weaker and more special versions of the conjecture are formulated and their plausibility is investigated. Great emphasis is given to the role of spurious poles of the approximants. A conjecture by Nuttall (1970) about the number and distribution of such poles is stated and its importance for the Padé Conjecture is analyzed.     

Padé approximations to certain power series

Summary An explicit identity involvingQ _n (q ⁱ z) (i = 0, 1,, 4) is shown, whereQ _n (z) is the denominator of thenth Padé approximant to the functionf(z) = _k=0 q ^1/2k(k–1 Z ^k . By using the Padé approximations, irrationality measures for certain values off(z) are also given.

     

The A-Stability of Methods with Padé and Generalized Padé Stability Functions

This paper surveys some stability results and suggests the use of order arrows as an alternative to order stars in studying questions about the possible A-stability of a numerical method. A discussion of the so-called Butcher–Chipman conjecture includes a proof of a partial result.     

Padé approximants to certain elliptic-type functions

Among all continua joining non-collinear points a ₁, a ₂, a ₃ ∈ ?, there exists a unique compact Δ ? ? that has minimal logarithmic capacity. For a complex-valued non-vanishing Dini-continuous function h on Δ, we define $${f_h}(z): = \frac{1}{{\pi i}}\int_\Delta {\frac{{h(t)}}{{t - z}}\frac{{dt}}{{{w^ + }(t)}}} $$ , where $w(z): = \sqrt {\prod\nolimits_{k = 0}^3 {(z - {a_k})} } $ and w ⁺ is the one-sided value according to some orientation of 1. In this work, we present strong asymptotics of diagonal Padé approximants to f _h and describe the behavior of the spurious pole and the regions of locally uniform convergence from a generic perspective.     

Padé approximants of special functions

Let $ {f_{\gamma }}(x) = \sum\nolimits_{{k = 0}}^{\infty } {{{{T_k (x)}} \left/ {{{{\left( \gamma \right)}_k}}} \right.}} $ , where (??) _k =??(??+1) ? (??+k?1) and T _k (x)=cos (k arccos x) are Padé?CChebyshev polynomials. For such functions and their Padé?CChebyshev approximations $ \pi_{n,m}^{ch}\left( {x;{f_{\gamma }}} \right) $ , we find the asymptotics of decreasing the difference $ {f_{\gamma }}(x) - \pi_{n,m}^{ch}\left( {x;{f_{\gamma }}} \right) $ in the case where 0 ? m ? m(n), m(n) = o (n), as n???? for all x ?? [?1, 1]. Particularly, we determine that, under the same assumption, the Padé?CChebyshev approximations converge to f _?? uniformly on the segment [?1, 1] with the asymptotically best rate.     

Basic asymptotic expansions of solutions to the sixth Painlevé equation

Asymptotic expansions of solutions of three types are found for the sixth Painlevé equation in its three singular points. Translated from Sovremennaya Matematika i Ee Prilozheniya (Contemporary Mathematics and Its Applications), Vol. 57, Suzdal Conference–2006, Part 3, 2008.