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;
theta=0.75;
mu_p=1.2;
mu_w=1.5;
gamma_g=0.2;
gamma_c=1-gamma_g-(delta*alpha)/(((1/beta - 1)+delta)*mu_p);


%% Steady State
pi = 0;
%%pi_star = 1;
p = 1;
r = 1/beta - 1;
rk = 1/beta - (1-delta);
q = 1;
%%n = (((1-alpha)/mu_p)*gamma_c^(-1))^(1/(1+psi));
n = (((1-alpha)/(mu_p*(1-(1-mu_w))))*gamma_c^(-1))^(1/(1+psi));
mc = (epsil-1)/epsil;
k = ((1/alpha)*mu_p*rk)^(1/(alpha-1))*n;
i = delta*k;
y = (k^(alpha))*(n^(1-alpha));
%%w = (1-alpha)*y/(mu_p*n);
g = 0.2*y;
c = (1-gamma_g-(delta*alpha/((r+delta)*mu_p)))*y;
b = 0;
t = g;
%%d_o = (y-1/(1-alpha)*w*n)/(1-lambda);
d_o = ((y-gamma_c*mu_w*n^(1+psi)*y-(alpha/mu_p)*y))/(1-lambda);
i_o = i/(1-lambda);
k_o = k/(1-lambda);
n_o = n;
n_r = n;
c_o = c;
c_r = c;
b_o = b/(1-lambda);
w = (mu_w*c*n^psi);

%%kappal = w/(c*n^psi);

lamdf = 1/c;
t_r = w*n-c;
t_o = (t - lambda*t_r)/(1-lambda);

%g1=c^(-1)*y*mc/(1-theta*beta);
%g2=pi_star*c^(-1)*y/(1-theta*beta);
ghat = 0;
b_norm = 0;
c_norm = c/y;
i_norm = i/y;
n_norm = n/y;
t_norm = t/y;
y_norm = y/y;





%% Non policy parameters
%%a_j = 1/0.23;
kappaq=4.5;

%% Policy parameters
%%phipi=1.5;
phipi=0.5;
rho=0.5;
rhog=0.9;
phig=0;
%%phib=0.33;
phiy=0.5;
phib=0;
%%phib=0.4;


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

%% Save parameters
save tesipar beta delta omega psi lambda epsil alpha lambdap theta phib phig phipi phiy rhog rho gamma_g gamma_c mu_p mu_w kappaq






