International Journal of Computational Mathematics

Volume 2014, Article ID 587430, 17 pages

http://dx.doi.org/10.1155/2014/587430

## A Numerical Test of Padé Approximation for Some Functions with Singularity

^{1}Yamada Physics Research Laboratory, Aoyama 5-7-14-205, Niigata 950-2002, Japan^{2}Department of Physics, Ritsumeikan University, Noji-higashi 1-1-1, Kusatsu 525, Japan

Received 17 July 2014; Accepted 14 October 2014; Published 20 November 2014

Academic Editor: Don Hong

Copyright © 2014 Hiroaki S. Yamada and Kensuke S. Ikeda. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The aim of this study is to examine some numerical tests of Padé approximation for some typical functions with singularities such as simple pole, essential singularity, brunch cut, and natural boundary. As pointed out by Baker, it was shown that the simple pole and the essential singularity can be characterized by the poles of the Padé approximation. However, it was not fully clear how the Padé approximation works for the functions with the branch cut or the natural boundary. In the present paper, it is shown that the poles and zeros of the Padé approximated functions are alternately lined along the branch cut if the test function has branch cut, and poles are also distributed around the natural boundary for some lacunary power series and random power series which rigorously have a natural boundary on the unit circle. On the other hand, Froissart doublets due to numerical errors and/or external noise also appear around the unit circle in the Padé approximation. It is also shown that the residue calculus for the Padé approximated functions can be used to confirm the numerical accuracy of the Padé approximation and quasianalyticity of the random power series.

#### 1. Introduction

Padé approximation was introduced in mathematics by Hermite and Padé and it has been used in physics for more than 40 years ago [1–3]. In particular, there have been several important examples in physics, to which the Padé approximation was applied, such as summation of the divergent Rayleigh-Schrodinger perturbation series in scattering theory [4], critical phenomena in statistical physics [5–9], denoising from noisy data of time-series [10–14], and detection of singularity of phase space trajectories of Hamiltonian dynamical systems [15–19].

Mathematically, the Padé approximation can be used to estimate analyticity of functions. Indeed, the Padé approximation is usually superior to the truncated Taylor expansions when the original function contains any singularity. Let us consider a simple example. The function has brunch points at and . The domain of convergence is . Nevertheless, we can obtain an exact solution for when we apply diagonal Padé approximation to the function (See Section 3.3 for more details on this example) Although the mathematical validity of the Padé approximation has not been exactly proved yet, the Padé approximation is practically very useful to continue a singular function beyond the domain of convergence.

Let us consider a critical phenomenon for Ising model as a simple example in the statistical physics [5, 7, 8]. We assume that at the critical point the exact magnetic susceptibility has a singularity as where is a function of temperature and interactions and so on. In this case we sometimes use the logarithmic derivative of when we estimate the critical point and the critical exponent as a pole-type singularity as where denotes the derivative with respect to the variable . In the low temperature expansion for the magnetic susceptibility with some coefficients , we assume that the approximated susceptibility and the logarithmic derivative can be obtained as follows: Then we can estimate the critical point and the critical exponent by the Padé approximation to the truncated expansion . Note that the singularity on will be infinitely differentiable if the coefficients fall off sufficiently rapidly.

In general, analytic continuation over the singular point is possible along any other path in complex plane even if the function diverges at the singular point determining the radius of convergence, as seen in the above singularity of the Ising model. Therefore, the Padé approximation is useful to improve the convergence of the power series and approximate the exact solution.

Furthermore, the Padé approximation has been used to investigate convergence of Fourier series [2, 3] and the breakdown of KAM curves in complex plane for Hamiltonian map systems, which is described as the analytic domains of Lindstedt series for standard map [15–19].

In addition, we can see an interesting example concerning the Padé approximation in noisy data analysis [10–14]. The power series with finite random coefficient “almost always” has a natural boundary on the unit circle in the complex plane [20–24]. In a finite time-series, Froissart has shown that a natural boundary generated by the random time-series is approximated by doublets of poles and zeros (Froissart doublets) of the Padé approximated function, which are surrounding the vicinity of the unit circle. Taking advantage of the characteristic, the Padé approximation has been used in order to remove the noise and extract the true poles associated with damping modes from the observed noisy time-series.

The main purpose of the present paper is to investigate whether the Padé approximation is numerically useful for detecting the singularity of some test functions. In particular, we examine the usefulness for the functions with a natural boundary such as a lacunary power series and a random power series [2, 13, 18].

The organization of the paper is as follows. In Section 2, we give a brief explanation of the Padé approximation and some important reminders in the numerical calculation. In Section 3, we present some numerical results of the Padé approximation for test functions with branch cut, essential singularity. We also try to apply the Padé approximation to an entire function. In Section 4, application of the Padé approximation to some lacunary series which is known to have a natural boundary is given. In Section 5, numerical results of the application to random power series with a natural boundary and some test functions with random noise are also shown. In Section 6, we discuss the residue calculus for the Padé approximated functions to confirm the numerical errors and quasianalyticity of the random power series. In the last section, we give summary and discussion.

