def tabulate_integral(self, grid, tabulate_bins):
eta_min = <NUM_LIT><EOL>eta_max = <NUM_LIT> * np.max(self.grid_to_elliptical_radii(grid))<EOL>minimum_log_eta = np.log10(eta_min)<EOL>maximum_log_eta = np.log10(eta_max)<EOL>bin_size = (maximum_log_eta - minimum_log_eta) / (tabulate_bins - <NUM_LIT:1>)<EOL>return eta_min, eta_max, minimum_log_eta, maximum_log_eta, bin_size<EOL>
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.