An Analytic Representation of the Second Symmetric Standard Elliptic Integral in Terms of Elementary Functions

We derive new convergent expansions of the symmetric standard elliptic integral RD(x,y,z)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_D(x,y,z)$$\end{document}, for x,y,z∈C\(-∞,0]\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x, y,z\in {\mathbb {C}}{\setminus }(-\infty ,0]$$\end{document}, in terms of elementary functions. The expansions hold uniformly for large and small values of one of the three variables x, y or z (with the other two fixed). We proceed by considering a more general parametric integral from which RD(x,y,z)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_D(x,y,z)$$\end{document} is a particular case. It turns out that this parametric integral is an integral representation of the Appell function F1(a;b,c;a+1;x,y)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_1(a;b,c;a+1;x,y)$$\end{document}. Therefore, as a byproduct, we deduce convergent expansions of F1(a;b,c;a+1;x,y)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_1(a;b,c;a+1;x,y)$$\end{document}. We also compute error bounds at any order of the approximation. Some numerical examples show the accuracy of the expansions and their uniform features.


Introduction
In a recent paper [5] we derived new representations of the first symmetric standard elliptic integral R F (x, y, z) in the form of convergent expansions whose terms are elementary functions. These expansions are uniformly valid in large (and unbounded) regions of the complex plane of a selected variable x, y or z. In this work we continue this line of research and investigate uniformly convergent expansions of the second symmetric standard elliptic integral R D (x, y, z), that is defined as follows [7,Sec. 19.16,eq. 19.16.5], where, for simplicity in the exposition, we assume at this moment that either, x, y ∈ C\(−∞, 0] and z > 0 or y, z ∈ C\(−∞, 0] and x > 0 1 . In Sect. 5 we extend the results derived in this paper for positive z or x to complex values of these variables. The square roots in the denominator of (1) are assumed to be positive for positive argument. It is also reasonable to assume that the three variables are different, because otherwise this integral is an elementary function; for example R D (x, x, x) = 1/ √ x 3 . The integral is normalized in the form R D (1, 1, 1) = 1, and it is a homogeneous function of degree −3/2 in its three variables [7,Sec. 19.20,eq. 19.20.18]. This means that, for z > 0, is indeed a function of only two variables x and y. For x > 0, is also a function of only two variables y and z. For convenience in the analysis, in the remaining of the paper we consider the functions G 1 (x, y) and G 2 (y, z) instead of R D (x, y, z). All the results that we are going to derive in this paper for G 1 (x, y) and G 2 (y, z) can be translated to R D (x, y, z) by means of the connection formulas RD(x, y, z RD(x, y, z) y ,z∈ C\(−∞, 0]; x > 0. (5) The standard elliptic integrals [7] are special functions that have several applications in a large number of mathematical and physical problems. With respect to mathematical applications, we highlight their connection to the famous Theta functions and Weierstrass' elliptic function [29,Sec. 12.3]. They also play an important role in certain problems of geometry and statistics [18,28]. In reference to physical applications, the first elliptic integral appears in the computation of the period of a simple pendulum in a constant Vol. 77 (2022) An Analytic Representation of the Second Symmetric Page 3 of 24 171 gravitational field [29, Sec. 12.1.1]; the classification of limit cycles of several hamiltonian systems is directly related to the zeros of these integrals [32]; several problems of electromagnetic waves are solved in terms of elliptic integrals [34]. For certain geometries, the electric capacity of a conductor is written in terms of the inverse of R F (x, y, z) [27]. For other mathematical and physical applications of the standard elliptic integrals the reader is referred to [20]. A large set of properties and formulas for the standard elliptic integrals may be found in [6,7] and [29,Chap. 12]. In particular, connections between the standard elliptic integrals and the symmetric standard elliptic integrals, that are very important, as Carlson showed that the symmetric standard elliptic integrals are more appropriate for numerical purposes than the standard elliptic integrals [8][9][10][11][12].
As we have mentioned above, in this paper we are interested in the approximation of the integral R D (x, y, z) by a convergent series of elementary functions that is uniformly valid for large and small values of a certain selected variable. In the literature we can find several attempts to represent the symmetric standard elliptic integrals in the form of a series of elementary functions. Regarding asymptotic approximations, the first results were derived by Carlson, Gustafson [13] and Wong [35,Chap. 6,Sec. 7]. On the one hand, Gustafson obtained the first term of the asymptotic expansion of R D when one of its variables tends to zero or infinity [19]. This approximation was later improved by Carlson and Gustafson in [14]. On the other hand, new convergent expansions of R D (x, y, z) have been obtained in [17,21] and [22].
The expansions mentioned in the above paragraph are valid for real positive values of the three variables x, y, z; and the expansions are accurate when one of the variables is large compared to the other two. They are not accurate when two variables are of the same order. Therefore, they cannot be used when we need an approximation simultaneously valid for large and small values of one of the variables (and fixed values of the other two). An expansion of R D (x, y, z) uniformly convergent in one of the variables when the other two are restricted to a certain bounded domain may be found in [33]. The result is valid for positive values of the variables; error bounds are not given. In this paper we extend and generalize the result derived in [33] for the integral R D (x, y, z). We derive our expansions as an application of the theory of uniformly convergent expansions of integral transforms developed in [25].
As an illustration of the type of approximations that we are going to obtain in this paper (see Corollaries 4.1 and 4.2 below), we show, for example, the following approximation that is valid for 0 ≤ x < 1, y > 0: with |θ 1 (x)| ≤ 0.0497 x 4 ≤ 0.0497 . And the following approximation valid for 0 ≤ y < 1, z > 0: with |θ 2 (y)| ≤ 0.0694y 3 ≤ 0.0694. The paper is organized as follows. In Sect. 2 we show some preliminary results that will be required in the later analysis. Section 3 is devoted to the main results of the paper: uniformly convergent expansions of a certain integral from which G 1 (x, y) and G 2 (y, z) are particular cases. It turns out that this integral is an integral representation of the Appell function F 1 (a; b, c; a+1; x, y) [1,Sec. 16.13]. Then, the results derived in Sect. 3 are applied in Sect. 4 to G 1 (x, y), G 2 (y, z) and F 1 (a; b, c; a + 1; x, y), deriving convergent expansions of these functions that are uniformly valid for one of their variables on large subsets of the complex plane. In Sect. 5 we eliminate the restriction z > 0 (or x > 0) for R D (x, y, z) and let z ∈ C\(−∞, 0] (or x ∈ C\(−∞, 0]) by reducing, on the other hand, the domain for the variables x and y (or z and y) to smaller sectors inside C\(−∞, 0]. We complete this paper by checking the accuracy of these expansions with some numerical examples. Throughout the paper, for any complex variable w, arg w ∈ (−π, π] denotes its main argument and square roots are assumed to take their principal value.

Preliminaries
With the exception of the expansion [33], the expansions mentioned above are derived from the integral definition (1) of the second symmetric standard elliptic integral R D (x, y, z). They follow by applying the standard techniques of the theory of asymptotic expansions of integrals to the integral (1), [23], [31,Chap 16], [35,Chaps. 3 and 6]. Then, whenever they are convergent or only asymptotic, those expansions are accurate when one of the variables is large compared to the other two. Therefore, they are not uniformly valid in large regions of the complex plane that include large and small values of any selected variable. This restriction may be avoided by using a new analytic technique for deriving uniform expansions of integral transforms, introduced in [25]. In fact, this idea has been previously used in [5] to derive uniform expansions of R F (x, y, z) in terms of elementary functions, and in [3,4,15,16,24] to deduce uniform expansions of several other special functions. In this section we summarize the uniform technique introduced in [25], and that we are going to apply to R D (x, y, z) and F 1 (a; b, c; a + 1; x, y) below. Consider the integral transform of a function g(t) with kernel h(t, y) of the form: where D is a certain unbounded region of the complex plane that contains the point y = 0, and with the following assumptions for the functions h and g: It has been shown in [25] that, when we replace g(t) in (8) by its Taylor expansion at an appropriate point w ∈ Ω, where r n (t) is the Taylor remainder, and interchange sum and integral in (8), we obtain an expansion of G(y), with the following three properties: the expansion is uniform for y ∈ D: for any order n of the approximation, the absolute error satisfies the bound |R n (y)| ≤ C n for any y ∈ D with C n independent of y, the expansion is convergent in D, and the terms of the expansion Φ k (y) are elementary functions of y.
The uniform technique described above requires the integration interval in (8) to be bounded. Consequently, in order to apply the above technique to G 1 (x, y) and G 2 (y, z), instead of (2) and (3), we need an integral representation of these functions defined on a bounded interval. With this aim, we introduce in (2) and (3) the change of variable s → t defined in the form 1 + s = 1/t to obtain For the sake of generality and convenience, we investigate in the next section uniformly convergent expansions of the more generalized integral that is indeed an integral representation of the Appell function F 1 (a; b, c; a . It only takes a little more effort, and the results that we derive for this integral may be applied, not only to the symmetric integral R D (x, y, z), but also to the first Appell function and other special functions, like for example the first symmetric elliptic integral R F (x, y, z). Therefore, from the uniformly convergent expansions of the integral (10) that we are going to derive in the next section, we will obtain as corollaries in Sect. 4, new uniformly convergent expansions of R D (x, y, z) and F 1 (a; b, c; a + 1; x, y). 2 The main results of the paper are given in Theorems 3.1 and 3.2 in the next section. But firstly we will give two preliminary lemmas that we need in the later analysis. These lemmas are similar to [5, Lemmas 1 and 2] (see also [16] for a Proof of Lemma 2.1). For that reason their proofs are omitted here.
Vol. 77 (2022) An Analytic Representation of the Second Symmetric Page 7 of 24 171 Figure 1. The region S(θ) given by (11) is marked in green Observation 2.1. Observe that, despite the overlapping in the regions defining w(x), the inequalities (13) hold whenever w(x) is given by any of the three lines in the right hand side of (12).

A Uniformly Convergent Expansion of F (a, b, c; x, y)
We now apply the theory of uniform expansions of integral transforms introduced in [25] and condensed in Sect. 2 to the integral (10). We select y as the uniform variable, corresponding to the exponent b in (10).
Theorem 3.1. For any fixed angle θ ∈ [π/2, π) consider the region S(θ) ⊂ C\(−∞, −1] given in (11). Then, for any a, b, c ∈ C with a, b ≥ 0, c > −1; x ∈ C\(−∞, −1], y ∈ S(θ), and n = 1, 2, 3, ..., the integral (10) admits the following representation with where (a) k denotes the Pochhamer's symbol [2] and 2 F 1 is the Gauss hypergeometric function [26]. Furthermore, w(x) is given in the first or second line of (12) and the remainder term is bounded in the form The right hand side of (14) is a uniform convergent expansion of F (a, b, c; x, y) with an exponential order of convergence: as n → ∞, Proof. We apply the uniform theory summarized in Sect. 2 for the integral (8) to the integral (10), with the identification and Also, from Lemma 2.1, for y ∈ S(θ) with H integrable on [0, 1]. For any w ∈ Ω and n = 1, 2, 3, ..., we have that We choose 3 w = w(x) given by the first or second line of (12).
3 There are other possible choices of w, but the value w = w(x) not only assures the convergence of expansion (14), it also minimizes the value of the factor x(w−t) After replacing g(t) in the integral (10) by the right hand side of (19) and exchanging summation and integration, we derive the right hand side of (14) with R n (a, b, c; x, y When we replace r n (t; x, a) by the right hand side of (20) into the above integral, we derive the bound (16) after using the standard series representation of the Gauss hypergeometric function, Lemma 2.2 and the bound |t − w| ≤ |w| for t ∈ [0, 1].
Using the asymptotic behavior of the quotient of two gamma functions [2, eq. 5.11.12] and the asymptotic behavior of the Gauss hypergeometric function [30, eq. (15)], we find (17), which proves the convergence of (14).
Finally, the uniform character of expansion (14) follows from the fact that the right hand side of (16) does not depend on y. In the following theorem, instead of considering the base point w = w(x) given in Lemma 2.2, we consider the base point w = 0. This alternative imposes a more demanding restriction on x (|x| < 1), but it gives the simplest possible expansion for the integral (10).

Theorem 3.2.
For any fixed angle θ ∈ [π/2, π) consider the region S(θ) ⊂ C\(−∞, −1] given in (11). Then, for a, b, c ∈ C with a, b ≥ 0, c > −1; |x| < 1, y ∈ S(θ) and n = 1, 2, 3, ..., the integral (10) can be written in the form a, b, c; x, y). The remainder term is bounded in the form 4 The right hand side of (21) is a uniform convergent expansion of F (a, b, c; x, y) with an exponential order of convergence: as n → ∞, We also have the following particular bounds: 1. For x ≥ 0 and 0 < a < 1: 2. For 0 < a < 1: Proof. The Proof of (21)  For 0 < a < 1 the hypergeometric 3 F 2 in (22) can be bounded by the same hypergeometric function but evaluated at |x| = 1. Then, bound (25) follows from the bound [4, pag. 1776] valid for α, β > 0 and n > β − α. Finally, using the asymptotic behavior of the quotient of two gamma functions and the asymptotic behavior of the 3 F 2 hypergeometric function, we find (23), which proves the convergence of (21). The uniformity feature follows from the fact that the right hand side of (22) does not depend on y.  F (a, b, c;

x, y)
In this section we consider some specific applications of the main results of the last section, Theorems 3.1 and 3.2. The first one is the derivation of two uniform expansions of the second symmetric standard elliptic integral R D (x, y, z). For this purpose we apply theorems 3.1 and 3.2 to the integrals G 1 (x, y) and G 2 (y, z) in the form given in (9), and related to R D (x, y, z) by (4) and (5). The second application is the derivation of a uniform expansion of the Appell function F 1 (a; b, c; a + 1; x, y). Finally, we comment the possibility of deriving uniform expansions of the functions R F (x, y, z) and R G (x, y, z).

The remainder term is bounded in the form
• For |x| < 1, The remainder term is bounded in the form The approximants in (26) and (27) are elementary functions: Either, the right hand side of (26) or of (27) is a uniform convergent expansion of G 1 (x, y) with an exponential order of convergence.
Proof. From (9) and (10) it is clear that Then, all the theses of this corollary are particular cases of Theorems 3.1 and 3.2. We note that the coefficients A n (1/2, 1/2; y, w) given in (28)
Proof. From (9) and (10) it is clear that Then, all the theses of this corollary are particular cases of Theorems 3.1 and 3.2. We note that the coefficients A n (3/2, 1/2; z, w) given in (34) are elementary functions as shown in (35), which can be proved by repeatedly integrating by parts on the integral representation of the Gauss hypergeometric function [26, eq. 15.6.1] on the left hand side of (35).

A Uniformly Convergent Expansion of the Appell Function
Recall the integral representation of Appell function (see [1,Sec 16.15, eq 16.15.1]) We have the following corollary.
Proof. Comparing (36) to (10) we have that Then, the results follow from Theorems 3.1 and 3.2. R F (x, y, z) in terms of elementary functions are given in [5], and we refer the reader to this reference for details. Alternatively, those expansions can be derived from Theorems 3.1 and 3.2 by using the relation

Note 4.1. Uniform expansions of the first symmetric standard elliptic integral
The function R G (x, y, z), [7,Sec 19.21,eq 19.21.11], may be written in the form

A Larger Domain of Applicability of Corollaries 4.1 and 4.2
An analytic representation of the first symmetric elliptic integral R F (x, y, z) is given in [5]. That expansion is first derived for z > 0, and later extended to complex values of z in [5,Sec. 4] by using analytic continuation arguments. Identical arguments may be used here to enlarge the range of applicability of Corollaries 4.1 and 4.2 and deduce that the expansions derived for G 1 (x, y) in Corollary 4.1 hold, not only for z > 0 and x, y ∈ C\(−∞, 0], but in the bigger domain:    (30) (32)   Similarly, the expansions derived for G 2 (y, z) in Corollary 4.2 hold, not only for x > 0 and y, z ∈ C\(−∞, 0], but in the bigger domain: Λ 2 := {(x, y, z) ∈ (C\(−∞, 0]) 3 ; | arg y − arg x| < π, | arg z 3 − arg x 3 | < π}.

Numerical Experiments
Finally, in Tables 1 and 2 and Fig. 3, we give some numerical experiments that illustrate the accuracy and uniform character of the expansions derived in Theorems 3.1 and 3.2 (and therefore in Corollaries 4.1-4.3). In the numerical Data Availibility Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

Declarations Competing Interests
The authors have no relevant financial or non-financial interests to disclose.
Open Access. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.