function [ys,check] = model11_steadystate(ys,exe);
global M_ lgy_ options_

if isfield(M_,'param_nbr') == 1
NumberOfParameters = M_.param_nbr;
for i = 1:NumberOfParameters
  paramname = deblank(M_.param_names(i,:));
  eval([ paramname ' = M_.params(' int2str(i) ');']);
end
check = 0;
end


ETA       =0.0125;

EPSILON   = 0.39;

OMEGA     = 0.105;

KAPPA_bar = .75; 

N         = 0.5;

SIGMA     = 2; 

RHO_a     = 0.9;

RHO_kappa = 0.9; 

SIGMANUi  = 0.006;

SIGMANUk  = 0.00819;

ri = 0.01;     




rb_s = (1+ri);

BETAs_s = 1/rb_s;

BETAi_s = 1/(1+ri+0.002);

a_s  = 1; 

kappa_t_s = KAPPA_bar;

ki_s = 1/(1.1*N);

ks_s = 0.10*(N*ki_s)/(1-N);

r_s  = (1-EPSILON)*N^(1-EPSILON)*(N*ki_s)^(-EPSILON);

w_s = (N*ki_s)^(1-EPSILON) - r_s*(N*ki_s);

q_s = ((BETAi_s*r_s)/(1-(KAPPA_bar*0.002*BETAi_s)-BETAi_s));

Z1_s = ((1-BETAs_s)*q_s)/(BETAs_s*OMEGA*(ks_s)^(OMEGA-1));

bi_s = KAPPA_bar*q_s*ki_s;

bs_s = -(N*bi_s)/(1-N);

ci_s = w_s + (r_s * ki_s) +(1-rb_s)*bi_s;

cs_s = w_s + Z1_s *(ks_s)^OMEGA + (1-rb_s)*bs_s;

%cs_s = (1/(1-N))*(a_s*(N*ki_s)^(1-EPSILON)+(1-N)*Z1_s*(ks_s)^OMEGA- N*ci_s);


mu_s = 0.002*(BETAi_s*(ci_s)^(-SIGMA));

ZETAs_s = (BETAs_s)/((1+(1-N)*cs_s)^(-ETA));

ZETAi_s = (BETAi_s)/((1+N*ci_s)^(-ETA));


ZETAs = ZETAs_s;
ZETAi = ZETAi_s;
BETAi = BETAi_s;
BETAs = BETAs_s;
Z1 = Z1_s;
ci = log(ci_s);
cs = log(cs_s);
ki = log(ki_s);
ks = log(ks_s);
r = log(r_s);
rb = log(rb_s);
w = log(w_s);
a = log(a_s);
kappa_t = log(kappa_t_s);
%log(KAPPA_bar);
bi = bi_s;
bs = bs_s;
%the value of bs_s is negative, therefore no log()
%mu = log(mu_s);
mu = mu_s;

q = log(q_s);
%ZETAs = ZETAs_s;


for iter = 1:length(M_.params)
  eval([ 'M_.params(' num2str(iter) ') = ' M_.param_names(iter,:) ';' ])
end

if isfield(M_,'param_nbr') == 1

if isfield(M_,'orig_endo_nbr') == 1
NumberOfEndogenousVariables = M_.orig_endo_nbr;
else
NumberOfEndogenousVariables = M_.endo_nbr;
end
ys = zeros(NumberOfEndogenousVariables,1);
if options_.linear
    
else
for i = 1:NumberOfEndogenousVariables
  varname = deblank(M_.endo_names(i,:));
  eval(['ys(' int2str(i) ') = ' varname ';']);
end
end

else
ys=zeros(length(lgy_),1);
for i = 1:length(lgy_)
    ys(i) = eval(lgy_(i,:));
end
check = 0;
end