code: def mass_integral(self, x, axis_ratio):
          
          r = x * axis_ratio
          return 2 * np.pi * r * self.convergence_func(x)
docs: Routine to integrate an elliptical light profiles - set axis ratio to 1 to compute the luminosity within a \
          circle