# SageMath version version()
y, z = var('y z') cauchy = y.substitute(solve(z*y^3 + y^2 - 2*z*y + 2, y)[0]) def F(u): return arg(cauchy).substitute(z=u).n() numerical_integral(F, 1, 2)