In Appendix A, the general result for Fibonacci generator used in Section 3.2 is given. In Appendix B, some theorems concerning zeros of polynomials are given. In Appendix C, some mathematical theorems for lacunary series, which are useful in reading the main text, are summarized. In Appendix D, exact Padé approximated functions to some lacunary power series with a natural boundary are given. Residue analysis for quasianalytic functions of Carleman class is given in Appendix E. Furthermore, some theorems concerning the random power series are given in Appendix F.

#### 2. Padé Approximation

In this section, we give preliminary instructions to the Padé approximation and some important reminders on the numerical calculation.

For a given function , a truncated Taylor expansion of order about zero is given as where denotes the coefficients of the Taylor expansion. Padé approximation is more accurate approximation for up to order than the Taylor expansion. The Padé approximation is a rational function, namely, a ratio of two polynomials, which agrees with the highest possible order with a truncated polynomial as follows: where is a polynomial of degree less than or equal to and is a polynomial of degree less than or equal to . Note that (normalized) here. A unique approximation can be specified for all choices of and such that when it exists. The coefficients can be obtained from the condition that the first terms vanish in the Taylor series in (4). The difference between the Padé approximation and the original function satisfies the following equation: In this paper, we sometimes use “ Padé approximation” or “ Padé approximated function” for . Solving the problem in (6) is called a linear Teoplitz problem which is generally ill-conditioned. We used full LU decomposition for the Teoplitz matrix in the problem as well as iterative improvement in order to eliminate the ill-posed problem [25, 26]. In addition, hereafter, we use the diagonal Padé approximation, that is, , of order to estimate the singularity of the test functions because of the convergence and limitation due to the round-off errors and other sources of numerical errors. Here, the singularity of the function is approximated by configuration of the poles and zeros of the order diagonal Padé approximated function. As mentioned in Introduction, in general, the Padé approximations are useful for representing unknown functions with possible poles. The application of the diagonal Padé approximation is insured for the functions with isolated singular points and rational-type functions. However, it is not fully clarified how the poles and zeros of the Padé approximated function describe brunch cut and natural boundary [27–29].

Generally, the magnitude of the residues associated with the spurious poles is much smaller than that with the true poles, and they are close to machine precision. Very recently, Gonnet et al. suggested an efficient algorithm for the Padé approximations [30, 31]. The algorithm detects and eliminates the spurious pole-zero pairs caused by the rounding errors by means of singular value decomposition for the Teoplitz matrix.

Before closing this section, we list up some important reminders when we numerically apply the Padé approximation to unknown functions as follows.(1)More accurate calculation becomes possible by scaling the expansion variable if there is a simple pole with large magnitude (≫1). That is, we should change the order of radius of convergence into by the scaling the expansion variable as in order to keep the numerical accuracy. This procedure is effective when we apply the Padé approximation to exponentially decaying coefficients with fluctuation.(2)Poles (i.e., roots of ) and zeros (i.e., roots of ) are sometimes cancelled (zero-pole ghost pairs). We can remove the effects of the ghost pairs and confirm the singularity of the functions by using the residue analysis of the Padé approximated function.(3)Poles and zeros of the Padé approximation to the truncated random power series accumulate around the unit circle as Froissart doublets. It is difficult to distinguish whether the poles of the Padé approximation originated from a natural boundary of the original function or the natural boundary generated by the numerical error and/or noise. Therefore, the numerical accuracy will be important to determine the coefficients of the Padé approximation.(4)In general, the denominators of the diagonal Padé approximated functions to the lacunary power series and the random power series become lacunary and random polynomials, respectively. Accordingly, the distribution of the poles and zeros of the approximated functions are similar to the distribution of the zeros corresponding to the original lacunary power polynomial and random power polynomial. In particular, it is well known that the zeros of the random polynomials uniformly distribute about the unit circle (see Appendix B).

#### 3. Examples of Padé Approximation for Some Functions

In this section we investigate the configuration of the poles and zeros of the Padé approximation to some test functions with singularity.

##### 3.1. Comparison of Padé Approximation with Taylor Expansion

First, we try the Padé approximation to the following test function with a brunch point at : The truncated Taylor expansion of the order around is The Padé approximation is given as

Figure 1 shows the approximated functions, the truncated Taylor series , and the original test function . The approximated function by the truncated Taylor expansion converges only within and deviates from the exact function for . On the other hand, it follows that the Padé approximated function well approximates the original function with very high precession even beyond the radius of convergence up to . Therefore, it is found that the divergent power series expansion (Taylor expansion) does still contain information about the original function outside the convergence radius, and rearranging the coefficients of the expansion into the Padé approximation recovers the information. As a result, the conversion from the Taylor form to the Padé form usually accelerates the convergence and often allows good accuracy even outside the radius of convergence of the power series.

##### 3.2. Padé Approximation for Fibonacci Generating Function

Let us consider the Fibonacci generating function , where Fibonacci sequence is encoded in the power series as the coefficients. The Fibonacci sequence is given by the following recursion relation: where and . Then the Fibonacci generating function becomes where and . The generating function has poles at and . Generating functions by more general recursion relation is given in Appendix A.

It should be noted that the diagonal Padé approximation for the truncated Fibonacci generation function has the following form: for even number . This means that the diagonal Padé approximation can detect the exact poles of the generating function irrespective of the order.

##### 3.3. Examples of Some Test Functions with Pole, Brunch Cut, and Essential Singularity

