diff --git a/FlowCal/plot.py b/FlowCal/plot.py
index 4236fd5..5698ca0 100644
--- a/FlowCal/plot.py
+++ b/FlowCal/plot.py
@@ -1780,7 +1780,8 @@ def density_and_hist(data,
 
     # Colors
     n_colors = n_plots - 1
-    colors = [cmap_default(i) for i in np.linspace(0, 1, n_colors)]
+    default_property_cycler = plt.rcParams['axes.prop_cycle']()
+    colors = [next(default_property_cycler)['color'] for i in range(n_colors)]
     # Histogram
     for i, hist_channel in enumerate(hist_channels):
         # Define subplot