numerator2 = math.sqrt(1.0 - 2 * rho * z + z * z) + z - rho