Related papers: Solving Kepler's equation with CORDIC double itera…
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…
It is argued that, for motion in a central force field, polar reciprocals of trajectories are an elegant alternative to hodographs. The principal advantage of polar reciprocals is that the transformation from a trajectory to its polar…
A Kepler solver is an analytical method used to solve a two-body problem. In this paper, we propose a new correction method by slightly modifying the Kepler solver. The only change to the analytical solutions is that the obtainment of the…
We present a simple method to obtain the solution of a few orbital problems: the Kepler problem, the modified Kepler problem by the addition of an inverse square potential and linear force.
A fundamental relation in celestial mechanics is Kepler's equation, linking an orbit's mean anomaly to its eccentric anomaly and eccentricity. Being transcendental, the equation cannot be directly solved for eccentric anomaly by…
Numerical solutions of Kepler's Equation are critical components of celestial mechanics software, and are often computation hot spots. This work uses symbolic regression and a genetic learning algorithm to find new initial guesses for…
An efficient geometric integrator is proposed for solving the perturbed Kepler motion. This method is stable and accurate over long integration time, which makes it appropriate for treating problems in astrophysics, like solar system…
In computer science, transforming spherical coordinates into Cartesian coordinates is an important mathematical operation. The CORDIC (Coordinate Rotation Digital Computer) iterative algorithm can perform this operation, as well as…
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…
A pure two-body problem has seven integrals including the Kepler energy, the Laplace vector, and the angular momentum vector. However, only five of them are independent. When the five independent integrals are preserved, the two other…
This work introduces a quantum algorithm for computing the function arcsine, with arbitrary accuracy. We leverage a technique from embedded computing and Field-Programmable Gate Arrays, called COordinate Rotation DIgital Computer (CORDIC).…
In order to approximate transandental functions, several algorithms were proposed.Historically, polynomial interpolation, infinite series, $\cdots$ and other$+,\times, -$ and $/$ based algorithms were studied for this purpose.The CORDIC…
To solve the spinor-spinor Bethe-Salpeter equation in Euclidean space we propose a novel method related to the use of hyperspherical harmonics. We suggest an appropriate extension to form a new basis of spin-angular harmonics that is…
In this paper we present a framework which provides an analytical (i.e., infinitely differentiable) transformation between spatial coordinates and orbital elements for the solution of the gravitational two-body problem. The formalism omits…
The first integrals of the Kepler problem are used to compute preliminary orbits starting from two short observed arcs of a celestial body, which may be obtained either by optical or radar observations. We write polynomial equations for…
The systems of nonlinear Volterra integral equations of the first kind with jump discontinuous kernels are studied. The iterative numerical method for such nonlinear systems is proposed. Proposed method employs the modified…
A test particle in a noncoplanar orbit about a member of a binary system can undergo Kozai-Lidov oscillations in which tilt and eccentricity are exchanged. An initially circular highly inclined particle orbit can reach high eccentricity. We…
We introduce a new iterative method for computing solutions of elliptic equations with random rapidly oscillating coefficients. Similarly to a multigrid method, each step of the iteration involves different computations meant to address…
The paper suggests a preconditioning type method for fast solving of elliptic equations with oscillating quasiperiodic coefficients $A_\epsilon$ specified by the small parameter $\epsilon>0$. We use an iteration method generated by an…
We propose the numerical methods for solution of the weakly regular linear and nonlinear evolutionary (Volterra) integral equation of the first kind. The kernels of such equations have jump discontinuities along the continuous curves…