code: def tabulate_integral(self, grid, tabulate_bins):
          
          eta_min = 1.0e-4
          eta_max = 1.05 * np.max(self.grid_to_elliptical_radii(grid))
  
          minimum_log_eta = np.log10(eta_min)
          maximum_log_eta = np.log10(eta_max)
          bin_size = (maximum_log_eta - minimum_log_eta) / (tabulate_bins - 1)
  
          return eta_min, eta_max, minimum_log_eta, maximum_log_eta, bin_size
docs: Tabulate an integral over the surface density of deflection potential of a mass profile. This is used in \
          the GeneralizedNFW profile classes to speed up the integration procedure.
  
          Parameters
          -----------
          grid : grids.RegularGrid
              The grid of (y,x) arc-second coordinates the potential / deflection_stacks are computed on.
          tabulate_bins : int
              The number of bins to tabulate the inner integral of this profile.