diff --git a/matlab/optimization/gmhmaxlik_core.m b/matlab/optimization/gmhmaxlik_core.m
index cf5fb7774b9ba676fd338cf584a1bd6b87cdace9..805516f1d1c14ae37ec38496e4afc97f8813f78e 100644
--- a/matlab/optimization/gmhmaxlik_core.m
+++ b/matlab/optimization/gmhmaxlik_core.m
@@ -149,6 +149,7 @@ j = 1;
 isux = 0;
 ilogpo2 = - feval(ObjFun,ix2,varargin{:});
 while j<= NumberOfIterations
+    j = j+1;
     proposal = iScale*dd*randn(npar,1) + ix2;
     if all(proposal > mh_bounds(:,1)) && all(proposal < mh_bounds(:,2))
         logpo2 = - feval(ObjFun,proposal,varargin{:});
@@ -175,7 +176,6 @@ while j<= NumberOfIterations
     MeanPar = oldMeanPar + (1/j)*(ix2-oldMeanPar);
     CovJump = CovJump + oldMeanPar*oldMeanPar' - MeanPar*MeanPar' + ...
               (1/j)*(ix2*ix2' - CovJump - oldMeanPar*oldMeanPar');
-    j = j+1;
 end
 dyn_waitbar_close(hh);
 PostVar = CovJump;