docs: def Integrate(self):
          """
          Integrate using the trapazoidal method
          """
          area = 0.
          X = self.GetX()
          Y = self.GetY()
          for i in range(self.GetN() - 1):
              area += (X[i + 1] - X[i]) * (Y[i] + Y[i + 1]) / 2.
          return area