In this subsection, we use some test functions in applying of the Padé approximation: Here, has no singularity for , has a brunch cut along a line on , and has an essential singularity at . has eight poles at points on the unit circle .

First, let us apply Padé approximation to . In this case the explicit form of the Padé approximated function can be obtained in the following form [33]: Note that the coefficients of the numerator have always alternative sign, and the zeros and poles of the Padé approximated function are symmetrical to the imaginary axis with each other because . Figure 2(a) shows the numerical results in the complex -plane. All poles are on the left half-plane , and all zeros are on the right half-plane . The poles and zeros of the Padé approximated functions for the regular function go infinite and disappear as because the function is an entire function.

In Figure 2(b), distribution of the zeros and poles of the Padé approximated function to is shown. The poles and zeros make a line alternately between the two branch points of the ; and .

Figure 2(c) shows the distribution of the zeros and poles of the Padé approximation to . The Padé approximation clusters the poles and zeros at the singular point for . As the poles and zeros approach the singular point that reflects the essential singularity of the original function .

Figure 3(a) shows the poles and zeros of the Padé approximated function to the test function . The poles and zeros are alternatively distributed in eight directions from the origin. It seems that the distance between the poles and/or zeros on the same line becomes small as they approach the location of the true poles. Some “spurious poles” appear around the unit circle with the increase of the order of the Padé approximation, as seen in Figure 3(b), which are irrelevant poles due to insufficient numerical accuracy. It is found that the numerical accuracy of the Padé approximation fails for the higher order of the Padé approximation. We discuss the spurious poles in Sections 4 and 5 again.

#### 4. Natural Boundary of Lacunary Power Series

In this section, we examine the applicability of the Padé approximation to investigate the analyticity of some well-known test functions with a natural boundary on the unit circle . This will provide a preliminary information about what occurs in the Padé approximated functions for some test functions with a natural boundary. Indeed, we do know only a very few numerical examples which have a natural boundary and allow an exact diagonal Padé approximation.

The following are famous lacunary power series with a natural boundary on the unit circle , , , and , where the , , and are called after Jacobi, Weierstrass, and Kronecker. Some theorems for the lacunary series with a natural boundary are given in Appendix C [20–22].

Here, we use polar-form for the function by changing the variable, that is, , in order to simply display the functions as Then, note that the modulus works as a convergence factor of the series because it well converges for . Typically we take on the unit circle or inside the circle in the following numerical calculations.

##### 4.1. Example 1: Jacobi Lacunary Series

We try to apply Padé approximation to the function with a natural boundary on the unit circle . The Padé approximated function exactly has the following form: where the explicit form of the numerator is given in Appendix D. Accordingly, the poles of the Padé approximated function are given by roots of the polynomial This is also just a lacunary polynomial. In Figure 4 the numerical result of the Padé approximation for is shown. The poles and zeros are plotted for the Padé approximation in Figure 4(a). Inside the circle some cancellations of the ghost pairs appear. The poles and zeros accumulate around with the increase of order of the Padé approximation. In the case of the , the poles accumulate around with making the zero-pole pairings. Figure 4(b) shows the Padé approximated functions in the polar-form with . It well approximates the original function when the order of the Padé approximation increases.

It is also shown that the complex zeros of the polynomial (17) cluster near unit circle and distribute uniformly on the circle as by Erdos-Turan-type theorem given in Appendices C and B [34–41].

##### 4.2. Example 2: Fibonacci Lacunary Series

As a second example, we would like to apply Padé approximation to the following lacunary series: where is th Fibonacci number. This function also has a natural boundary on . The Padé approximated function exactly has the following form: The explicit form of the numerator is given in Appendix D. The poles of the Padé approximated function are given by zeros of the lacunary polynomial

In Figure 5 the numerical result of the Padé approximation to is shown. The poles and zeros are plotted for the Padé approximation in Figure 5(a). The poles and zeros accumulate around with the increase of the order of the Padé approximation. No pole appears inside the unit circle. The original function is also well approximated by the Padé approximation (see Figure 5(b)).

#### 5. Natural Boundary of Random Power Series and the Noise Effect on Padé Approximation

In this section, we apply Padé approximation to the random power series with a natural boundary with probability 1 and investigate how the approximation detect the singularity of the series. In addition, we examine the effect of noise on the coefficients of the power expansion for some test functions. Some related theorems for the natural boundary of the function generated by the random power series are given in Appendix F.

##### 5.1. Random Power Series and Natural Boundary

Let us consider a random power series: Here the coefficients, , are i.i.d. random variables which take a value within , and is the strength of the randomness. It is shown that, in general, the random power series has a natural boundary on the unit circle with probability one. Figure 6(a) shows distribution of poles and zeros of the Padé approximated function for . Some pairs of poles and zeros are perfectly cancelled inside the circle . On the other hand, almost all the poles and zeros of the Padé approximated function assemble around the circle and not cancelled. The pair of poles and zeros around the circle is called “Froissart doublets,” and it well corresponds to the natural boundary of . The original function is also well approximated by the Padé approximation (see Figure 6(b)).

Figure 7 shows an example of the coefficients of the random power series and the coefficients and of the Padé approximated function. The fluctuation of the coefficient that determines the poles of the Padé approximated function is smaller than that of the numerator.

