code: def cumulative_distribution(self, X):
          
          self.check_fit()
  
          low_bounds = self.model.dataset.mean() - (5 * self.model.dataset.std())
  
          result = []
          for value in X:
              result.append(self.model.integrate_box_1d(low_bounds, value))
  
          return np.array(result)
docs: Computes the integral of a 1-D pdf between two bounds
  
          Args:
              X(numpy.array): Shaped (1, n), containing the datapoints.
  
          Returns:
              numpy.array: estimated cumulative distribution.