Related papers: Computation of parabolic cylinder functions having…
Some sharp two-sided Tur\'an type inequalities for parabolic cylinder functions and Tricomi confluent hypergeometric functions are deduced. The proofs are based on integral representations for quotients of parabolic cylinder functions and…
An approximation Ansatz for the operator solution, $U(z',z)$, of a hyperbolic first-order pseudodifferential equation, $\d_z + a(z,x,D_x)$ with $\Re (a) \geq 0$, is constructed as the composition of global Fourier integral operators with…
We consider the asymptotic expansion of the Wright function \[W_{\lambda,\mu}(z)=\sum_{n=0}^\infty\frac{z^n}{n! \Gamma(\lambda n+\mu)}\qquad (\lambda>-1)\] for large (positive and negative) variable and large parameter $\mu$. The analysis…
Uniform asymptotic expansions are derived for Whittaker's confluent hypergeometric functions $M_{\kappa,\mu}(z)$ and $W_{\kappa,\mu}(z)$, as well as the numerically satisfactory companion function $W_{-\kappa,\mu}(ze^{-\pi i})$. The…
In a recent paper \cite{Temme:2021:AKH} new asymptotic expansions are given for the Kummer functions $M(a,b,z)$ and $U(a,b+1,z)$ for large positive values of $a$ and $b$, with $z$ fixed and special attention for the case $a\sim b$. In this…
A new algorithm to calculate Coulomb wave functions with all of its arguments complex is proposed. For that purpose, standard methods such as continued fractions and power/asymptotic series are combined with direct integrations of the…
In this paper we present a unified framework for asymptotic analysis and computation of the finite Hankel transform. This framework enables us to derive asymptotic expansions of the transform, including the cases where the oscillator has…
The Brent-McMillan algorithm is the fastest known procedure for the high-precision computation of Euler's constant $\gamma$ and is based on the modified Bessel functions $I_0(2x)$ and $K_0(2x)$. An error estimate for this algorithm relies…
The cylindrical algebraic covering method was originally proposed to decide the satisfiability of a set of non-linear real arithmetic constraints. We reformulate and extend the cylindrical algebraic covering method to allow for checking the…
In this work, we extend the fractional linear multistep methods in [C. Lubich, SIAM J. Math. Anal., 17 (1986), pp.704--719] to the tempered fractional integral and derivative operators in the sense that the tempered fractional derivative…
We propose a new asymptotic expansion for the fractional $p$-Laplacian with precise computations of the errors. Our approximation is shown to hold in the whole range $p\in(1,\infty)$ and $s\in(0,1)$, with errors that do not degenerate as…
A discrete rotation algorithm can be apprehended as a parametric application $f\_\alpha$ from $\ZZ[i]$ to $\ZZ[i]$, whose resulting permutation ``looks like'' the map induced by an Euclidean rotation. For this kind of algorithm, to be…
Some mathematical models of applied problems lead to the need of solving boundary value problems with a fractional power of an elliptic operator. In a number of works, approximations of such a nonlocal operator are constructed on the basis…
We consider the numerical approximation of parabolic-elliptic interface problems by the non-symmetric coupling method of MacCamy and Suri [Quart. Appl. Math., 44 (1987), pp. 675--690]. We establish well-posedness of this formulation for…
A recently developed method for the calculation of Lyapunov exponents of dynamical systems is described. The method is applicable whenever the linearized dynamics is Hamiltonian. By utilizing the exponential representation of symplectic…
Periodic orbit theory provides two important functions---the dynamical zeta function and the spectral determinant for the calculation of dynamical averages in a nonlinear system. Their cycle expansions converge rapidly when the system is…
In the article we outline the set of Matlab functions that enable the computation of elliptic Integrals and Jacobian elliptic functions for real arguments. Correctness, robustness, efficiency and accuracy of the functions are discussed in…
Integral representations are considered of solutions of the inhomogeneous Airy differential equation $w''-z w=\pm1/\pi$. The solutions of these equations are also known as Scorer functions. Certain functional relations for these functions…
Bayesian probabilistic numerical methods are a set of tools providing posterior distributions on the output of numerical methods. The use of these methods is usually motivated by the fact that they can represent our uncertainty due to…
We construct a new type of convergent, and asymptotic, representations, dyadic expansions. Their convergence is geometric and the region of convergence often extends from infinity down to $0^+$. We show that dyadic expansions are…