Source author record

Amparo Gil

Amparo Gil appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
6topics
2close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2020arXiv

Asymptotic expansions of Jacobi polynomials and of the nodes and weights of Gauss-Jacobi quadrature for large degree and parameters in terms of elementary functions

Asymptotic approximations of Jacobi polynomials are given in terms of elementary functions for large degree $n$ and parameters $α$ and $β$. From these new results, asymptotic expansions of the zeros are derived and methods are given to obtain the coefficients in the expansions. These approximations can be used as initial values in iterative methods for computing the nodes of Gauss--Jacobi quadrature for large degree and parameters. The performance of the asymptotic approximations for computing the nodes and weights of these Gaussian quadratures is illustrated with numerical examples.

preprint2020arXiv

Numerical evaluation of Airy-type integrals arising in uniform asymptotic analysis

We describe a method to evaluate integrals that arise in the asymptotic analysis when two saddle points may be close together. These integrals, which appear in problems from optics, acoustics or quantum mechanics as well as in a wide class of special functions, can be transformed into Airy-type integrals and we use the trapezoidal rule to compute these integrals numerically. The quadrature method, which remains valid when two saddle points coalesce, is illustrated with numerical examples.

preprint2013arXiv

An improved algorithm and a Fortran 90 module for computing the conical function $P^m_{-1/2+iτ}(x)$

In this paper we describe an algorithm and a Fortran 90 module ({\bf Conical}) for the computation of the conical function $P^m_{-\tfrac12+iτ}(x)$ for $x>-1$, $m \ge 0$, $τ>0$. These functions appear in the solution of Dirichlet problems for domains bounded by cones; because of this, they are involved in a large number of applications in Engineering and Physics. In the Fortran 90 module, the admissible parameter ranges for computing the conical functions in standard IEEE double precision arithmetic are restricted to $(x,m,τ) \in (-1,1) \times [0,\,40] \times [0,\,100]$ and $(x,m,τ) \in (1,100) \times [0,\,100] \times [0,\,100]$. Based on tests of the three-term recurrence relation satisfied by these functions and direct comparison with Maple, we claim a relative accuracy close to $10^{-12}$ in the full parameter range, although a mild loss of accuracy can be found at some points of the oscillatory region of the conical functions. The relative accuracy increases to $10^{-13}\,-\,10^{-14}$ in the region of the monotonic regime of the functions where integral representations are computed ($-1<x<0$).

preprint2013arXiv

Efficient and accurate algorithms for the computation and inversion of the incomplete gamma function ratios

Algorithms for the numerical evaluation of the incomplete gamma function ratios $P(a,x)=γ(a,x)/Γ(a)$ and $Q(a,x)=Γ(a,x)/Γ(a)$ are described for positive values of $a$ and $x$. Also, inversion methods are given for solving the equations $P(a,x)=p$, $Q(a,x)=q$, with $0<p,q<1$. Both the direct computation and the inversion of the incomplete gamma function ratios are used in many problems in statistics and applied probability. The analytical approach from earlier literature is summarized and new initial estimates are derived for starting the inversion algorithms. The performance of the associated software to our algorithms (the Fortran 90 module {\bf IncgamFI}) is analyzed and compared with earlier published algorithms.

preprint2004arXiv

The ABC of Hyper Recursions

Each family of Gauss hypergeometric functions $$ f_n={}_2F_1(a+ε_1n, b+ε_2n ;c+ε_3n; z), $$ for fixed $ε_j=0,\pm1$ (not all $ε_j$ equal to zero) satisfies a second order linear difference equation of the form $$ A_nf_{n-1}+B_nf_n+C_nf_{n+1}=0. $$ Because of symmetry relations and functional relations for the Gauss functions, many of the 26 cases (for different $ε_j$ values) can be transformed into each other. We give a set of basic equations from which all other equations can be obtained. For each basic equation, we study the existence of minimal solutions and the character of $f_n$ (minimal or dominant) as $n\to \pm\infty$. A second independent solution is given in each basic case which is dominant when $f_n$ is minimal and vice-versa. In this way, satisfactory pairs of linearly independent solutions for each of the 26 second order linear difference equations can be obtained.