function  [ys,params1,check] = dsge11_steadystate(ys_init, exo_ss,M_,options_);

% function [ys,params,check] = NK_baseline_steadystate(ys,exo,M_,options_)
% computes the steady state for the NK_baseline.mod and uses a numerical
% solver to do so
% Inputs:
%   - ys        [vector] vector of initial values for the steady state of
%                   the endogenous variables
%   - exo       [vector] vector of values for the exogenous variables
%   - M_        [structure] Dynare model structure
%   - options   [structure] Dynare options structure
%
% Output:
%   - ys        [vector] vector of steady state values for the the endogenous variables
%   - params    [vector] vector of parameter values
%   - check     [scalar] set to 0 if steady state computation worked and to
%                    1 of not (allows to impose restrictions on parameters)

% Copyright © 2013-2024 Dynare Team
%
% This file is part of Dynare.
%
% Dynare is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% Dynare is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with Dynare.  If not, see <https://www.gnu.org/licenses/>.
alpha=0;
beta=0;
% read out parameters to access them with their name
NumberOfParameters = M_.param_nbr;
for ii = 1:NumberOfParameters
    paramname = M_.param_names{ii};
    eval([ paramname ' = M_.params(' int2str(ii) ');']);
end
% initialize indicator
check = 0;

%% Calculate steady state equations
N= 0.33;
q = 1;
rd = (1/betah) - 1;
rl = (1/betab) - 1;
K = ((1/betae + delta - 1)/alpha)^(1/(alpha-1));
I = delta * K;
Y = Abar * K^alpha * N^(1-alpha);
gbar = 0.2 * Y;  
t = gbar;
W = (1-alpha) * (K/N)^alpha; 
R = 1/betae;
L = q * K * theta;
d = eta * L;
ce = Y - t - I - wn - rl * L; 
cb = rl * L - rd * d - t; 
ch = wn + rd * d - t; 
me = (xi_e * ce) / (1 - betae);
mb = (xi_b * cb) / (1 - betab);
mh = (xi_h * ch) / (1 - betah);
m = me + mh + mb;
tb = 0.1 * t;
th = 0.5 * t;
te = 0.4 * t;
Abar =1;
%% end own model equations

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
