code: def sigma_r2(self, r, kwargs_profile, kwargs_anisotropy, kwargs_light):
          
          return self.jeans_solver.sigma_r2(r, kwargs_profile, kwargs_anisotropy, kwargs_light)
docs: computes radial velocity dispersion at radius r (solving the Jeans equation
          :param r:
          :return: