I absolutely need to find out the problem of my simulation for this weekend. Could you help me to fix it ?
Here is my model :
- Code: Select all
// Declaration of variables
var c d k f r p V g;
varexo R;
parameters a wy wo V0;
// Calibration of parameters
a = 0.2;
wy = 8;
wo = 4;
V0=2;
// The model
model;
f = (1+R)*k(-1)^a;
k = (d(+1) / (c * (1 + R)*a))^(1 / (a-1));
k = wy - c - g*V;
d(+1) = wo + k * (1+r) + g* (V(+1) + p(+1));
V + p = (1+r) * V(-1);
d(+1) / c = (1+r);
p = f - (1+r) * k(-1);
c + d + k = wy + wo + f(-1);
end;
initval;
c = 4;
d = wo + V0;
k = 2;
g=1;
end;
check;
steady;
and the result :
- Code: Select all
>> dynare simulation_test1
Configuring Dynare ...
[mex] Generalized QZ.
[mex] Sylvester equation solution.
[mex] Kronecker products.
[mex] Sparse kronecker products.
Starting Dynare ...
Starting preprocessing of the model file ...
8 equation(s) found
Processing derivation ...
Processing Order 1... done
Processing Order 2... done
Processing outputs ...
Preprocessing completed.
Starting Matlab computing ...
??? Error using ==> mrdivide
Matrix dimensions must agree.
Error in ==> lnsrch1 at 64
alamin = tolx/test ;
Error in ==> solve1 at 126
[x,f,fvec,check]=lnsrch1(xold,fold,g,p,stpmax,func,j1,j2,varargin{:});
Error in ==> dynare_solve at 110
[x,info]=solve1(func,x,j1(r(i):r(i+1)-1),j2(r(i):r(i+1)-1),jacobian_flag,varargin{:});
Error in ==> resol at 75
[dr.ys,check1] =
dynare_solve(fh,dr.ys,options_.jacobian_flag,...
Error in ==> check at 48
[dr, info] = resol(oo_.steady_state,1);
Error in ==> simulation_test1 at 103
check;
Error in ==> dynare at 102
evalin('base',fname) ;
Thanks a lot,
Adrien