diff --git a/getdist/densities.py b/getdist/densities.py
index 551e41b..35ad9da 100644
--- a/getdist/densities.py
+++ b/getdist/densities.py
@@ -287,7 +287,7 @@ def Prob(self, x, y, grid=False):
         :param y: y value or array
         :param grid: whether to make a grid, see :class:`~scipy:scipy.interpolate.RectBivariateSpline`. Default False.
         """
-        self.__call__(x, y, grid=grid)
+        return self.__call__(x, y, grid=grid)
 
     def __call__(self, *args, **kwargs):
         if self.spl is None: