code: def luminosity_integral(self, x, axis_ratio):
          
          r = x * axis_ratio
          return 2 * np.pi * r * self.intensities_from_grid_radii(x)
docs: Routine to integrate the luminosity of an elliptical light profile.
  
          The axis ratio is set to 1.0 for computing the luminosity within a circle