I am trying to solve a simple integral in Python using different methods for numerical integration. However, when I try to calculate the integrals, I obtain wrong results for small values of my x coordinate. This is a messy because I am trying to solve some system of equations in which one of the variables I am evolving on time is the integral of some quantities. Of course, I don’t have a functional form for the expression I want to integrate but instead I have an array of data (one f[i] for each x[i]). Do you know a way to obtain more accurate results of my integrals?