%close all;

// variables
var pie y ir;
varexo ey ep ei;

// parameters

parameters theta kappa sigma mu nu;

theta = 0.6;
kappa = 1;
sigma = 1/5;
mu = 1.5;
nu = 1;

// model
model(linear);
pie = 2 + (theta)*(pie(+1)-2) + (1-theta)*(pie(-1)-2) + kappa*y + ep;
y = y(+1) - 1/sigma*(ir-pie(+1)-3) + ey;
ir = 5 + mu*(pie-2) + nu*y + ei;
end;

// initial values
initval;
ir = 3;
pie = 3;
y=0;
ep=0;
ey=0;
ei=0;
end;

steady;

shocks;
var ey; periods 1; values -5;
var ep; periods 1; values 5; 
//var ei; periods 1; values 5; 
end;


simul(periods=20);
rplot pie y ir;