Note that the truncated random series is a random polynomial. As for the random polynomial, it is well known that the distribution of the zeros converges on the unit circle when the order of the random polynomial increases (Erdos-Turan-type theorem) [34, 35, 40]. Accordingly, we can generally interpret that in the Padé approximated function to the random power series the distribution of poles and zeros also accumulates around the unit circle when the order of the Padé approximation increases. The dependence of the zeros of the random polynomial and the zeros and poles of the Padé approximation has been studied by Gilewicz and Kryakin [42] and Ding and Xiao [43].

##### 5.2. Effect of Noise on a Function with a Simple Pole

In the following subsections, we investigate influences of noise on the Padé approximation for some constructed noisy test functions as follows; where and . The is i.i.d. random variables within , where is the noise strength. Essentially, is the same as the random power series . First of all, in this subsection, we consider a truncated function with a simple pole. Note that if (constant) and , that is, in noise-free case, with a simple pole at . In [2] by Baker Jr., the noise effect is summarized as follows: the Padé approximation has an unstable zero at the distance of order from the origin, and the other zeros make Froissart doublets (zero-pole pairs).

Next, we consider a function with the noise strength . Note that with a simple pole at to clearly show the shift of the poles of the approximated function due to the noisy series.

Figure 8 shows distribution of the poles and zeros of the Padé approximated functions. It clearly shows the pole shift by the noise effect. In the noise-free case (), a pole of the Padé approximation appears at and the other poles are cancelled with zeros (zero-pole ghost pairs). In a case when the relatively small noise () is added, the poles and zeros move toward with making Froissart doublets, although a pole at is quite stable. It becomes impossible to detect the true pole at when the noise strength is relatively large (), not shown in Figure 8.

As a result, it is found that the locations of the ghost pairs are unstable for noise, and the residues for the poles are much smaller than one corresponding to the true pole. We can guess that the proximity of the nonmodal poles and zeros of the Padé approximated function can be understood in a sense that the poles due to the noise need zeros to cancel with each other as .

##### 5.3. Effect of Noise on a Function with a Branch Cut

We investigate the effect of the noise on functions with a branch cut. First, let us consider a function with an algebraic branch points at and and with the branch cut in . Distribution of the poles and zeros of the Padé approximated function is shown in Figure 9. In a case with relatively small noise (), some poles make a line on the branch cut, and some poles and zeros move toward the unit circle . It is impossible to detect the branch cut when the noise strength is relatively large ().

Next, let us consider a function with a logarithmic branch point at and with a brunch cut from to . The distribution of the poles and zeros of the Padé approximated function for the with the noisy perturbation is shown in Figure 10. Some poles and zeros are making a line alternatively on the branch cut in the noise-free case (). It assembles around the unit circle with making Froissart doublets when the noise with strength is added.

##### 5.4. Effect of Noise on a Function with a Natural Boundary

Figure 11 shows distribution of the poles and zeros of the Padé approximated function for which has a natural boundary on .

In the noise-free case, the pairs of poles and zeros of the Padé approximated function are perfectly cancelled inside the unit circle . The other poles and zeros of the Padé approximated function assemble around the circle without cancellation. In the relatively small noise case (), the location of the poles is not significantly changed compared with the zeros shifted outside the unit circle due to the noise effect. And, again, the poles and zeros move toward with making Froissart doublets when the noise strength is relatively large (). It is closely related to a fact that fluctuation of the coefficients of the numerator of the Padé approximated function is much larger than those in the denominator, as seen in Padé approximation to the random power series in Figure 7. As a result, the singularity of the Padé approximated function for the function with a natural boundary is more sensitive to the noisy perturbation than that in the functions with the other type singularity such as simple poles and branch points.

It is very difficult to effectively distinguish whether the poles of the Padé approximation originated from the natural boundary on of the original function or from the other natural boundary on generated by noisy series or numerical errors. Actually the round-off error affects the distribution of the poles and zeros of the Padé approximated function. Accordingly, to determine the expansion coefficients with adequate accuracy becomes very important in the numerical calculation. This is a drawback of the Padé approximation when we use it for functions with unknown singularities.

##### 5.5. Numerical Accuracy and Spurious Poles

As we observed in the last subsection, the effect of rounding error and accuracy limit of computers work in the numerical results of the Padé approximation. As the result of accumulation of the round-off error, the “spurious poles” appear around the unit circle as the pole-zero pairs when the order of Padé approximation increases (we used a term “Froissart doublets” for the poles-zero pairs generated by random power series, conveniently, although we cannot numerically distinguish it from the spurious poles due to the round-off errors; in the next section, we will discuss the Froissart doublets again).

However, we can roughly distinguish between true poles and the spurious poles by “residue analysis” of the Padé approximated function because the spurious poles-zero pairs are unstable for the change of the order. In this subsection, we try to investigate the residues of the Padé approximation for some test functions. Up to now, the residue analysis has been mainly used for performance comparison between the different algorithms of the Padé approximation of the same order [30, 31]. On the other hand, it seems that the study by using the information of the residue analysis is still rare in the Padé approximation [10, 12].

Generally, the rational polynomials of the diagonal Padé approximation can be uniquely identified by the poles and the corresponding residues as follows: where the residues are given by

Here, we investigate the convergence property of the magnitude of residues arranged in descending order.

