Adding new equations of motion accounting for multiple bodies, gravitating and luminous. 3D kepler example in geostationary orbit works.