function [ys,check] = CC1_steadystate(ys,exe)
global M_

%% 
%%
NumberOfParameters = M_.param_nbr;                           
for i = 1:NumberOfParameters                                 
  paramname = deblank(M_.param_names(i,:));                
  eval([ paramname ' = M_.params(' int2str(i) ');']);         
end                                                             
check = 0;
%%
%%
w=-0.1461;
N= -0.4771; 
C=0.6021; 
q=0; 
D=6;
pi_d=0; 
e= 0.01; 
R_n=1.0101;
pi_c=0; 
R_star_n=1.0101; 
psi=-2;
b=1; 
b_star= 0; 
b_star_tilde= 0;
T= 0;
N_tilde=-0.4771;
C_tilde=0.699; 
D_tilde=0.8451; 
b_tilde=1.0414; 
N_c=0.9031;
N_d=0.9542; 
p_h= 0;
p_f= 0; 
Ytot=-0.4771; 
epsilon_m=0; 
durshare=0.699;
shadow_dur=1.20236;
shadow_dur_tilde=0.674256;
S= 0;
pi_c_h= 0;

S = 1;            
q = 1;            
pi_d = 1;            
%e = 1;             
pi_c = 1;           
p_h = 1;            
p_f = 1;         
pi_c_h = 1;           

%w=log((epsilon_c-1)/epsilon_c);
N= 1/3; 
C=4; 
q=0; 
D=6;
pi_d=0; 
e= 1;
%R_n=log(1/beta);
pi_c=0; 
%R_star_n=log(1/zeta);
%psi=zeta-beta;
b=10; 
b_star= 0;
b_star_tilde= 0;
T= 0;
N_tilde=1/3;
C_tilde=5; 
D_tilde=7; 
b_tilde=11; 
N_c=8;
N_d=9; 
p_h= 0;
p_f= 0;
Ytot=1/3; 
epsilon_m=0; 
durshare=log(0.2);
Dtot=log(omega*exp(D)+(1-omega)*exp(D_tilde));
Ctot=log(omega*exp(C)+(1-omega)*exp(C_tilde));
CDtot=log(omega*exp(D)+(1-omega)*exp(D_tilde));

%%
NumberOfEndogenousVariables = M_.endo_nbr;                 
ys = zeros(NumberOfEndogenousVariables,1);                
for i = 1:NumberOfEndogenousVariables                      
  varname = deblank(M_.endo_names(i,:));                                    
  eval(['ys(' int2str(i) ') = ' varname ';']);             
end                                                        
%%