numerical_integral calls routines in the GNU Scientific Library. As described in the relevant introduction, the library returns a numerical integration result and an absolute error that satisfy
|RESULT - I| <= ABSERR <= max(epsabs, epsrel |I|)

where I is the exact value of the integral. The absolute and relative error inputs on the right both default to 1e-6 in SageMath.