diff --git a/matlab/UnivariateSpectralDensity.m b/matlab/UnivariateSpectralDensity.m
index 46fdf895da906f29579593700cb6060b21c3c42c..2757da15b010e8518d0473217c316c4da6168435 100644
--- a/matlab/UnivariateSpectralDensity.m
+++ b/matlab/UnivariateSpectralDensity.m
@@ -33,10 +33,9 @@ if options_.order > 1
     disp('Set order = 1.')
     return
 end
-
-pltinfo  = 1;%options_.SpectralDensity.Th.plot;
-cutoff   = 150;%options_.SpectralDensity.Th.cutoff;
-sdl      = 0.01;%options_.SepctralDensity.Th.sdl;
+pltinfo  = options_.SpectralDensity.plot; 
+cutoff   = options_.SpectralDensity.cutoff; 
+sdl      = options_.SepctralDensity.sdl; 
 omega    = (0:sdl:pi)';
 GridSize = length(omega);
 exo_names_orig_ord  = M_.exo_names_orig_ord;
diff --git a/matlab/global_initialization.m b/matlab/global_initialization.m
index 24e0bdcaec2edc120a3467629b6ac65ddced964b..01203e97b69cae3444db0de30bf59b3cfc21f08b 100644
--- a/matlab/global_initialization.m
+++ b/matlab/global_initialization.m
@@ -125,7 +125,10 @@ options_.nomoments = 0;
 options_.nocorr = 0;
 options_.periods = 0;
 options_.noprint = 0;
-options_.SpectralDensity = 0;
+options_.SpectralDensity.trigger = 0;
+options_.SpectralDensity.plot  = 1; 
+options_.SpectralDensity.cutoff  = 150; 
+options_.SpectralDensity.sdl = 0.01; 
 
 % Extended path options
 %
diff --git a/matlab/stoch_simul.m b/matlab/stoch_simul.m
index 62edee0618afa6e540e4ae6ed689733769cea389..2295c6d0f13c661bf113b67a44cf66024c47aaf6 100644
--- a/matlab/stoch_simul.m
+++ b/matlab/stoch_simul.m
@@ -323,7 +323,7 @@ if options_.irf
     end
 end
 
-if options_.SpectralDensity == 1
+if options_.SpectralDensity.trigger == 1
     [omega,f] = UnivariateSpectralDensity(oo_.dr,var_list);
 end
 
diff --git a/matlab/stoch_simul_sparse.m b/matlab/stoch_simul_sparse.m
index fe5eab4bdb6684f164c741df6e95b17b939d4d7f..1eef7ba2b32ef6c0b3c83af54f289ca871edea25 100644
--- a/matlab/stoch_simul_sparse.m
+++ b/matlab/stoch_simul_sparse.m
@@ -286,7 +286,7 @@ if options_.irf
     end
 end
 
-if options_.SpectralDensity == 1
+if options_.SpectralDensity.trigger == 1
     [omega,f] = UnivariateSpectralDensity(oo_.dr,var_list);
 end