diff --git a/include/evalhyd/detail/probabilist/intervals.hpp b/include/evalhyd/detail/probabilist/intervals.hpp
index 51ab71ec8f2953e90a0ea0707ebdd374c5576736..2052631946b8bbfb2150fe9823ba593ddd2c8c40 100644
--- a/include/evalhyd/detail/probabilist/intervals.hpp
+++ b/include/evalhyd/detail/probabilist/intervals.hpp
@@ -67,11 +67,11 @@ namespace evalhyd
                     for (std::size_t i = 0; i < n_itv; i++)
                     {
                     auto a = xt::broadcast(xt::view(quantiles, i), std::vector<std::size_t>({q_lvl.size(), 2}));
-                    auto b = xt::broadcast(q_lvl, std::vector<std::size_t>({2, q_lvl.size()}));
+                    auto b = xt::broadcast(q_lvl / 100., std::vector<std::size_t>({2, q_lvl.size()}));
                     auto res = xt::where(xt::equal(a, xt::transpose(b)));
                     std::cout << "res: " << res.size() << std::endl;
-                    std::cout << "res: " << res[0][1] << std::endl;
-                    std::cout << "res: " << res[1][1] << std::endl;
+                    std::cout << "res: " << res[0] << std::endl;
+                    std::cout << "res: " << res[1] << std::endl;
                     if (res.size() != 2) 
                     {
                         throw std::runtime_error(