calc_zeta_inf

luminet.black_hole_math.calc_zeta_inf(p: float, bh_mass: float) → float

Calculate \(\zeta_\infty\)

This is used in the Jacobi incomplete elliptic integral \(F(\zeta_\infty, k)\)

\[\zeta_\infty = \arcsin \left( \sqrt{\frac{Q - P + 2M}{Q - P + 6M}} \right)\]

Parameters:

• p (float) – periastron distance

• bh_mass (float) – Black hole mass

Returns:

float – \(\zeta_\infty\)