diff --git a/matlab/estimation/draw_prior_density.m b/matlab/estimation/draw_prior_density.m
index 348231a225a468966aca05aacb7b47f749d7b1a0..f835a1a71d7a2e6c69048fe6adb541896d039ca4 100644
--- a/matlab/estimation/draw_prior_density.m
+++ b/matlab/estimation/draw_prior_density.m
@@ -86,10 +86,9 @@ end
 if pshape(indx) ~= 5
     [~,k1] = max(dens);
     if k1 == 1 || k1 == length(dens)
-        k = find(dens > 10);
-        dens(k) = NaN;
-        end
+        dens(dens > 10) = NaN;
     end
+end
 binf = abscissa(1);
 bsup = abscissa(end);
 x = abscissa;