diff --git a/plot_scripts/contraction_benefits_whiskers.py b/plot_scripts/contraction_benefits_whiskers.py
index 4934cf95452d47194a3fbba10ff5101000bbe55f..77af0a94568aa6e1e6dcc3887612a413948f8378 100755
--- a/plot_scripts/contraction_benefits_whiskers.py
+++ b/plot_scripts/contraction_benefits_whiskers.py
@@ -40,13 +40,13 @@ for i in [0, 1, 2]:
     axs[i].axhline(y=60, color='lightblue', linestyle='--')
     axs[i].axhline(y=80, color='lightblue', linestyle='--')
 
-bplot1 = axs[0].boxplot(vars_datas, showfliers=False, patch_artist=True)
+bplot1 = axs[0].boxplot(vars_datas, meanline=True, showmeans=True, whis=(0, 100), patch_artist=True)
 axs[0].set(xlabel=None, ylabel='removed variables')
 axs[0].set_xticklabels([])
-bplot2 = axs[1].boxplot(constr_datas, showfliers=False, patch_artist=True)
+bplot2 = axs[1].boxplot(constr_datas, meanline=True, showmeans=True, whis=(0, 100), patch_artist=True)
 axs[1].set(xlabel=None, ylabel='removed constraints')
 axs[1].set_xticklabels([])
-bplot3 = axs[2].boxplot(elems_datas, showfliers=False, patch_artist=True)
+bplot3 = axs[2].boxplot(elems_datas, meanline=True, showmeans=True, whis=(0, 100), patch_artist=True)
 axs[2].set(xlabel="percentage of improvable arcs", ylabel='removed entries')
 axs[2].set_xticklabels([x if x % 10 == 0 else "" for x in x_datas])
 
@@ -54,8 +54,6 @@ for bplot in (bplot1, bplot2, bplot3):
     for patch in bplot['boxes']:
         patch.set_facecolor('white')
 
-# legend = plt.legend(loc='lower right', shadow=True, fontsize='medium')
-
 plt.tight_layout()
 fig.subplots_adjust(hspace=0.08)