diff --git a/apps/probability/distribution/distribution.cpp b/apps/probability/distribution/distribution.cpp
index fafdcdad4..549c36177 100644
--- a/apps/probability/distribution/distribution.cpp
+++ b/apps/probability/distribution/distribution.cpp
@@ -105,10 +105,10 @@ double Distribution::evaluateAtDiscreteAbscissa(int k) const {
 
 double Distribution::cumulativeDistributiveInverseForProbabilityUsingIncreasingFunctionRoot(double * probability, double ax, double bx) {
   assert(ax < bx);
-  if (*probability >= 1.0) {
+  if (*probability > 1.0 - DBL_EPSILON) {
     return INFINITY;
   }
-  if (*probability <= 0.0) {
+  if (*probability < DBL_EPSILON) {
     return -INFINITY;
   }
   Poincare::Coordinate2D<double> result = Poincare::Solver::IncreasingFunctionRoot(