Figure 12 shows the absolute value of the residues of some Padé approximated functions for the test function , which are arranged in descending order (note that they are noise-free cases.) The distribution of the poles and zeros of the Padé approximated functions is given in Figure 3. In a case of , the magnitude of the all residues is larger than , which correspond to the relevant poles arranged radially in eight directions from the true poles. On the other hand, in a case of , the spurious poles appear and distribute around the unit circle (see Figure 3(b)). It is found that the absolute values of the residues corresponding the spurious poles are several order of magnitude smaller than the relevant poles.

Distribution of poles and zeros of the Padé approximated function for the test function is shown in Figure 13. The stable poles and zeros are lined on , and the spurious poles appear around . The magnitude of the residues of the spurious poles is also enormously small compared with that of the stable poles remaining with the increase of the order of the Padé approximation.

Figure 14 is also the result of the residues analysis for the Padé approximated function for the test function with a natural boundary on the unit circle . In the Padé approximated function, the magnitude of the residues is shown in changing the noise strengths corresponding to poles-zeros distribution in Figure 11.

In the small noise case (), the results of the residue analysis for is almost the same as the noise-free case (), and in the case with relatively strong noise (), the noise shifts the magnitude of the residues with larger value. In addition, the result of the residue analysis of the noise-free cases for some different orders of the Padé approximation is shown in Figure 14(b). We should have in mind that the order is important when we apply Padé approximation to the lacunary power series because we should not take the order of the approximation in the gap of the series.

#### 6. Froissart Doublets

The problem of constructing the -transform of a finite time-series is a standard problem in mathematics [10–14]. For example, it is shown that for a sum of oscillating damped signals, the -transform associated with the time-series can be characterized by a sum of the poles of the Padé approximated function. The position of each pole is simply linked to the damping factor and the frequency of each of the oscillators. Also, it is important to note that all these poles lie strictly outside the unit circle because it corresponds to the damping [10–13]. In addition, we will consider quasianalyticity property of the random power series by the residue analysis of the Padé approximation.

##### 6.1. Noise Attractor

In signal processing, we can use the fact that the poles and zeros of the Padé approximated function to the noisy series distribute around the unit circle when we remove the noise from the observed data through -transform and/or Fourier transform of the data. Let a sequence be a sample signal without noise. Then we define the -transform of the sequence as The function is analytic interior of if the number of signals is finite [44]. Note that discrete Fourier transform is a special case of the -transform.

Next, let us consider a signal sequence in consisting of the superimposed damping oscillators as where is the amplitude of the th oscillator and . Here, and are the frequency and the damping factor of the th oscillator. Then, the -transform is where we take a limit keeping and . Accordingly, the singularity of appears as the poles at outside the unit circle , and the residue is .

On the other hand, let us consider a noise-added sequence . Then, the Froissart pointed out that there are two types of the poles: stable poles and unstable poles when we apply the diagonal Padé approximation to the unknown data set. In general, the -transform of the noisy sequence has a natural boundary on the unit circle with probability 1. In fact, the poles and zeros (Froissart doublets) of the Padé approximated function often distribute around the unit circle when the numerical error and/or noise are mixed into the Taylor series of the analytic functions, as seen in the last section. That is to say, we sometimes call the unit circle * noise attractor* in a sense that the poles and zeros are attracted to the circle as the Froissart doublets [45]. Accordingly, it is found that Padé approximated function for the function has stable poles associated with the damping modes and unstable spurious poles associated with the noisy fluctuation. After elimination of the spurious poles around the noise attractor from the noisy sequence we can reconstruct the noise-free sequence consisting of the stable poles located in the domain . Another remarkable feature of the nonmodal poles is that the absolute values of the Cauchy residues associated with them are usually much smaller than those associated with true poles.

##### 6.2. Random Power Series and Quasianalytic Function

Weierstrass defined the analytic function by direct analytic continuation of function. Then, apparently the analytic continuation is impossible beyond the natural boundary even if we can uniquely define the function and it is analytic outside the analytic domain. Borel and Gammel extended the narrow condition for the analyticity and gave a definition of* quasianalytic functions* [46, 47]. Gammel conjectured the following for the random power series [10, 45].

*Gammel Conjecture (1973).* The random power series belongs to the Borel class of quasianalytic functions as the following form:
where and are real numbers in the interval in and decreases rapidly with . Then, the natural boundary in the Weierstrass sense can be crossed.

The function (33) is a simple example that poles are densely distributed on the unit circle. Then the convergence property of the sequence is important for the analyticity of the function. Carleman proved that is quasianalytic if satisfies the following condition: This is Carleman class of quasianalytic functions. See Gammel’s paper [45] for the details. Moreover, Gammel and Nuttall proved that the quasianalytic functions can be exactly approximated by the Padé approximation [45].

*Gammel-Nuttall Theorem (1973).* If in (33) satisfies the condition (34) and , then the sequence of Padé approximation to the converges in measure to the function as in any closed, bounded region of the complex plane, where is a natural number that equals or less.

*Is the Gammel Conjecture True?* We try to examine the validity of the Gammel conjecture by applying residue analysis of the Padé approximated function to the random power series . Figure 15 shows the absolute values of the residues of the Padé approximated functions for three different samples in descending order. roughly exponentially decreases with respect to as
where is the decay exponent. It shows exponential decay (or faster), and, on the surface, supports the Gammal conjecture.

