Related papers: Numerical methods and arbitrary-precision computat…
We consider a unique continuation problem where the Dirichlet trace of the solution is known to have finite dimension. We prove Lipschitz stability of the unique continuation problem and design a finite element method that exploits the…
Laguerre and Laguerre-type polynomials are orthogonal polynomials on the interval $[0,\infty)$ with respect to a weight function of the form $w(x) = x^{\alpha} e^{-Q(x)}, Q(x) = \sum_{k=0}^m q_k x^k, \alpha > -1, q_m > 0$. The classical…
We describe some "unrestricted" algorithms which are useful for the computation of elementary and special functions when the precision required is not known in advance. Several general classes of algorithms are identified and illustrated by…
We present an efficient numerical scheme based on Monte Carlo integration to approximate statistical solutions of the incompressible Euler equations. The scheme is based on finite volume methods, which provide a more flexible framework than…
We find two convergent series expansions for Legendre's first incomplete elliptic integral $F(\lambda,k)$ in terms of recursively computed elementary functions. Both expansions are valid at every point of the unit square $0<\lambda,k<1$.…
We develop new algorithms for approximating extremal toric K\"ahler metrics. We focus on an extremal metric on $\mathbb{CP}^{2}\sharp2\overline{\mathbb{CP}}^{2}$, which is conformal to an Einstein metric (the Chen-LeBrun-Weber metric). We…
We propose an empirical Bayes estimator based on Dirichlet process mixture model for estimating the sparse normalized mean difference, which could be directly applied to the high dimensional linear classification. In theory, we build a…
Partial differential equations with highly oscillatory input terms are hardly ever solvable analytically and their numerical treatment is difficult. Modulated Fourier expansion used as an {\it ansatz} is a well known and extensively…
We study the expansions of the completed Riemann zeta function and completed Dirichlet L-functions in Meixner-Pollaczek polynomials. We give the proof of the uniform convergence, the multiplicative structure for the coefficients of these…
By use of the threshold expansion we develop an algorithm for analytical evaluation, within dimensional regularization, of arbitrary terms in the expansion of the (two-loop) sunset diagram with general masses m_1, m_2 and m_3 near its…
For a given transcendental number $\xi$ and for any polynomial $P(X)=: \lambda_0+\cdots+\lambda_k X^k \in \mathbb{Z}[X]$, we know that $ P(\xi) \neq 0.$ Let $k \geq 1$ and $\omega (k, H)$ be the infimum of the numbers $r > 0$ satisfying the…
To decide upon the arithmetic nature of some numbers may be a non-trivial problem. Some cases are well know, for example exp(1) and W(1), where W is the Lambert function, are transcendental numbers. The Tsallis q-exponential, e_q (z), and…
We obtain the asymptotic expansion for the Gauss hypergeometric function \[F(a-\lambda,b+\lambda;c+i\alpha\lambda;z)\] for $\lambda\rightarrow+\infty$ with $a$, $b$ and $c$ finite parameters by application of the method of steepest…
In a recent work [JNT \textbf{129}, 2154 (2009)], Gun and co-workers have claimed that the number $\,\log{\Gamma(x)} + \log{\Gamma(1-x)}\,$, $x$ being a rational number between $0$ and $1$, is transcendental with at most \emph{one} possible…
We derive explicit formulas for the matroidal mixed Eulerian numbers. We resolve a question posed by Berget, Spink, and Tseng, demonstrating that the invariant defined by matroidal mixed Eulerian numbers is precisely equivalent to Derksen's…
In this paper we consider the strong asymptotic behavior of Laguerre polynomials in the complex plane. The leading behavior is well known from Perron and Mehler-Heine formulas, but higher order coefficients, which are important in the…
Loop calculations involve the evaluation of divergent integrals. Usually [1] one computes them in a number of dimensions different than four where the integral is convergent and then one performs the analytical continuation and considers…
In this project, I examine the lowest Dirichlet eigenvalue of the Laplacian within the ellipse as a function of eccentricity. Two existing analytic expansions of the eigenvalue are extended: Close to the circle (eccentricity near zero) nine…
We introduce a numerical method for the approximation of functions which are analytic on compact intervals, except at the endpoints. This method is based on variable transforms using particular parametrized exponential and…
We consider the application of multilevel Monte Carlo methods to elliptic PDEs with random coefficients. We focus on models of the random coefficient that lack uniform ellipticity and boundedness with respect to the random parameter, and…