Although Kepler’s laws can be empirically proven by applying Newton’s laws to the dynamics of two particles attracted by gravitational interaction, an explicit formula for the motion as a function of time remains undefined. This paper proposes a quasi-analytical solution to address this challenge. It approximates the real dynamics of celestial bodies with a satisfactory degree of accuracy and minimal computational cost. This problem is closely related to Kepler’s equation, as solving the equations of motion as a function of time also provides a solution to Kepler’s equation. The results are presented for each planet of the solar system, including Pluto, and the solution is compared against real orbits.

Quantum systems whose states are tightly distributed among several invariant subspaces (variable spin systems) can be described in terms of distributions in a four-dimensional phase-space $T^{\ast }{\mathcal{S}}_2$ in the limit of large average angular momentum. The cotangent bundle $T^{\ast }{\mathcal{S}}_2$ is also the classical manifold for systems with E(3) symmetry group with appropriately fixed Casimir operators. This allows us to employ the asymptotic form of the star-product proper for variable (integer) spin systems to develop a deformation quantization scheme for a particle moving on the two-dimensional sphere, whose observables are elements of e(3) algebra and the corresponding phase-space is $T^{\ast }{\mathcal{S}}_2$. We show that the standard commutation relations of the $e(3$) algebra are recovered from the corresponding classical Poisson brackets and the explicit expressions for the eigenvalues and eigenfunctions of some quantized classical observables (such as the angular momentum operators and their squares) are obtained.

Although the empirical Kepler’s laws can be proven by applying Newton’s laws to the dynamics of two particles attracted due to the gravitational interaction, there is no explicit formula for the motion as a function of time. In this paper, a quasi-analytical solution for this problem is proposed. It approximates the real behavior of celestial bodies with an acceptable degree of accuracy. All calculations involved require a low computational cost. This problem is closely related to Kepler’s equation, since, the solution for the equations of motion as a function of time gives us the solution to Kepler’s equation as well. The results are presented for each planet of the solar system (including Pluto), and the solution is compared against the real orbits.