* Áydano Ribeiro Leite/URCA
* Modelo Novo Keynesiano/ Dominância fiscal 
* Código no Dynare
* Áydano Ribeiro Leite/URCA
//Variáveis Endógenas
var Y,c,I,R,K,W,N,A,P,Pj,Pa,m,b,RT,CM,SX,SM,D;
//Variáveis Exógenas
varexo e e_G e_tauc e_tau e_Zr;
//Parâmetros
parameters alpha, beta,delta,gamma,gammar,gammapi,gammay,gammab,phik,eta,psim,psiy,xi,rhoA,rhoG,rhotauc,rhotau,rhoZr,sigmaa,sigmag,sigmatauc,sigmatau,sigmar;
//Calibração
alpha = 0.448;
beta = 0.989;
delta = 0.015;
gamma = 2.4;
gammar= 0.79;
gammapi= 2.43;
gammay = 0.16;
gammab= 0.50;
phik = 6.599;
eta = 1.6553;
psim = 0.25;
psiy = 8;
xi = 0.74;
rhoA = 0.8;
rhog = 0.8;
rhotauc = 0.8;
rhotau = 0.8;
rhoZr = 0.5;
sigmaa = 0.01;
sigmag = 0.01;
sigmatauc = 0.02;
sigmatau = 0.02;
sigmar = 0.01;
model;
//Equações do Modelo não-linear
(1+tauc)*c+I+(phik/2)*(I/K(-1)-delta)^2*K(-1)+m+b =(1-tau)*(W*N+R*K(-1))+tau*delta*K+m(-1)*PI+b(-1)*PI*RT(-1);
//Condições de Primaira Ordem p/Problema da Família
lambda=1/(1+tauc)*c;
lambda=gamma*(m)^(-1/psim)+beta*(lambda(+1)/PI(+1));
lambda=eta/(1-tau)*(1-N)*W;
lambda*(1+phik*(I/K(-1)-delta)=beta*(lambda(+1)*(1+(1-tau(+1))*(R(+1)-delta)+phik*(I(+1)/K-delta)+(phik/2)*(I(+1)/K-delta)^2);
m=(gamma*(1+tauc)*c*(RT/RT-1))^(psim);
c=(1-tau)*(W*N+R*K(-1))+tau*delta*K+m(-1)/PI-m+b(-1)/PI*RT(-1)-b-I-(phik/2)*(I/K(-1)-delta)^2/(1+tauc) 
//Lei do Movimento do Capital
K=(1-delta)*k(-1)+I;
//Firmas
Y=A*(K^alpha)*(N^(1-alpha));
log(A)=rho*log(A(-1))+e;
W=(1-alpha)*A*(K^alpha)*(N^alpha);
R=alpha*A*(K^(alpha-1))*(N^(1-alpha));
N=(1-alpha)*CM*(Y/W);
K=alpha*CM*(Y/R);
CM=(W^(1-alpha))/A*(1-alpha)^(1-alpha))*(R^alpha)/alpha^alpha));
//Fixação do Preço ótimo por Parte das Firmas e Regra de Rigidez
Pj=(psiy/psiy-1)*(1/1-beta*xi)*CM;
Pa=((1-xi)*(Pj^(1-psiy))+xi*(Pa(-1)^(1-psiy))^(1/1-psiy));
//16-Taxa de inflação bruta                                               
PI = P - P(-1);
//Governo
//Restrição Orçamentária do Governo
G+b(-1)*(RT(-1)-1)=tauc*c+tau*(W*N+R*K(-1)-delta*K)+(b-b(-1)+(m-m(-1);
b(-1)=P*tauc*c+tau(W*N+R*K(-1)-delta*K)+(b-b(-1))+(m-m(-1))-G/(RT(-1)-1);
//Processos Estocásticos para as Variáveis Fiscais
log(G/Gss)=rhog*log(G(-1)/Gss)+e_G;
log(tauc/taucss)=rhotauc*log(tauc(-1)/taucss)+e_tauc;
log(tau/tauss)=rhotau*log(tau(-1)/tauss)+e_tau;
//Superávit Primário e Receita de Senhoriagem
SX=tauc*c+tau*(W*N+R*K(-1)-delta*K)-G;
SM=(m-m(-1));
//Regra de Política Fiscal de Curto Prazo
//Evolução da Dívida para a Regra de Política de Curto Prazo
D(+1)=b*(RT-1)+G-tauc*c+tau*(W*N+R*K(-1)-delta*K);
//Regra de Taylor
RT=(RT(-1))^(gammar)*((1/beta+(PI(+4))^(gammapi/4)*(Y/Yss)^(gammay)*(b/bss)^(-gammab))^(1-gammar)*Z_r;
//Condições de Equilíbrio
Y=c+I+G+(phik/2)*(I/K(-1)-delta)^2*K(-1);
m=M;
b=B;
end;
initval;
//Equações de Steady State
Pss=1;
Ass=1;
PIss=1;
Rss=(1/beta-1)*(1/1-tauss)+delta;
RTss=(1/beta);
mss= (gamma*(1+taucss)*css*(1/1-beta))^psym;
Pjss=(psiy/psiy-1)*(1/1-beta*xi)*CM;
Pass=((1-xi)*Pjss+xi*Pass^(1-psiy))^(1/1-psiy);
CMss=(psiy/psiy-1)*(1-beta*xi)*Pjss;
Wss= (1-alpha)*CMss^(1/1-alpha)*(alpha/Rss)^alpha/1-alpha;
Nss= (1-alpha)*CMss*Yss/Wss;
Kss= alpha*CMss*Yss/Rss;
Iss=delta*Kss;
css=(1-beta)*(1-taucss)/(1+beta*(-1+taucss)*(1-RTss)*((1-tauss)*(Wss*Nss+Rss*Kss)/(1+taucss)-(1-tauss)*(delta*Kss/1+taucss)-(1-RTss)/1+taucss)*(beta*taucss*(Wss*Nss+Kss*(Rss-delta)-beta*G)/(1-beta);
bss=beta*taucss*css+tauss*(Wss*Nss+Kss*(Rss-delta)-G/1-beta;
G=Gss;
tauc=tacss;
tau=tauss;
SXss=taucss*css+tau*(Wss*Nss+Kss*(Rss-delta)-G;
SMss=0;
Dss=bss*(RTss-1)+Gss-taucss*css+tauss*(Wss*Nss+Kss*(Rss-delta);
Yss=css+Iss+Gss;
end;

steady;
check;

shocks;
var e= 0.01;
var e_G=0.01; 
var e_tau_c=0.02; 
var e_tau=0.02;
var e_Zr=0.01;
end;

stoch_simul Y,c,I,R,K,W,N,A,P,Pj,Pa,m,b,RT,CM,SX,SM,D;