Related papers: Sinc integrals and tiny numbers
The value of a highly oscillatory integral is typically determined asymptotically by the behaviour of the integrand near a small number of critical points. These include the endpoints of the integration domain and the so-called stationary…
Classical multiple zeta values can be viewed as iterated integrals of the differentials $\frac{dt}{t}, \frac{dt}{1-t}$ from $0$ to $1$. In this paper, we reprove Brown's theorem: For $a_i, b_i, c_{ij}\in \mathbb{Z}$, the iterated integral…
We obtain integral representations of the $n$-th derivatives of the Bessel functions with respect to the order. The numerical evaluation of these expressions is very efficient using a double exponential integration strategy. Also, from the…
The hyperharmonic numbers h_{n}^{(r)} are defined by means of the classical harmonic numbers. We show that the Euler-type sums with hyperharmonic numbers: {\sigma}(r,m)=\sum_{n=1}^{\infty}((h_{n}^{(r)})/(n^{m})) can be expressed in terms of…
We give an example of a formula involving the sinc function that holds for every N = 0, 1, 2, ..., up to about 10^102832732165, then fails for all larger N. We give another example that begins to fail after about N ~…
We consider the numerical integration of the matrix Hill's equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, the Hill's equations originate from a Hamiltonian…
We consider the computation of Bernoulli, Tangent (zag), and Secant (zig or Euler) numbers. In particular, we give asymptotically fast algorithms for computing the first n such numbers in O(n^2.(log n)^(2+o(1))) bit-operations. We also give…
We consider the sum $\sum 1/\gamma$, where $\gamma$ ranges over the ordinates of nontrivial zeros of the Riemann zeta-function in an interval $(0,T]$, and consider the behaviour of the sum as $T \to\infty$. We show that, after subtracting a…
A numerical approach to compute tensor integrals in one-loop calculations is presented. The algorithm is based on a recursion relation which allows to express high rank tensor integrals as a function of lower rank ones. At each level of…
Maclaurin Integration is a new series-based technique for solving infamously difficult integrals in terms of elementary functions. It has fairly liberal conditions for sound use, making it one of the most versatile integration techniques.…
We present a new form of the Machin-like formula for $\pi$ that can be generated by using iteration. This form of the Machin-like formula may be promising for computation of the constant $\pi$ due to rapidly increasing integers at each step…
From a global series for the alternating zeta function, we derive an infinite product for pi that resembles the product for $e^\gamma$ ($\gamma$ is Euler's constant) in math.CA/0306008. (An alternate derivation accelerates Wallis's product…
We prove an Euler-Maclaurin formula for double polygonal sums and, as a corollary, we obtain approximate quadrature formulas for integrals of smooth functions over polygons with integer vertices. Our Euler-Maclaurin formula is in the spirit…
We derive an analytic expression for the scalar one-loop pentagon and hexagon functions which is convenient for subsequent numerical integration. These functions are of relevance in the computation of next-to-leading order radiative…
Integrals involving the kernel function $sech (\pi x)$ over a semi-infinite range are of general interest in the study of Riemann's function $\zeta(s)$ and Hurwitz' function $\zeta(s,a)$. Such integrals that include the $arctan$ and $log$…
We develop two classes of composite moment-free numerical quadratures for computing highly oscillatory integrals having integrable singularities and stationary points. The first class of the quadrature rules has a polynomial order of…
We introduce the Wick integral $\int_s^t p(X_u) \Diamond \mathrm{d} X_u$ for a class of stochastic processes $X$ which are not necessarily Gaussian, in the regime of bounded $2> q$-variation. The integral is defined for polynomial…
A general framework of Numerical Singular Integrals (NSI) method based on the Integration By Parts (IBP) has been developed for integrals involving singular and nearly singular integrands, or NSI-IBP. Through a general integration by parts…
Among the several proofs known for $\sum_{n=1}^\infty{1/n^2} = {\pi^2/6}$, the one by Beukers, Calabi, and Kolk involves the evaluation of $\,\int_0^1 {\int_0^1{1/(1-x^2 y^2) \, dx} \, dy}$. It starts by showing that this double integral is…
An algorithm is devised for computing $\Phi(n) = \phi(1) + \phi(2) + \cdots + \phi(n)$ in time $\widetilde{\Theta}(n^{2/3})$ and space $\widetilde{\Theta}(n^{1/3})$. The starting point is an existing algorithm based on the Dirichlet…