Expectation shock codes

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: ideas.repec.org/c/dge/qmrbcd/163.html

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


It looks like this code has been written for Dynare version 3.

This is an old version that we no longer support.


I tried Version 3, the following errors appear:

??? Error using ==> dynare at 19
syntax error at line 29

Can someone hlep me to fix it?
Thanks a lot!