However, it is not nearly so simple. We should check the stability of the exponential-like decay of the magnitude of the residues by changing the order of the Padé approximation. Figure 16 shows the result for the three different orders: , , and . It expresses an indication that the decay exponent does not converge to a positive certain value. It seems that the exponent behaves as a limit . On the other hand, if we directly apply the Padé approximation to the quasianalytic function with , the exponent is stable for changing the order of the Padé approximation (see Appendix E). These facts suggest that the random power series does not belong to Carleman class of quasianalytic functions although it has a natural boundary on the unit circle, and it has the form (33). As a result, we can say that no optimism is warranted on the Gammel conjecture.

How does the residue analyses of the Padé approximation for the analyticity and/or quasianalyticity of unknown function work? It is an interesting and future problem.

#### 7. Summary and Discussion

In the present paper we numerically examined the effectiveness of the Padé approximation for some test functions with branch point, essential singularity, and natural boundary by watching the singularities of the Padé approximated functions. For the functions with a branch cut, the poles and zeros of the Padé approximated function are lined along the true branch cut. The poles and zeros are distributed around the true natural boundary if the original test function has a natural boundary. In addition, we gave the explicit Padé approximated functions for some lacunary power series which are useful to check the numerical result. It was shown that, in particular, the distribution of poles and zeros of the Padé approximated function for lacunary power series and the random power series accumulated around the unit circle when the order of the approximation increases.

We often suffer from the difficulty to distinguish whether or not the poles of the Padé approximation are intrinsically originated from the natural boundary of the original power series because the numerical errors contained in the expansion coefficients also yield a false natural boundary. Therefore, the expansion coefficients with adequate numerical accuracy are necessary when we apply the Padé approximation to functions with unknown singularities.

Furthermore, the residue calculus of the Padé approximated function is useful when we detect the singularity of the original power series from the asymptotic behavior of the truncated series. It is useful also for estimating the accuracy of the approximation. As a result, the residue calculus suggested that the random power series does not obey Gammel’s conjecture; that is, it does not belong to Borel class of the quasianalytic functions.

We finally remark that the most serious problem to be improved is the numerical accuracy due to the limitation of the order in the Padé approximation when we use it for detecting unknown singularities of wave functions in quantum physics [32].

#### Appendices

#### A. General Recursion Relation

We can construct a power series that has some pole-type singularities in the following form: where , , , , , and are real and for simplicity. Then the coefficients can be obtained by rearranging and comparing with the coefficients of the both sides in the same order as follows:

As a result, the power series with the pole-type singularities can be constructed by the recursion relation with , , and .

It becomes Fibonacci sequence when we set , , and .

#### B. Random Polynomial

The following theorems concerning the random power series are well known.

*Erdos-Turan-Type Theorem (1950).* Let us define a polynomial
where coefficients are randomly distributed and for simplicity. Then the zeros of the random polynomial cluster uniformly around the unit circle if “size of the truncated series” is small compared to the order of the polynomial, where

Note that this theorem also holds for the polynomials with deterministic coefficients such as Newman-type polynomial having coefficients in the sets or .

*Peres-Virag Theorem (2005).* Let be i.i.d. Gaussian-type random variables; then the distribution of the complex zeros of the power series
is

#### C. Some Gap Theorems of Lacunary Power Series

Weierstrass considered the analyticity of the power series where is a positive number. In the main text, we set , for . Then, it is proved that the function (C.1) has a natural boundary on the unit circle if the convergence radius of the function is unity based on the following theorems for the lacunary power series.

*Hadamard-Barck Gap Theorem (1892).* Let
where is a positive number and denote a strictly increasing sequence of the natural numbers, satisfying an inequality for . Then the function has a natural boundary on the unit circle .

*Fabry’s Gap Theorem (1899).* Power series
with radius of convergence has a natural boundary on the unit circle , provided that it is Fabry series; that is,

#### D. Numerators of Diagonal Padé Approximations for and

The diagonal Padé approximation for the truncated lacunary power functions and can be exactly executed as given in the main text. The numerators and of the Padé approximated functions can be given as follows: where .

Numerator of the diagonal Padé approximated function for is where , . means th Fibonacci number, and we set .

We have inductively obtained above results by means of Mathematica.

#### E. Residue Analysis for Carleman Class of Quasianalytic Functions

In this appendix, we give a direct result of residue analysis for “Carleman class” of the quasianalytic functions, for comparison with the other residue analyses in the main text. We apply the Padé approximation to the quasiperiodic function of the Carleman class, which is artificially constructed by a set of the poles as follows: where we set the poles at () on the unit circle. are i.i.d random variables in the interval and we take . Figure 17 shows the absolute values of the residues of the Padé approximated functions of order , , and for . They are arranged in descending order.

As a result, it seems that exponentially decreases with a stable exponent, regardless of the order of the Padé approximation. This supports that certainly, the Padé approximation is applicable to the quasianalytic functions in the Gammel conjecture as given in Gammel-Nuttall theorem. The Padé approximation for the quasianalytic function converges to the function even outside the unit circle. It should be also noted that in all cases the tails of are rapidly decay because the “truncated” series are essentially analytic functions.

#### F. Some Results for Natural Boundary in Noisy Series

In this appendix some theorems for the random power series are given. See, for example, [21] for the proofs.

