First time here? Check out the FAQ!
x = var('x') a = 5.0 + I b = 3.0 + I numerical_integral((1/sqrt(x^3 + a*x + b)).real(), 0, 3)