diff --git a/apps/probability/calculation/left_integral_calculation.cpp b/apps/probability/calculation/left_integral_calculation.cpp
index 4dd9056a6..f0ccc8198 100644
--- a/apps/probability/calculation/left_integral_calculation.cpp
+++ b/apps/probability/calculation/left_integral_calculation.cpp
@@ -47,6 +47,10 @@ void LeftIntegralCalculation::compute(int indexKnownElement) {
   if (indexKnownElement == 0) {
     m_result = m_law->cumulativeDistributiveFunctionAtAbscissa(m_upperBound);
   } else {
+    double currentResult = m_law->cumulativeDistributiveFunctionAtAbscissa(m_upperBound);
+    if (std::fabs(currentResult - m_result) < std::pow(10.0, - Constant::LargeNumberOfSignificantDigits)) {
+      return;
+    }
     m_upperBound = m_law->cumulativeDistributiveInverseForProbability(&m_result);
   }
 }
diff --git a/apps/probability/calculation/right_integral_calculation.cpp b/apps/probability/calculation/right_integral_calculation.cpp
index 48858f584..83a82436f 100644
--- a/apps/probability/calculation/right_integral_calculation.cpp
+++ b/apps/probability/calculation/right_integral_calculation.cpp
@@ -47,6 +47,10 @@ void RightIntegralCalculation::compute(int indexKnownElement) {
   if (indexKnownElement == 0) {
     m_result = m_law->rightIntegralFromAbscissa(m_lowerBound);
   } else {
+    double currentResult = m_law->rightIntegralFromAbscissa(m_lowerBound);
+    if (std::fabs(currentResult - m_result) < std::pow(10.0, - Constant::LargeNumberOfSignificantDigits)) {
+      return;
+    }
     m_lowerBound = m_law->rightIntegralInverseForProbability(&m_result);
   }
 }