%Steady State File (two sector model with labor
%market frictions, exogenous separations and on-the-job search)

function[ys,check]=twosecexonew_steadystate(ys,exe)
global M_
 
% Here, the values of the deep parameters are loaded in a loop.
%
NumberOfParameters = M_.param_nbr;                            % Number of deep parameters.
for i = 1:NumberOfParameters                                  % Loop...
  paramname = deblank(M_.param_names(i,:));                   %    Get the name of parameter i. 
  eval([ paramname ' = M_.params(' int2str(i) ');']);         %    Get the value of parameter i.
end                                                           % End of the loop.  
check = 0;




% Enter the steady state model equations here


x0=[0.0846659 0.524426 0.0609509 0.0524426 0.514111 0.860408 0.468662 0.59752 0.445913 0.0353341 0.355574 0.0844513 0.027534 0.779247 0.326034 0.456579 2.39008 0.438494 0.0173404 0.0487675 0.0101306 0.721345 0.88 0.145402 0.104885 0.442215 0.4];
options=optimset('Display','iter','MaxIter',2500,'MaxFunEvals',10000);
[x,fval,exitflag]=fsolve(@ss_fun,x0,options);


% HIgh wage sector

ug_ss=x(1);
ng_ss=x(2);
vg_ss=x(3);
mg_ss=x(4);
pg_ss=x(5);
qg_ss=x(6);
epsilong_ss=x(7);
thetag_ss=x(8);
wg_ss=x(9);

% Low wage sector

ub_ss=x(10);
nb_ss=x(11);
vb_ss=x(12);
mb_ss=x(13);
pb_ss=x(14);
qb_ss=x(15);
epsilonb_ss=x(16);
thetab_ss=x(17);
wb_ss=x(18);
e_ss=x(19);
s_ss=x(20);
qr_ss=x(21);

% Aggregates, unemployment benefits, other variables

q_ss=x(22);
n_ss=x(23);
v_ss=x(24);
m_ss=x(25);
c_ss=x(26);
b=x(27);

Yb_ss=Ass*nb_ss-cb*vb_ss;
Yg_ss=Ass*ng_ss-cg*vg_ss;
yg_ss=Ass*ng_ss;
yb_ss=Ass*nb_ss;
y_ss=yb_ss^alphha*yg_ss^(1-alphha);


% End of the steady state model block

for iter = 1:length(M_.params) %update parameters set in the file
eval([ 'M_.params(' num2str(iter) ') = ' M_.param_names(iter,:) ';' ])
end




% The steady state values of the endogenous variables of the model are
% defined here.

NumberOfEndogenousVariables = M_.endo_nbr;                    % Number of endogenous variables.
ys = zeros(NumberOfEndogenousVariables,1);                    % Initialization of ys (steady state).
for i = 1:NumberOfEndogenousVariables                         % Loop...
  %varname = deblank(M_.endo_names(i,:));                      %    Get the name of endogenous variable i.                     
  %eval(['ys(' int2str(i) ') = ' varname ';']);                %    Get the steady state vZNue of this variable.
    ys(i)=0;
end                                                           % End of the loop.