diff --git a/matlab/set_prior.m b/matlab/set_prior.m
index cbbd6cc72..3d3aa458f 100644
--- a/matlab/set_prior.m
+++ b/matlab/set_prior.m
@@ -192,14 +192,12 @@ end
 % truncation parameters by default for normal distribution
 k  = find(bayestopt_.pshape == 3);
 k1 = find(isnan(bayestopt_.p3(k)));
+k2 = find(isnan(bayestopt_.p4(k)));
 bayestopt_.p3(k(k1)) = -Inf*ones(length(k1),1);
-k1 = find(isnan(bayestopt_.p4(k)));
-bayestopt_.p4(k(k1)) = Inf*ones(length(k1),1);
-for i=1:length(k)
-    bayestopt_.p6(k(i)) = bayestopt_.p1(k(i)) ; 
-    bayestopt_.p7(k(i)) = bayestopt_.p2(k(i)) ;
-    bayestopt_.p5(k(i)) = bayestopt_.p1(k(i)) ;
-end
+bayestopt_.p4(k(k2)) = Inf*ones(length(k2),1);
+bayestopt_.p6(k) = bayestopt_.p1(k);
+bayestopt_.p7(k) = bayestopt_.p2(k);
+bayestopt_.p5(k) = bayestopt_.p1(k);
 
 % inverse gamma distribution (type 1)
 k = find(bayestopt_.pshape == 4);