The OsculatingOrbitalElements package solves the forced geodesic motion in either Schwarzshild or Kerr spacetime using the Method of Osculating Orbital Elements (also known as the Method of Osculating Geodeiscs).
The user can specify the function for the force or make use of the following models included with the package:
SchwarzFastGSF[p, e, ξ]: Fast first order Gravitational Self Force Model for p<12 and e<0.2.
SchwarzGasDrag[p, e, ξ]: Relativistic gas drag model for Schwarzschild spacetime.
KerrGasDrag[En, L, K, ψr, ψθ]: Relativistic gas drag model for Kerr spacetime.
The latest development version will always be available from the project git repository:
git clone https://github.com/ucd-relativity/OsculatingOrbitalElements.git
Mathematica: OsculatingOrbitalElements requires a recent version of Mathematica. It is typically tested with only the latest available version.
Clone the repository and place it somewhere on Mathematica's $Path.
Typical locations are inside ${HOME}/.Mathematica/Applications/ for Linux,
${HOME}/Library/Mathematica/Applications/ for Mac OSX, or in
C:\Users\yourusername\AppData\Roaming\Mathematica\Applications
for Windows.
This package requires the KerrGeodesics Package.
This package, along with the KerrGeodesics package may be loaded into Mathematica using the command:
<< OsculatingOrbitalElements`Examples are included in tutorial.nb.
25 March 2019: Initial version created.
19 September 2019: Added code for Kerr Spacetime
Known bugs are recorded in the project bug tracker:
https://github.com/ucd-relativity/OsculatingOrbitalElements/issues
This code is distributed under the MIT Open Source License. Details can be found in the LICENSE file.
Philip Lynch