I have simulated a stochastic model using the second order approximation of the model and now want to see how the model reacts to a determinstic shock. I have used the forecast function for that following the manuals instructions. I receive the response below from Matlab (the steady state values are fine). I have attached the file in case you have the time to take a look. The derministic shock is introduced in line 3 (called tau) and the shocks are calibrated in line 91 to 97. Your help would be much appreciated!
??? Attempted to access x(2,3); index out of bounds because size(x)=[3,2].
Error in ==> stoch1_dynamic at 110
rhs =x(it_, 1)+x(it_, 3);
Error in ==> dr1 at 177
[junk,jacobia_] = feval([M_.fname ‘_dynamic’],z,tempex);
Error in ==> resol at 64
[dr,info] = dr1(dr,check_flag);
Error in ==> check at 23
[dr, info] = resol(oo_.steady_state,1);
Error in ==> stoch1 at 240
Error in ==> dynare at 49
stoch1.mod (2.05 KB)