code: def _integrand_A15(self, r, R, kwargs_mass, kwargs_light, kwargs_anisotropy):
          
          k_r = self.anisotropy.K(r, R, kwargs_anisotropy)
          l_r = self.lightProfile.light_3d_interp(r, kwargs_light)
          m_r = self.massProfile.mass_3d_interp(r, kwargs_mass)
          out = k_r * l_r * m_r / r
          return out
docs: integrand of A15 (in log space) in Mamon&Lokas 2005
  
          :param r: 3d radius
          :param R: 2d projected radius
          :param kwargs_mass: mass model parameters (following lenstronomy lens model conventions)
          :param kwargs_light: deflector light parameters (following lenstronomy light model conventions)
          :param kwargs_anisotropy: anisotropy parameters, may vary according to anisotropy type chosen.
              We refer to the Anisotropy() class for details on the parameters.
          :return: