diff --git a/doc/manual/source/the-model-file.rst b/doc/manual/source/the-model-file.rst
index 9e25601ba4387f28b4614b66964c325fdf1dafab..68acc1c0dbdbf589bd11afc2981caaa3a895a93d 100644
--- a/doc/manual/source/the-model-file.rst
+++ b/doc/manual/source/the-model-file.rst
@@ -3407,8 +3407,8 @@ Computing the stochastic solution
 
     *Output*
 
-    This command sets ``oo_.dr``, ``oo_.mean``, ``oo_.var`` and
-    ``oo_.autocorr``, which are described below.
+    This command sets ``oo_.dr``, ``oo_.mean``, ``oo_.var``, ``oo_.var_list``,
+    and ``oo_.autocorr``, which are described below.
 
     If the ``periods`` option is present, sets ``oo_.skewness``,
     ``oo_.kurtosis``, and ``oo_.endo_simul`` (see
@@ -3465,6 +3465,10 @@ Computing the stochastic solution
     an approximation thereof for ``order=2``), and simulated variance
     otherwise. The variables are arranged in declaration order.
 
+.. matvar:: oo_.var_list
+
+    The list of variables for which results are displayed.
+
 .. matvar:: oo_.skewness
 
     |br| After a run of ``stoch_simul`` contains the skewness
diff --git a/matlab/stoch_simul.m b/matlab/stoch_simul.m
index 9efb10530b3121f737973a843b73dd4c3ccb830a..190796b28ecef648d94c4768bbb04c33fb902857 100644
--- a/matlab/stoch_simul.m
+++ b/matlab/stoch_simul.m
@@ -65,6 +65,7 @@ end
 
 [i_var, nvar, index_uniques] = varlist_indices(var_list, M_.endo_names);
 var_list=var_list(index_uniques);
+oo_.var_list = var_list;
 
 iter_ = max(options_.periods,1);
 if M_.exo_nbr > 0