Multivariate Frobenius–Padé approximants: Properties and algorithms

The aim of this paper is to construct rational approximants for multivariate functions given by their expansion in an orthogonal polynomial system. This will be done by generalizing the concept of multivariate Padé approximation. After defining the multivariate Frobenius–Padé approximants, we will be interested in the two following problems: the first one is to develop recursive algorithms for the computation of the value of a sequence of approximants at a given point. The second one is to compute the coefficients of the numerator and denominator of the approximants by solving a linear system. For some particular cases we will obtain a displacement rank structure for the matrix of the system we have to solve. The case of a Tchebyshev expansion is considered in more detail.     

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.     

On thermal stability of a reactive third-grade fluid in a channel with convective cooling the walls

In this paper, the thermal stability of a reactive third-grade liquid flowing steadily between two parallel plates with symmetrical convective cooling at the walls is investigated. The system is assumed to exchange heat with the ambient following Newton’s cooling law and the reaction is exothermic under Arrhenius kinetics, neglecting the consumption of the material. Approximate solutions are constructed for the governing nonlinear boundary value problem using regular perturbation techniques together with a special type of Hermite-Padé approximants and important properties of the velocity and temperature fields including bifurcations and thermal criticality conditions are discussed. It is observed that a combined increase in non-Newtonian parameter and convective cooling enhances the thermal stability of the material.     

Anti-Gaussian Padé Approximants

This paper gives a synthesis of Padé approximants and anti-Gaussian quadratures. New rational approximants for Stieltjes series have been constructed. In addition, a three term recurrence relation is given for the numerator and denominator, which is useful when the given functional is not defin ite positive.We give the different algebraic properties of these new polynomials, which are similar to those obtained with the Gaussian quadrature formula. We find an easy definition and several relations with Padé approximants. Finally, some numerical results are given in the last section.     

Multivariate Padé approximants revisited

Several definitions of multivariate Padé approximants have been introduced during the last decade. We will here consider all types of definitions based on the choice that the coefficients in numerator and denominator of the multivariate Padé approximant are defined by means of a linear system of equations. In this case a determinant representation for the multivariate Padé approximant exists. We will show that a general recursive algorithm can be formulated to compute a multivariate Padé approximant given by any definition of this type. Here intermediate results in the recursive computation scheme will also be multivariate Padé approximants. Up to now such a recursive computation of multivariate Padé approximants only seemed possible in some special cases.     

Quasipower and hypergeometric series — construction and evaluation

It is proved that some power series converging very slowly in a neighbourhood of the point 1 can be transformed intoquasipower series. The latter converge faster but are more complicated because they contain some hypergeometric series₂ F ₁. Standard methods of the values evaluation for needed hypergeometric series with the aid of recurrence relations are not sufficiently efficient for some variable values. Therefore a new method, formally similar to Levin's transforms, is proposed. More generally, this is a method of approximative evaluating of such a solution of an inhomogeneous recurrence relation of order one which has some particular asymptotic properties.The efficacity of the proposed methods is analyzed in detail for Euler's dilogarithm. This is a typical function whose power series is approached with difficulties ifz1. In particular, its Padé approximants are sufficiently accurate only for, sayx[–1, 1/2]. Hermite-Padé approximation is more effective. Resulting irrational approximants generalize in some sense partial sums of the quasipower series introduced here.     

Defects and the Convergence of Padé Approximants

A selective survey is given of convergence results for sequences of Padé approximants. Various approaches for dealing with the convergence problems due to `defects" are discussed. Attention is drawn to the close relationship between analyticity properties of a function and the `smoothness" of its Taylor series coefficients. A new theorem on the convergence of horizontal sequences of Padé approximants to functions in the Baker–Gammel–Wills conjecture function class is presented.     

On the convergence of type I Hermite–Padé approximants

Padé approximation has two natural extensions to vector rational approximation through the so-called type I and type II Hermite–Padé approximants. The convergence properties of type II Hermite–Padé approximants have been studied. For such approximants Markov and Stieltjes type theorems are available. To the present, such results have not been obtained for type I approximants. In this paper, we provide Markov and Stieltjes type theorems on the convergence of type I Hermite–Padé approximants for Nikishin systems of functions.     

General order Newton-Padé approximants for multivariate functions

Summary Padé approximants are a frequently used tool for the solution of mathematical problems. One of the main drawbacks of their use for multivariate functions is the calculation of the derivatives off(x ₁, ...,x _p ). Therefore multivariate Newton-Padé approximants are introduced; their computation will only use the value off at some points. In Sect. 1 we shall repeat the univariate Newton-Padé approximation problem which is a rational Hermite interpolation problem. In Sect. 2 we sketch some problems that can arise when dealing with multivariate interpolation. In Sect. 3 we define multivariate divided differences and prove some lemmas that will be useful tools for the introduction of multivariate Newton-Padé approximants in Sect. 4. A numerical example is given in Sect. 5, together with the proof that forp=1 the classical Newton-Padé approximants for a univariate function are obtained.     

An algebraic model for the center problem

We construct an algebraic model for the Center Problem for equation . This problem is related to the classical Poincaré Center-Focus problem for polynomial vector fields.     

A multivariate QD-like algorithm

