Related papers: Solving Kepler's equation via Smale's $\alpha$-the…
We provide an approximate zero $\widetilde{S}(g,L)$ for the hyperbolic Kepler's equation $S-g\, {\rm asinh} (S)-L=0$ for $g\in(0,1)$ and $L\in[0,\infty)$. We prove, by using Smale's $\alpha$-theory, that Newton's method starting at our…
Context. Many algorithms to solve Kepler's equations require the evaluation of trigonometric or root functions. Aims. We present an algorithm to compute the eccentric anomaly and even its cosine and sine terms without usage of other…
In a recent MNRAS article, Raposo-Pulido and Pelaez (RPP) designed a scheme for obtaining very close seeds for solving the elliptic Kepler Equation with the classical and the modified Newton-Rapshon methods. This implied an important…
Denote by $E_\epsilon$ the Ginzburg-Landau functional in the plane and let $\tilde u_\varepsilon$ be the radial solution to the Euler equation associated to the problem $\min \left\{E_\varepsilon(u,B_1): \>\left. u\right\vert _{\partial…
The minimum Kullback entropy principle (mKE) is a useful tool to estimate quantum states and operations from incomplete data and prior information. In general, the solution of a mKE problem is analytically challenging and an approximate…
Numerically obtaining the inverse of a function is a common task for many scientific problems, often solved using a Newton iteration method. Here we describe an alternative scheme, based on switching variables followed by spline…
We present novel approaches for solving Milne's equation, which was introduced in 1930 as an efficient numerical scheme for the Schr\"odinger equation. Milne's equation appears in a wide class of physical problems, ranging from astrophysics…
In a previous work, we developed the idea to solve Kepler's equation with a CORDIC-like algorithm, which does not require any division, but still multiplications in each iteration. Here we overcome this major shortcoming and solve Kepler's…
The two-parameter Mittag-Leffler function $E_{\alpha, \beta}$ is of fundamental importance in fractional calculus. It appears frequently in the solutions of fractional differential and integral equations. Nonetheless, this vital function is…
We reexamine Smale's alpha theory as a way to certify a numerical solution to an analytic system. For a given point and a system, Smale's alpha theory determines whether Newton's method applied to this point shows the quadratic convergence…
A class of bivariate infinite series solutions of the elliptic and hyperbolic Kepler equations is described, adding to the handful of 1-D series that have been found throughout the centuries. This result is based on an iterative procedure…
The complete elliptic integral of the first and second kind, K(k) and E(k), appear in a multitude of physics and engineering applications. Because there is no known closed-form, the exact values have to be computed numerically. Here,…
The Bethe-Salpeter equation (BSE) for bound states in scalar theories is reformulated and solved in terms of a generalized spectral representation directly in Minkowski space. This differs from the conventional approach, where the BSE is…
We consider the dissipative spin-orbit problem in Celestial Mechanics, which describes the rotational motion of a triaxial satellite moving on a Keplerian orbit subject to tidal forcing and "drift". Our goal is to construct quasi-periodic…
In his 1981 Fundamental Theorem of Algebra paper Steve Smale initiated the complexity theory of finding a solution of polynomial equations of one complex variable by a variant of Newton's method. In this paper we reconsider his algorithm in…
A quasi-Toeplitz $M$-matrix $A$ is an infinite $M$-matrix that can be written as the sum of a semi-infinite Toeplitz matrix and a correction matrix. This paper is concerned with computing the square root of invertible quasi-Toeplitz…
We present a sublinear query algorithm for outputting a near-optimal low-rank approximation to any positive semidefinite Toeplitz matrix $T \in \mathbb{R}^{d \times d}$. In particular, for any integer rank $k \leq d$ and $\epsilon,\delta >…
In quantum theory, the so-called "spinless Salpeter equation," the relativistic generalization of the nonrelativistic Schroedinger equation, is used to describe both bound states of scalar particles and the spin-averaged spectra of bound…
Submodular functions are set functions mapping every subset of some ground set of size $n$ into the real numbers and satisfying the diminishing returns property. Submodular minimization is an important field in discrete optimization theory…
We consider autonomous stochastic ordinary differential equations (SDEs) and weak approximations of their solutions for a general class of sufficiently smooth path-dependent functionals f. Based on tools from functional It\^o calculus, such…