News (expectation) shock routine

Hi all,

Heedon Kang of the Bank of Korea revised our code for news shock (a shock examined in Beaudry and Portier, 2004) routine with Dynare. We hope there will be no bugs any more!

You can download the codes and manual at:


Dear Ippei, I am trying to make a new shock to the Lubik and Schorfheide (2003), but I have an error, this this error also happen in some models that I have, Could you tell me what to do?

??? Error using ==> expectation_shock at 131
A12 must be zero matrix

Thanks very much in advance


I find it is useless for Dynare version 4.1, anybody offers some ideas?

I use Dynare 4.2.2 to run the model in Fujiwara and Kang (2008) “Expectation Shock Simulation with Dynare: Mark II”:
(1) type “dynare pigoucycle”
(2) Then run the file “expectation_shock.m”.

The following errors appear:

??? Error using ==> feval
Undefined function or method ‘_fff’ for input arguments of type ‘double’.

Error in ==> resol_es at 39
if max(abs(feval(fh,ys))) > dynatol_ & options_.olr == 0

Error in ==> stoch_simul_es at 47
[dr_,info,jacobiamat] = resol_es(ys_,0);

Error in ==> expectation_shock at 25
[info,jacobimat]=stoch_simul_es(var_list_); % To obtain Jacobian matrix

I appreciate very much if the above problems can be solved. Please find the attached key files for reference. The full toolkit files can be downloaded from:

Thanks for your attention.
resol_es.m (1.82 KB)
stoch_simul_es.m (12.2 KB)
expectation_shock.m (22.7 KB)
pigoucycle.mod (2.19 KB)