code: def sigma_r2(self, r, kwargs_profile, kwargs_anisotropy, kwargs_light):
          
          if self._mass_profile == 'power_law':
              if self._anisotropy_type == 'r_ani':
                  if self._light_profile == 'Hernquist':
                      sigma_r = self.power_law_anisotropy(r, kwargs_profile, kwargs_anisotropy, kwargs_light)
                  else:
                      raise ValueError('light profile %s not supported for Jeans solver' % self._light_profile)
              else:
                  raise ValueError('anisotropy type %s not implemented in Jeans equation modelling' % self._anisotropy_type)
          else:
              raise ValueError('mass profile type %s not implemented in Jeans solver' % self._mass_profile)
          return sigma_r
docs: solves radial Jeans equation