diff --git a/matlab/dynare_sensitivity.m b/matlab/dynare_sensitivity.m
index efc847a5ecc91e258a619523f927af9318d91c41..ed9af2e1d2a97d1d46fb30571f6c9a7ba1f9da5c 100644
--- a/matlab/dynare_sensitivity.m
+++ b/matlab/dynare_sensitivity.m
@@ -214,6 +214,7 @@ for ii=1:length(options_gsa.var_rmse)
     temp_name = M_.endo_names_tex{strmatch(options_gsa.var_rmse{ii}, M_.endo_names, 'exact')};
     options_gsa.var_rmse_tex = vertcat(options_gsa.var_rmse_tex, temp_name);
 end
+options_gsa.var_rmse_tex = cellfun(@(x) horzcat('$', x, '$'), options_gsa.var_rmse_tex, 'UniformOutput', false);
 options_gsa = set_default_option(options_gsa,'pfilt_rmse', 0.1);
 options_gsa = set_default_option(options_gsa,'istart_rmse', options_.presample+1);
 options_gsa = set_default_option(options_gsa,'alpha_rmse', 0.001);
diff --git a/tests/gsa/ls2003.mod b/tests/gsa/ls2003.mod
index 47b806e34f41505a279a52c2e57aa36923061413..25e02deffe02f9b97a50305bc10b452d59b3682c 100644
--- a/tests/gsa/ls2003.mod
+++ b/tests/gsa/ls2003.mod
@@ -1,4 +1,4 @@
-var y y_s R pie dq pie_s de A y_obs pie_obs R_obs ${R^{obs}}$;
+var y y_s R pie $\pi$ dq pie_s de A y_obs pie_obs R_obs ${R^{obs}}$;
 varexo e_R ${e^R}$ e_q e_ys e_pies e_A;
 
 parameters psi1 psi2 psi3 rho_R ${\rho_R}$  tau alpha rr k rho_q rho_A rho_ys rho_pies;