| Hosted by CoCalc | Download
Kernel: SageMath (Development, Py3)
# SageMath version version()
'SageMath version 8.8.beta6, Release Date: 2019-05-22'
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)
(0.6116545934235456, 6.790730127365591e-15)