No astrodynamics libraries are used. No astropy, no poliastro. T' = (3Tx − 2 + 2λ³x/y) / (1 − x²) T'' = (3T + 5xT' + 2(1−λ²)λ³/y³) / (1 − x²) T ...