clc; clear all; close all

%% Structural parameter in the steady state
beta=0.99;
delta=0.025;
alpha=1/3;
epsil=6;
lambda=0.5;
psi=0.2;


%% Steady State
pi = 0;
p = 1;
r = 1/beta - 1;
rk = 1/beta - (1-delta);
q = 1;
n = 1/3;
mc = (epsil-1)/epsil;
k = ((rk/(alpha*mc))^(1/(alpha-1)))*n;
i = delta*k;
y = (k^(alpha))*(n^(1-alpha));
w = (1-alpha)*mc*(k/n)^alpha;
g = 0.2*y;
c = y-i-g;
b = 0;
t = g;
d_o = y*(1-(epsil-1)/epsil);
i_o = i/(1-lambda);
k_o = k/(1-lambda);
n_o = n;
n_r = n;
c_o = c;
c_r = c;
kappal = w/(c*n^psi);
lamdf = 1/c;
t_r = w*n-c;
t_o = (g - lambda*t_r)/(1-lambda);




%% Non policy parameters
theta=0.75;
kappaq=4;

%% Policy parameters
phipi=1.5;
rhog=0.9;
phig=0.1;
phib=0.33;


%% Parameters values depending on other parameters
omega=(alpha^(alpha))*(1-alpha)^(1-alpha);
lambdap=(1-beta*theta)*(1-theta)*(theta)^(-1);

%% Save parameters
save tesipar beta delta omega psi lambda epsil alpha lambdap theta phib phig phipi rhog kappaq kappal






