Open Access
The swallowtail integral in the highly oscillatory region III
(Taylor & Francis, 2021) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística, Informática y Matemáticas; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
We consider the swallowtail integral Ψ(x,y,z):=∫∞−∞ei(t5+xt3+yt2+zt)dt for large values of |z| and bounded values of |x| and |y|. The integrand of the swallowtail integral oscillates wildly in this region and the asymptotic analysis is subtle. The standard saddle point method is complicated and then we use the modified saddle point method introduced in López et al., A systematization of the saddle point method application to the Airy and Hankel functions. J Math Anal Appl. 2009;354:347–359. The analysis is more straightforward with this method and it is possible to derive complete asymptotic expansions of Ψ(x,y,z) for large |z| and fixed x and y. The asymptotic analysis requires the study of three different regions for argz separated by three Stokes lines in the sector −π
Open Access
The swallowtail integral in the highly oscillatory region II
(Kent State University, 2020) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística, Informática y Matemáticas; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
We analyze the asymptotic behavior of the swallowtail integral R ∞ −∞ e i(t 5+xt3+yt2+zt)dt for large values of |y| and bounded values of |x| and |z|. We use the simpli ed saddle point method introduced in [López et al., 2009]. With this method, the analysis is more straightforward than with the standard saddle point method and it is possible to derive complete asymptotic expansions of the integral for large |y| and xed x and z. There are four Stokes lines in the sector (−π, π] that divide the complex y−plane in four sectors in which the swallowtail integral behaves di erently when |y| is large. The asymptotic approximation is the sum of two asymptotic series whose terms are elementary functions of x, y and z. One of them is of Poincaré type and is given in terms of inverse powers of y 1/2 . The other one is given in terms of an asymptotic sequence of the order O(y −n/9 ) when |y| → ∞, and it is multiplied by an exponential factor that behaves di erently in the four mentioned sectors. Some numerical experiments illustrate the accuracy of the approximation.
Open Access
Uniform convergent expansions of the error function in terms of elementary functions
(Springer, 2023) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
We derive a new analytic representation of the error function erfz in the form of a convergent series whose terms are exponential and rational functions. The expansion holds uniformly in z in the double sector | arg (±z) | <π/4. The expansion is accompanied by realistic error bounds.
Open Access
On a modifcation of Olver's method: a special case
(Springer US, 2016) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Ingeniería Matemática e Informática; Matematika eta Informatika Ingeniaritza
We consider the asymptotic method designed by Olver (Asymptotics and special functions. Academic Press, New York, 1974) for linear differential equations of the second order containing a large (asymptotic) parameter : xm y −2 y = g(x)y, with m ∈ Z and g continuous. Olver studies in detail the cases m = 2, especially the cases m = 0, ±1, giving the Poincaré-type asymptotic expansions of two independent solutions of the equation. The case m = 2 is different, as the behavior of the solutions for large is not of exponential type, but of power type. In this case, Olver’s theory does not give many details. We consider here the special case m = 2. We propose two different techniques to handle the problem: (1) a modification of Olver’s method that replaces the role of the exponential approximations by power approximations, and (2) the transformation of the differential problem into a fixed point problem from which we construct an asymptotic sequence of functions that converges to the unique solution of the problem. Moreover, we show that this second technique may also be applied to nonlinear differential equations with a large parameter.
Open Access
Orthogonal basis with a conicoid first mode for shape specification of optical surfaces
(Optical Society of America, 2016) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Navarro, Rafael; Ingeniería Matemática e Informática; Matematika eta Informatika Ingeniaritza
A rigorous and powerful theoretical framework is proposed to obtain systems of orthogonal functions (or shape modes) to represent optical surfaces. The method is general so it can be applied to different initial shapes and different polynomials. Here we present results for surfaces with circular apertures when the first basis function (mode) is a conicoid. The system for aspheres with rotational symmetry is obtained applying an appropriate change of variables to Legendre polynomials, whereas the system for general freeform case is obtained applying a similar procedure to spherical harmonics. Numerical comparisons with standard systems, such as Forbes and Zernike polynomials, are performed and discussed.
Open Access
Uniform convergent expansions of the Gauss hypergeometric function in terms of elementary functions
(Taylor & Francis, 2018) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Matematika eta Informatika Ingeniaritza; Institute for Advanced Materials and Mathematics - INAMAT2; Ingeniería Matemática e Informática
We consider the hypergeometric function 2F1(a, b; c; z) for z ∈ C \ [1,∞). For Ra ≥ 0, we derive a convergent expansion of 2F1(a, b; c; z) in terms of the function (1 − z)−a and of rational functions of z that is uniformly valid for z in any compact in C \ [1,∞). When a ∈ N, the expansion also contains a logarithmic term of the form log(1 − z). For Ra ≤ 0, we derive a convergent expansion of (1 − z)a 2F1(a, b; c; z) in terms of the function (1 − z)−a and of rational functions of z that is uniformly valid for z in any compact in C \ [1,∞) in the exterior of the circle |z − 1| = r for arbitrary r > 0. The expansions are accompanied by realistic error bounds. Some numerical experiments show the accuracy of the approximation.
Open Access
An asymptotic expansion of the hyberbolic umbilic catastrophe integral
(Springer, 2022) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
We obtain an asymptotic expansion of the hyperbolic umbilic catastrophe integral Ψ(H) (x,y,z) := ∫∞−∞∫∞−∞exp(i(s3+t3+zst +yt+xs))ds dt for large values of |x| and bounded values of |y| and |z|. The expansion is given in terms of Airy functions and inverse powers of x. There is only one Stokes ray at argx=π . We use the modified saddle point method introduced in (López et al. J Math Anal Appl 354(1):347–359, 2009). The accuracy and the asymptotic character of the approximations are illustrated with numerical experiments.
Open Access
A convergent version of Watson’s lemma for double integrals
(Taylor & Francis, 2022) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística, Informática y Matemáticas; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
A modification of Watson’s lemma for Laplace transforms ∞ 0 f(t) e−zt dt was introduced in [Nielsen, 1906], deriving a new asymptotic expansion for large |z| with the extra property of being convergent as well. Inspired in that idea, in this paper we derive asymptotic expansions of two-dimensional Laplace transforms F(x, y) := ∞ 0 ∞ 0 f(t,s) e−xt−ys dt ds for large |x| and |y| that are also convergent. The expansions of F(x, y) are accompanied by error bounds. Asymptotic and convergent expansions of some specialfunctions are given as illustration.
Open Access
Convergent and asymptotic methods for second-order difference equations with a large parameter
(Springer, 2018) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Matematika eta Informatika Ingeniaritza; Institute for Advanced Materials and Mathematics - INAMAT2; Ingeniería Matemática e Informática; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
We consider the second-order linear difference equation y(n+2)−2ay(n+1)−Λ2y(n)=g(n)y(n)+f(n)y(n+1) , where Λ is a large complex parameter, a≥0 and g and f are sequences of complex numbers. Two methods are proposed to find the asymptotic behavior for large |Λ|of the solutions of this equation: (i) an iterative method based on a fixed point method and (ii) a discrete version of Olver’s method for second-order linear differential equations. Both methods provide an asymptotic expansion of every solution of this equation. The expansion given by the first method is also convergent and may be applied to nonlinear problems. Bounds for the remainders are also given. We illustrate the accuracy of both methods for the modified Bessel functions and the associated Legendre functions of the first kind.
Embargo
Convergent and asymptotic expansions of the displacement elastodynamic integral in terms of known functions
(Elsevier, 2025-05-01) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2
The integral [Formula presented] plays an essential role in the study of several phenomena in the theory of elastodynamics (Ceballos and Prato, 2014). But an exact evaluation of this integral in terms of known functions is not possible. In this paper, we derive an analytic representation of this integral in the form of convergent series of elementary functions and hypergeometric functions. This series have an asymptotic character for either, small values of the variable s, or for small values of the variables r and R. It is derived by using the asymptotic technique designed in Lopez (2008) for Mellin convolution integrals. Some numerical experiments show the accuracy of the approximation supplied by the first few terms of the expansion.