The problem of constructing a univariate rational interpolant or Padé approximant for given data can be solved in various equivalent ways: one can compute the explicit solution of the system of interpolation or approximation conditions, or one can start a recursive algorithm, or one can obtain the rational function as the convergent of an interpolating or corresponding continued fraction.In case of multivariate functions general order systems of interpolation conditions for a multivariate rational interpolant and general order systems of approximation conditions for a multivariate Padé approximant were respectively solved in [6] and [9]. Equivalent recursive computation schemes were given in [3] for the rational interpolation case and in [5] for the Padé approximation case. At that moment we stated that the next step was to write the general order rational interpolants and Padé approximants as the convergent of a multivariate continued fraction so that the univariate equivalence of the three main defining techniques was also established for the multivariate case: algebraic relations, recurrence relations, continued fractions. In this paper a multivariate qd-like algorithm is developed that serves this purpose.     

Rational approximation to the Thomas-Fermi equations

We show that a simple and straightforward rational approximation to the Thomas-Fermi equation provides the slope at origin with unprecedented accuracy and that Padé approximants of relatively low order are far more accurate than more elaborate approaches proposed recently by other authors. We consider both the Thomas-Fermi equation for isolated atoms and for atoms in strong magnetic fields.     

On a closed orbit of some constrained system

We study on what one calls a constrained system of ODEs on It consists of two ordinary differential equations and an algebraic equation with respect to three unknown functions. We seek closed orbits of such a system. A necessary and sufficient condition for the system to have non-trivial closed orbits is given. Elementary tools such as Lyapunov functions and Poincaré’s index theory are used in the proof of the result. As an application we consider a constrained system associated with a non-constraint system of ODEs called the modified Bonhöffer-van der Pol system.     

Classification of quadratic systems admitting the existence of an algebraic limit cycle

In the paper we find a set of necessary conditions that must be satisfied by a quadratic system in order to have an algebraic limit cycle. We find a countable set of ?5 parameter families of quadratic systems such that every quadratic system with an algebraic limit cycle must, after a change of variables, belong to one of those families. We provide a classification of all the quadratic systems which can have an algebraic limit cycle based on geometrical properties of the embedding of the system in the Poincaré compactification of R². We propose names for all the classes we distinguish and we classify all known examples of quadratic systems with algebraic limit cycle. We also prove the integrability of certain classes of quadratic systems.     

Analytic continuation of Taylor series and the boundary value problems of some nonlinear ordinary differential equations

We compare and discuss the respective efficiency of three methods (with two variants for each of them), based respectively on Taylor (Maclaurin) series, Padé approximants and conformal mappings, for solving quasi-analytically a two-point boundary value problem of a nonlinear ordinary differential equation (ODE). Six configurations of ODE and boundary conditions are successively considered according to the increasing difficulties that they present. After having indicated that the Taylor series method almost always requires the recourse to analytical continuation procedures to be efficient, we use the complementarity of the two remaining methods (Padé and conformal mapping) to illustrate their respective advantages and limitations. We emphasize the importance of the existence of solutions with movable singularities for the efficiency of the methods, particularly for the so-called Padé-Hankel method. (We show that this latter method is equivalent to pushing a movable pole to infinity.) For each configuration, we determine the singularity distribution (in the complex plane of the independent variable) of the solution sought and show how this distribution controls the efficiency of the two methods. In general the method based on Padé approximants is easy to use and robust but may be awkward in some circumstances whereas the conformal mapping method is a very fine method which should be used when high accuracy is required.     

A Rational Approximant for the Digamma Function

Power series representations for special functions are computationally satisfactory only in the vicinity of the expansion point. Thus, it is an obvious idea to use Padé approximants or other rational functions constructed from sequence transformations instead. However, neither Padé approximants nor sequence transformation utilize the information which is avaliable in the case of a special function – all power series coefficients as well as the truncation errors are explicitly known – in an optimal way. Thus, alternative rational approximants, which can profit from additional information of that kind, would be desirable. It is shown that in this way a rational approximant for the digamma function can be constructed which possesses a transformation error given by an explicitly known series expansion.     

Convergence Rates of Padé and Padé-Type Approximants

A comparison is made between Padé and Padé-type approximants. LetQ_nbe thenth orthonormal polynomial with respect to a positive measureμwith compact support inC. We show that for functions of the form[formula]wherewis an analytic function on the support ofμ, Padé-type approximants with denominatorQ_ngive a successful and, in general, better approximation procedure than Padé approximation.     

Generalized moment representations and Padé-Chebyshev approximations

An approach to the application of Dzyadyk's generalized moment representations in problems of construction and investigation of the Padé-Chebyshev approximants is developed. With its help, certain properties of the Padé-Chebyshev approximants of a class of functions that is a natural analog of the class of Markov functions are studied. In particular, it is proved that the poles of the Padé-Chebyshev approximants of these functions lie outside their domain of analyticity.Translated from Ukrainskii Matematicheskii Zhurnal, Vol. 42, No. 6, pp. 762–766, June, 1990.     

Two-point Padé approximants for the expansions of Stieltjes functions in real domain

The fundamental inequalities for the sequences of subdiagonal and diagonal one-point Padé approximants to Stieltjes function has been extended to the case of certain two-point Padé approximants. The results can be applied to the theory of inhomogeneous media for calculating the bounds for the effective transport coefficients of two-components heterogeneous materials.