Related papers: An asymptotic expansion for the error term in the …
Brent and McMillan introduced in 1980 a new algorithm for the computation of Euler's constant $\gamma$, based on the use of the Bessel functions I\_0(x) and K\_0(x). It is the fastest known algorithm for the computation of $\gamma$. The…
The Brent-McMillan algorithm B3 (1980), when implemented with binary splitting, is the fastest known algorithm for high-precision computation of Euler's constant. However, no rigorous error bound for the algorithm has ever been published.…
A fast and numerically stable algorithm is described for computing the discrete Hankel transform of order $0$ as well as evaluating Schl\"{o}milch and Fourier--Bessel expansions in $\mathcal{O}(N(\log N)^2/\log\!\log N)$ operations. The…
In this paper, we reconsider the large-argument asymptotic expansions of the Hankel, Bessel and modified Bessel functions and their derivatives. New integral representations for the remainder terms of these asymptotic expansions are found…
We apply the Euler--Maclaurin formula to find the asymptotic expansion of the sums $\sum_{k=1}^n (\log k)^p / k^q$, ~$\sum k^q (\log k)^p$, ~$\sum (\log k)^p /(n-k)^q$, ~$\sum 1/k^q (\log k)^p $ in closed form to arbitrary order ($p,q…
We design an energy-stable and asymptotic-preserving finite volume scheme for the compressible Euler system. Using the relative energy framework, we establish rigorous error estimates that yield convergence of the numerical solutions in two…
A recent asymptotic expansion for the positive zeros $x=j_{\nu,m}$ ($m=1,2,3,\ldots$) of the Bessel function of the first kind $J_{\nu}(x)$ is studied, where the order $\nu$ is positive. Unlike previous well-known expansions in the…
We examine the sum of modified Bessel functions with argument depending non-linearly on the summation index given by \[S_{\nu,p}(a)=\sum_{n\geq 1} (an^p/2)^{-\nu} K_\nu(an^p)\qquad (a>0,\ 0\leq\nu<1)\] as the parameter $a\to 0+$, where $p$…
By application of the theory for second-order linear differential equations with two turning points developed in [Olver F.W.J., Philos. Trans. Roy. Soc. London Ser. A 278 (1975), 137-174], uniform asymptotic approximations are obtained in…
Uniform asymptotic expansions are derived for the zeros of the reverse generalized Bessel polynomials of large degree $n$ and real parameter $a$. It is assumed that $-\Delta_{1} n+\frac{3}{2} \leq a \leq \Delta_{2} n$ for fixed arbitrary…
The paper presents the derivation of the asymptotic behavior of $\nu$-zeros of the modified Bessel function of imaginary order $K_{{\rm i}\nu}(z)$. This derivation is based on the quasiclassical treatment of the exponential potential on the…
The asymptotic expansion of digamma function is a starting point for the derivation of approximants for harmonic sums or Euler-Mascheroni constant. It is usual to derive such approximations as values of logarithmic function, which leads to…
Bessel and modified Bessel functions of imaginary order $i\nu$ ($\nu >0$) are studied. Asymptotic expansions are derived as $\nu \to \infty$ that are uniformly valid in unbounded complex domains, with error bounds provided. Coupled with…
Over the last several decades, improvements in the fields of analytic combinatorics and computer algebra have made determining the asymptotic behaviour of sequences satisfying linear recurrence relations with polynomial coefficients largely…
We study the asymptotic expansion for the Landau constants $G_n$ $$\pi G_n\sim \ln N + \gamma+4\ln 2 + \sum_{s=1}^\infty \frac {\beta_{2s}}{N^{2s}},~~n\rightarrow \infty, $$ where $N=n+3/4$, $\gamma=0.5772\cdots$ is Euler's constant, and…
We give an asymptotic development of the maximum likelihood estimator (MLE), or any other estimator defined implicitly, in a way which involves the limiting behavior of the score and its higher-order derivatives. This development, which is…
Computable and sharp error bounds are derived for asymptotic expansions for linear differential equations having a simple turning point. The expansions involve Airy functions and slowly varying coefficient functions. The sharpness of the…
Reformulated uniform asymptotic expansions are derived for ordinary differential equations having a large parameter and a simple turning point. These involve Airy functions, but not their derivatives, unlike traditional asymptotic…
We consider the uniform asymptotic expansion for the Gauss hypergeometric function \[F(a+\epsilon\lambda,m;c+\lambda;x),\qquad \lambda\to+\infty\] for $x<1$ and positive integer $m$ when the parameter $\epsilon>1$ and the constants $a$ and…
The Conway-Maxwell-Poisson distribution is a two-parameter generalisation of the Poisson distribution that can be used to model data that is under- or over-dispersed relative to the Poisson distribution. The normalizing constant…