diff --git a/dynare++/dynare_simul/dynare_simul.m b/dynare++/dynare_simul/dynare_simul.m
index 095d417b7a3489fa2c4e5390209c4d00a9df889c..408010b72d924f0ad4ddbdc01f861e6f91bd3af5 100644
--- a/dynare++/dynare_simul/dynare_simul.m
+++ b/dynare++/dynare_simul/dynare_simul.m
@@ -42,6 +42,11 @@
 % All examples suppose that the prefix is 'dyn' and that your_model.mat
 % has been loaded into Matlab.
 %
+% You could e.g. use
+%       load('example1.mat');
+%       shocks=randn(size(dyn_shocks,1),1000)'*chol(dyn_vcov_exo);
+%       dynare_simul('example1',shocks')
+%
 % 1. response to permanent negative shock to the third exo var EPS3 for
 %    100 periods
 %
@@ -168,9 +173,5 @@ end
 seed = ceil(10000*rand(1,1));
 
 % call dynare_simul_
-[err,r]=dynare_simul_(order-1,nstat,npred,nboth,nforw,...
-    nexog,ystart,shocks,vcov_exo,seed,ss,dr);
-
-if err
-    error('Simulation failed')
-end
\ No newline at end of file
+r=dynare_simul_(order-1,nstat,npred,nboth,nforw,...
+    nexog,ystart,shocks,vcov_exo,seed,ss,dr);
\ No newline at end of file