*Steinhaus’s Theorem (1929).* Suppose that the power series
has radius of convergence . Let be a sequence of i.i.d. random variables in the interval . Then, with probability one, the random power series
has a natural boundary on , where .

*Paley-Zygmund Theorem (1932).* Suppose that the power series (F.1) has the radius of convergence 1. Let be a sequence of binary stochastic variables taking or with equal probability. Then, with probability one, the random power series
has a natural boundary on the unit circle .

The similar theorems can hold for random power series with a sequence of stochastic variables obeying i.i.d. in the interval or [48].

*Kahane’s Theorem (1985).* The circle of convergence is the natural boundary for random Taylor series (F.1) if the coefficients are independent and symmetric random variables.

The more generalized version has been given in the following form [22].

*Breuer-Simon Theorem (2011).* Suppose that the power series (F.1) has the convergence radius 1. Then for a.e. , has a strong natural boundary on if the is a stationary, ergodic, bounded, and nondeterministic process.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This paper was partially written for “International Symposium of Complexified Dynamics, Tunnelling and Chaos” held on 2005 in Kusatsu. This work is partly supported by Japanese people’s tax via MEXT, and the authors would like to acknowledge them. They are also very grateful to Dr. T. Tsuji and to Koike Memorial House for using the facilities during this study.

#### References

