Beginning with a Lagrangian, we derived an approximate relativistic orbit equation which describes relativistic corrections to Keplerian orbits. The critical angular moment to guarantee the existence of periodic orbits is determined. An approximate relativistic Kepler’s elliptic orbit is illustrated by numerical simulation via a second-order perturbation method of averaging.