diff --git a/include/evalhyd/detail/probabilist/intervals.hpp b/include/evalhyd/detail/probabilist/intervals.hpp
index 4e4d29aff4528104ee132715f4a45042c2eda024..3bfe578c5eff0c5f74e83add4fb3086bddcc29cb 100644
--- a/include/evalhyd/detail/probabilist/intervals.hpp
+++ b/include/evalhyd/detail/probabilist/intervals.hpp
@@ -88,9 +88,9 @@ namespace evalhyd
                     } else 
                     {
                         xt::view(itv_bnds, xt::all(), xt::all(), i, 0, xt::all()) =
-                                xt::view(q_prd, xt::all(), xt::all(), 0, xt::all());
+                                xt::view(q_prd, xt::all(), xt::all(), std::min(res[0][0], res[0][1]), xt::all());
                         xt::view(itv_bnds, xt::all(), xt::all(), i, 1, xt::all()) = 
-                                xt::view(q_prd, xt::all(), xt::all(), 1, xt::all());
+                                xt::view(q_prd, xt::all(), xt::all(), std::max(res[0][0], res[0][1]), xt::all());
                     }
                     }
                 }