diff --git a/R/custom_function.R b/R/custom_function.R
index 66c2256..704ecf4 100644
--- a/R/custom_function.R
+++ b/R/custom_function.R
@@ -59,11 +59,17 @@ distribution = function(value, coverage_fraction, breaks = 10, constrain = FALSE
   x1 = value*coverage_fraction
   x1 = x1[!is.na(x1)]
 
-  if(constrain | length(breaks) > 1){
+  if(constrain & length(breaks) > 1){
 
     breaks_tmp = c(breaks[1],breaks[2])
 
-    breaks = breaks[!breaks > max(x1, na.rm = TRUE)]
+    ulimit = max(x1, na.rm = TRUE)
+    
+    if (ulimit < max(breaks, na.rm = TRUE)){
+       ulimit = min(breaks[breaks >= ulimit])
+    }
+
+    breaks = breaks[breaks <= ulimit]
 
     if (length(breaks) == 1){
       breaks = breaks_tmp