%close all;

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

// parameters

parameters beta kappa sigma mu nu;

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

// model
model(linear);
pie = 2 + (beta)*(pie(+1)-2) + (1-beta)*(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;

// values of lagged variables

histval;
pie(0)=2;
y(0)=0;
ir(0)=5;
end;

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;

