

function outc = Code_constant_calibration(inc)

% Function that contains the equations for the steady state

global delta gamma zss ess eta xi phi beta mu;
global  h1;


n      = inc(1); 
vs     = inc(2);
w      = inc(3);
s      = inc(4);
m      = inc(5);
qv     = inc(6);
qs     = inc(7);
y      = inc(8);
v      = inc(9);
C      = inc(10);



%Calibration targets
h1       = inc(11);
phi      = inc(12);
gamma    = inc(13);
mu       = inc(14);





%Law of motion for employment
outc(1,1)  = -n + (1 - delta)*(n) + m;

%Job creation condition
outc(2,1)  = -gamma/(qv)+zss-w+(1-delta)*beta*(gamma/(qv));

%Nash wage
%outc(3,1)  = -w+eta*zss+(1-eta)*phi+(1-eta)*(1-delta)*beta*(hprime/(1/C)-phi);

outc(3,1)  = -w+eta*zss+(1-eta)*phi-(1-eta)*(1-delta)*phi;


outc(4,1)   = -((1 - delta)*n + s) + 1;


%Matching function
outc(5,1)  = -m + mu*(v)^(1-xi)*(ess*s)^xi;

%Job-filling probability
outc(6,1)  = -qv + m/(v);

%Job-finding probability
outc(7,1)  = -qs + m/(s);

%Output
outc(8,1)  = -y + zss*n;



%Vacancies
outc(9,1) = -vs + v/(s);



%Resource constraint
outc(10,1) = -(C + gamma*v) + y;




outc(11,1)=-h1+(qs/ess)*(w-phi + (1-delta)*beta*((-(1-qs)*phi)/(qs)));

%Calibrated parameters




%Value of unemployment (chi)
outc(12,1) = -phi + 0.50*w;

%Cost of posting vacancies (gammag): match (gamma*v) = 0.05*y
outc(13,1) = -qv + 0.70;

%Matching efficiency (M): match ftheta = 0.80
outc(14,1) = -qs + 0.60;

end