stable_trajectory

kerrgeopy.stable.stable_trajectory(a, p, e, x, initial_phases=(0, 0, 0, 0), M=None, distance_units='natural', time_units='natural')

Computes the time, radial, polar, and azimuthal coordinates of the orbit with the given parameters as a function of mino time.

Parameters:

• a (double) – dimensionless spin parameter

• p (double) – semi-latus rectum

• e (double) – orbital eccentricity

• x (double) – cosine of the orbital inclination

• initial_phases (tuple, optional) – tuple of initial phases \((q_{t_0},q_{r_0},q_{\theta_0},q_{\phi_0})\), defaults to (0,0,0,0)

• M (double, optional) – mass of the primary in solar masses

• distance_units (str, optional) – units to compute the radial component of the trajectory in (options are “natural”, “mks”, “cgs”, “au” and “km”), defaults to “natural”

• time_units (str, optional) – units to compute the time component of the trajectory in (options are “natural”, “mks”, “cgs”, and “days”), defaults to “natural”

Returns:

tuple of functions \((t(\lambda), r(\lambda), \theta(\lambda), \phi(\lambda))\)

Return type:

tuple(function, function, function, function)