Kepler's orbits with corrections due to Special Relativity are explored using the Lagrangian formalism. A very simple model includes only relativistic kinetic energy by defining a Lagrangian that is consistent with both the relativistic momentum of Special Relativity and Newtonian gravity. The corresponding equations of motion are solved in a Keplerian limit, resulting in an approximate relativistic orbit equation that has the same form as that derived from General Relativity in the same limit and clearly describes three characteristics of relativistic Keplerian orbits: precession of perihelion; reduced radius of circular orbit; and increased eccentricity. The prediction for the rate of precession of perihelion is in agreement with established calculations using only Special Relativity. All three characteristics are qualitatively correct, though suppressed when compared to more accurate general-relativistic calculations. This model is improved upon by including relativistic gravitational potential energy. The resulting approximate relativistic orbit equation has the same form and symmetry as that derived using the very simple model, and more accurately describes characteristics of relativistic orbits. For example, the prediction for the rate of precession of perihelion of Mercury is one-third that derived from General Relativity. These Lagrangian formulations of the special-relativistic Kepler problem are equivalent to the familiar vector calculus formulations. In this Keplerian limit, these models are supposed to be physical based on the likeness of the equations of motion to those derived using General Relativity. The derivation of this orbit equation is approachable by undergraduate physics majors and nonspecialists whom have not had a course dedicated to relativity.

An alternative derivation of the first-order relativistic contribution to perihelic precession is presented. Orbital motion in the Schwarzschild geometry is considered in the Keplerian limit, and the orbit equation is derived for approximately elliptical motion. The method of solution makes use of coordinate transformations and the correspondence principle, rather than the standard perturbative approach. The form of the resulting orbit equation is similar to that derived from Newtonian mechanics and includes first-order corrections to Kepler's orbits due to general relativity. The associated relativistic contribution to perihelic precession agrees with established first-order results. The reduced radius for the circular orbit is in agreement to first-order with that calculated from the Schwarzschild effective potential. The method of solution is understandable by undergraduate students.

The problem of a test body in the Schwarzschild geometry is investigated in a Keplerian limit. Beginning with the Schwarzschild metric, a solution to the limited case of approximately elliptical (Keplerian) motion is derived in terms of trigonometric functions. This solution is similar in form to that derived from Newtonian mechanics, and includes first-order corrections describing three effects due to general relativity: precession; reduced radial coordinate; and increased eccentricity. The quantitative prediction of increased eccentricity may provide an additional observational test of general relativity. By analogy with Keplerian orbits, approximate orbital energy parameters are defined in terms of a relativistic eccentricity, providing first-order corrections to Newtonian energies for elliptical orbits. The first-order relativistic equation of orbit is demonstrated to be a limiting case of a very accurate self-consistent solution. This self-consistent solution is supported by exact numerical solutions to the Schwarzschild geometry, displaying remarkable agreement. A more detailed energy parameterization is investigated using the relativistic eccentricity together with the apsides derived from the relativistic effective potential in support of the approximate energy parameters defined using only first-order corrections. The methods and approximations describing this Keplerian limit are applied to more general static spherically-symmetric geometries. Specifically, equations of orbit and energy parameters are also derived in this Keplerian limit for the Reissner-Nordstr\"{o}m and Schwarzschild-de Sitter metrics.