function [ys,params,check] = test_250624_steadystate(ys,exo,M_,options_)% for dynare 4.5.6

 
%% DO NOT CHANGE THIS PART.
%%
%% Here we load the values of the deep parameters in a loop.
%%
NumberOfParameters = M_.param_nbr;
for ii = 1:NumberOfParameters
  paramname = M_.param_names{ii};
  eval([ paramname ' = M_.params(' int2str(ii) ');']);
end
% initialize indicator
check = 0;
%
%%


%% THIS BLOCK IS MODEL SPECIFIC.
%%
%% Here the user has to define the steady state.
%%

pie = 1;
A = 1;
q = 1;
R = 1/betap-1;
Rb = (1+R)*mul-1;
pstar = 1;                                                    
pe = (epsilonp-1)/epsilonp;                                  
vp = 1;   
mb = mb_ss;
Lambdap = betap;

% solve for eight vars using eq. 2,4,5,6,8,9,10,21
fun = @(x)root210624(x, chi, etap, pe, deltak, Rb, alphae, etapen, mb, betae );
x0 = [2, 6, 1.5, 1, 3, 0.1, 0.1, 0.6];
x = fsolve(fun,x0)

Y = x(1);
k = x(2);
w = x(3);
l = x(4);
b = x(5);
psib = x(6);
ce = x(7);
cp = x(8);

I = deltak*k;
lambdap = 1/cp;
lambdae = 1/ce;
x1 = (pe*Y)/(1-phip*betap);                                
x2 = (Y)/(1-phip*betap);                                   
ye = Y;
C = ce + cp;
slope = (1+Rb)*psib;

ss_leverage = ((1+Rb)*b)/((1-deltak)*k);

%% END OF THE MODEL SPECIFIC BLOCK.

%% DO NOT CHANGE THIS PART.
%%
%% Here we define the steady state values of the endogenous variables of
%% the model.
%%
% %%
params=NaN(NumberOfParameters,1);
for iter = 1:length(M_.params) %update parameters set in the file
  eval([ 'params(' num2str(iter) ') = ' M_.param_names{iter} ';' ])
end

NumberOfEndogenousVariables = M_.orig_endo_nbr; %auxiliary variables are set automatically
for ii = 1:NumberOfEndogenousVariables
  varname = M_.endo_names{ii};
  eval(['ys(' int2str(ii) ') = ' varname ';']);
end

                                                           % End of the loop.
%%