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.