diff --git a/src/skmatter/decomposition/_kernel_pcovr.py b/src/skmatter/decomposition/_kernel_pcovr.py
index 520e542ae..68a088770 100644
--- a/src/skmatter/decomposition/_kernel_pcovr.py
+++ b/src/skmatter/decomposition/_kernel_pcovr.py
@@ -159,7 +159,7 @@ class KernelPCovR(_BasePCA, LinearModel):
            [-0.18992219,  0.82064368],
            [ 1.11923584, -1.04798016],
            [-1.5635827 ,  1.11078662]])
-    >>> round(kpcovr.score(X, Y), 5)
+    >>> print(round(kpcovr.score(X, Y), 5))
     -0.52039
     """  # NoQa: E501
 
diff --git a/src/skmatter/neighbors/_sparsekde.py b/src/skmatter/neighbors/_sparsekde.py
index 7220eb8bb..a4f040e38 100644
--- a/src/skmatter/neighbors/_sparsekde.py
+++ b/src/skmatter/neighbors/_sparsekde.py
@@ -100,7 +100,7 @@ class SparseKDE(BaseEstimator):
 
     The total log-likelihood under the model
 
-    >>> round(estimator.score(result), 3)
+    >>> print(round(estimator.score(result), 3))
     -759.831
     """
 
diff --git a/src/skmatter/utils/_sparsekde.py b/src/skmatter/utils/_sparsekde.py
index 077131312..483fa29fe 100644
--- a/src/skmatter/utils/_sparsekde.py
+++ b/src/skmatter/utils/_sparsekde.py
@@ -22,7 +22,7 @@ def effdim(cov):
     >>> import numpy as np
     >>> from skmatter.utils import effdim
     >>> cov = np.array([[1, 1, 0], [1, 1, 0], [0, 0, 1]])
-    >>> round(effdim(cov), 3)
+    >>> print(round(effdim(cov), 3))
     1.89
 
     References