Calculation of definite integrals implies to replace the function by another one (usually a polynomial form) which is a good approximation and easier to compute.
Error estimation depends on parameters of the method and refinement of spatial grid discretization. It is interesting to check these errors using various algorithm and mesh resolution.