x_numerator = math.sqrt(1.0 - 2.0 * rho * z + z * z) + z - rho