diff --git a/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh b/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
index 9b7a19caeab3cd1726bc01ae8de0a5c784dfbfc4..2324200b5d62d791e425a0bcea57bfee7da45015 100644
--- a/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
+++ b/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
@@ -634,7 +634,7 @@ public:
         if (swe <= pcLowSwe_)
             return pcLowSwePcValue_ + pcDerivativeLowSw_*(swe - pcLowSwe_);
 
-        else if (swe > 1.0)
+        else if (swe > 1.0 - std::numeric_limits<Scalar>::epsilon())
             return pcDerivativeHighSwEnd_*(swe - 1.0) + entryPressure_;
 
         else
@@ -649,7 +649,7 @@ public:
         if (swe < pcLowSwe_)
             return pcDerivativeLowSw_;
 
-        else if (swe > 1.0)
+        else if (swe > 1.0 - std::numeric_limits<Scalar>::epsilon())
             return pcDerivativeHighSwEnd_;
 
         else