diff --git a/matlab/gsa/mc_moments.m b/matlab/gsa/mc_moments.m
index 4b06241453edea86dddd7f702448761652bf39ce..5c2cd034abc202e3d72218457c88391e069b2982 100644
--- a/matlab/gsa/mc_moments.m
+++ b/matlab/gsa/mc_moments.m
@@ -24,8 +24,8 @@ global options_ M_ estim_params_ oo_
   disp('Computing theoretical moments ...')
   h = dyn_waitbar(0,'Theoretical moments ...');
   vdec = zeros(nobs,M_.exo_nbr,nsam);
-  cc = zeros(size(options_.varobs,1),size(options_.varobs,1),nsam);
-  ac = zeros(size(options_.varobs,1),size(options_.varobs,1)*options_.ar,nsam);
+  cc = zeros(nobs,nobs,nsam);
+  ac = zeros(nobs,nobs*options_.ar,nsam);
   
   for j=1:nsam,
     oo_.dr.ghx = mm(:, [1:(nc1-M_.exo_nbr)],j);