- G. A. Baker and J. L. Gammel,
*The Pade Approxima tion in Theoretical Physics*, Academic Press, New York, NY, USA, 1970. View at MathSciNet - G. A. Baker Jr.,
*Essentials of Pade Approximants*, Academic Press, New York, NY, USA, 1975. - J. Baker and P. Graves-Morris,
*Padé Approximants*, Cambridge University Press, Cambridge, UK, 2nd edition, 1996. View at Publisher · View at Google Scholar · View at MathSciNet - F. Sasagawa,
*Scattering Theory*, Syouk&bou, 1991, (Japanese). - H. Stanly,
*Introduction to Phase Transitions and Critical Phenomena*, Clarendon Press, Oxford, UK, 1971. - C. N. Yang and T. D. Lee, “Statistical theory of equations of state and phase transitions. I. Theory of condensation,”
*Physical Review*, vol. 87, p. 404, 1952. View at Publisher · View at Google Scholar · View at MathSciNet - B. Nickel, “On the singularity structure of the 2D Ising model susceptibility,”
*Journal of Physics A*, vol. 32, no. 21, pp. 3889–3906, 1999. View at Publisher · View at Google Scholar - R. Kubo, M. Toda, and N. Hashitsume,
*Statistical Physics II*, Springer, Berlin, Germany, 1985. View at Publisher · View at Google Scholar · View at MathSciNet - B. M. McCoy, “Do hard spheres have natural boundaries?” http://arxiv.org/abs/cond-mat/0103556.
- D. Bessis, “Padé approximations in noise filtering,”
*Journal of Computational and Applied Mathematics*, vol. 66, no. 1-2, pp. 85–88, 1996. View at Publisher · View at Google Scholar - H. Stahl, “The convergence of diagonal Padé approximants and the Padé conjecture,”
*Journal of Computational and Applied Mathematics*, vol. 86, no. 1, pp. 287–296, 1997. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - J. Gilewicz and M. Pindor, “Padé approximants and noise: a case of geometric series,”
*Journal of Computational and Applied Mathematics*, vol. 87, no. 2, pp. 199–214, 1997. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - D. Bessis and L. Perotti, “Universal analytic properties of noise: introducing the $J$-matrix formalism,”
*Journal of Physics A*, vol. 42, no. 36, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - L. A. Barbosa Coelho and L. A. Baccala, “Pade approximations as a modal identification technique,” in
*Proceedings of the 27th IMAC*, Orlando, Fla, USA, February 2009. - A. Berretti and L. Chierchia, “On the complex analytic structure of the golden invariant curve for the standard map,”
*Nonlinearity*, vol. 3, no. 1, pp. 39–44, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - C. Falcolini and R. de la Llave, “Numerical calculation of domains of analyticity for perturbation theories in the presence of small divisors,”
*Journal of Statistical Physics*, vol. 67, no. 3-4, pp. 645–666, 1992. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - R. de la Llave and S. Tompaidis, “Computation of domains of analyticity for some perturbative expansions of mechanics,”
*Physica D: Nonlinear Phenomena*, vol. 71, no. 1-2, pp. 55–81, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - A. Berretti and S. Marmi, “Scaling, perturbative renormalization and analyticity for the standard map and some generalizations,”
*Chaos, Solitons and Fractals*, vol. 5, no. 2, pp. 257–269, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - A. Berretti, C. Falcolini, and G. Gentile, “Shape of analyticity domains of Lindstedt series: the standard map,”
*Physical Review E*, vol. 64, no. 1, Article ID 015202, 2001. View at Publisher · View at Google Scholar · View at MathSciNet - T. W. Korner,
*Exercises for Fourier Analysis*, Cambridge University Press, Cambridge, UK, 1993. View at MathSciNet - R. Remmert,
*Classical Topics in Complex Function Theory*, Springer, New York, NY, USA, 1st edition, 1998. - J. Breuer and B. Simon, “Natural boundaries and spectral theory,”
*Advances in Mathematics*, vol. 226, no. 6, pp. 4902–4920, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - O. Knill and J. Lesieutre, “Analytic continuation of Dirichlet series with almost periodic coefficients,”
*Complex Analysis and Operator Theory*, vol. 6, no. 1, pp. 237–255, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - O. Costin and M. Huang, “Behavior of lacunary series at the natural boundary,”
*Advances in Mathematics*, vol. 222, no. 4, pp. 1370–1404, 2009. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery,
*Numerical Recipes in C*, Cambridge University Press, 1988. - W. H. Press and S. A. Teukolsky, “Padé approximants,”
*Computers in Physics*, vol. 6, p. 82, 1982. View at Google Scholar - J. Nuttall, “The convergence of Padé approximants of meromorphic functions,”
*Journal of Mathematical Analysis and Applications*, vol. 31, no. 1, pp. 147–153, 1970. View at Publisher · View at Google Scholar · View at MathSciNet - C. Pommerenke, “Padé approximants and convergence in capacity,”
*Journal of Mathematical Analysis and Applications*, vol. 41, no. 3, pp. 775–780, 1973. View at Publisher · View at Google Scholar · View at MathSciNet - H. Stahl, “Spurious poles in Padé approximation,”
*Journal of Computational and Applied Mathematics*, vol. 99, no. 1-2, pp. 511–527, 1998. View at Publisher · View at Google Scholar - P. Gonnet, R. Pachon, and L. N. Trefethen, “Robust rational interpolation and least-squares,”
*Electronic Transactions on Numerical Analysis*, vol. 38, pp. 146–167, 2011. View at Google Scholar · View at MathSciNet · View at Scopus - P. Gonnet, S. Güttel, and L. N. Trefethen, “Robust Padé Approximation via SVD,”
*SIAM Review*, vol. 55, no. 1, pp. 101–117, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - H. S. Yamada and K. S. Ikeda, “Analyticity of quantum states in one-dimensional tight-binding model,”
*The European Physical Journal B*, 2014. View at Publisher · View at Google Scholar - E. B. Saff and R. S. Varga, “On the zeros and poles of Padé approximants toez,”
*Numerische Mathematik*, vol. 25, no. 1, pp. 1–14, 1975. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - M. Kac, “On the average number of real roots of a random algebraic equation,”
*Bulletin of the American Mathematical Society*, vol. 49, pp. 314–320, 1943. View at Publisher · View at Google Scholar · View at MathSciNet - P. Erdos and P. Turan, “On the distribution of roots of polynomials,”
*Annals of Mathematics*, vol. 51, pp. 105–119, 1950. View at Publisher · View at Google Scholar · View at MathSciNet - F. Amoroso and M. Mignotte, “On the distribution of the roots of polynomials,”
*Annales de l'Institut Fourier*, vol. 46, no. 5, pp. 1275–1291, 1996. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A. Odlyzko and B. Poonen, “Zeros of polynomials with 0, 1 coefficients,”
*L'Enseignement Mathématique*, vol. 39, pp. 317–348, 1993. View at Google Scholar - B. Simon,
*Orthogonal Polynomials on the Unit Circle Part l :Classical Theory*, American Mathematical Society, 2004. - B. Simon,
*Orthogonal Polynomials on the Unit Circle Part 2 :Spectral Theory*, American Mathematical Society, Providence, RI, USA, 2004. View at MathSciNet - Y. Peres and B. Virág, “Zeros of the i.i.d. Gaussian power series: a conformally invariant determinantal process,”
*Acta Mathematica*, vol. 194, no. 1, pp. 1–35, 2005. View at Publisher · View at Google Scholar · View at MathSciNet - B. Simon,
*Szego's Theorem and Its Descendants: Spectral Theory for L2 Perturbations of Orthogonal Polynomials*, Princeton University Press, 2010. - J. Gilewicz and Y. Kryakin, “Froissart doublets in Padé approximation in the case of polynomial noise,”
*Journal of Computational and Applied Mathematics*, vol. 153, no. 1-2, pp. 235–242, 2003. View at Publisher · View at Google Scholar - X. Ding and Y. Xiao, “Natural boundary of random Dirichlet series,”
*Ukrainian Mathematical Journal*, vol. 58, no. 7, pp. 1129–1138, 2006. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - “We can also define the $Z$-transform by negative power ${z}^{-n}$. Then the function $Z(z)$ is analytic in outer domain of $\left|z\right|=1$, the poles corresponding to damping oscilla tions appear in the inside the unit circle $\left|z\right|<1$”.
- J. L. Gammel and J. Nuttall, “Convergence of Padé approximants to quasianalytic functions beyond natural boundaries,”
*Journal of Mathematical Analysis and Applications*, vol. 43, no. 3, pp. 694–696, 1973. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A. Shenitzer and N. Luzin, “Function: part I,”
*The American Mathematical Monthly*, vol. 105, no. 1, pp. 59–67, 1998. View at Publisher · View at Google Scholar · View at MathSciNet - N. Luzin, “Function: Part II,”
*The American Mathematical Monthly*, vol. 105, no. 3, pp. 263–270, 1998. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - J.-P. Kahane,
*Some Random Series of Functions*, vol. 5 of*Cambridge Studies in Advanced Mathematics*, Cambridge University Press, Cambridge, UK, 2nd edition, 1985. View at MathSciNet