diff --git a/mapie/conformity_scores/residual_conformity_scores.py b/mapie/conformity_scores/residual_conformity_scores.py
index 584ee2284..ee8ae630a 100644
--- a/mapie/conformity_scores/residual_conformity_scores.py
+++ b/mapie/conformity_scores/residual_conformity_scores.py
@@ -804,7 +804,7 @@ def get_bounds(
             X, ensemble
         )
         signed = -1 if self.sym else 1
-        conformity_scores = conformity_scores * y_std_multi
+        conformity_scores = conformity_scores * np.maximum(self.eps, y_std_multi)
         if method == "plus":
             alpha_low = alpha_np if self.sym else alpha_np / 2
             alpha_up = 1 - alpha_np if self.sym else 1 - alpha_np / 2
diff --git a/mapie/regression/test.py b/mapie/regression/test.py
index d3b0a60c8..12b72c888 100644
--- a/mapie/regression/test.py
+++ b/mapie/regression/test.py
@@ -1,3 +1,7 @@
+import sys
+
+sys.path.append("/Users/edgarjaber/Desktop/Thèse/Conférences/CEMRACS-2023/projet/MAPIE")
+
 import numpy as np
 from sklearn.gaussian_process import GaussianProcessRegressor
 from sklearn.model_selection import train_test_split