# Issues with numerically integrating a complex value function. Any way of getting around the 'Unable to convert 5.0 + 1.0*I to float; use abs() or real_part() as desired' error?

x = var('x')
a = 5.0 + I

b = 3.0 + I

numerical_integral((1/sqrt(x^3 + a*x + b)).real(), 0, 3)