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

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



%%%% Model equations to be entered here
    alphaa  = 0.32;
		rhoo    = 0.42;
		phii    = 0.028;
		r_w    = 0.04;		
        gammaa  = 2;
		omegaa  = 1.455;
		psii    = 0.000742;
		deltaa  = 0.1;
		sigmae = 0.0129;
		betaa   = 1/(1+r_w);
		h_ss   = ((1-alphaa)*(alphaa/(r_w+deltaa))^(alphaa/(1-alphaa)))^(1/(omegaa-1)); 
		k_ss   = h_ss/(((r_w+deltaa)/alphaa)^(1/(1-alphaa)));
		i_ss   = deltaa*k_ss;                                                     
		y_ss   = (k_ss^alphaa)*(h_ss^(1-alphaa));                                   
		d_bar  = 0.7442;
		d_ss   = d_bar;                                                        
		c_ss   = y_ss-i_ss-r_w*d_ss;
		tb_ss  = y_ss-c_ss-i_ss;
        
% vars
    d     = d_ss;
    c     = log(c_ss);
    h     = log(h_ss); 
    y     = log(y_ss);
    i     = log(i_ss);
    k     = log(k_ss);
    a=0;
    lambdaa= log((exp(c)-((exp(h)^omegaa)/omegaa))^(-gammaa));
    tb    =  1-((exp(c)+exp(i))/exp(y));
    ca=0;
    riskpremium=0;
    r     = log((1-betaa)/betaa);
    
    
% from dynare's solved ss     
d  =         		 0.744200000000000; 
c  =         		0.110602456369385; % 0.110602456369385
h  =         		 0.007390615600776; % 0.007390615600776
y  =         		 0.396415826511114; % 0.396415826511114
i  =         		-1.079490693298464; % -1.079490693298464
k  =         		 1.223094399695582; % 1.223094399695582
a  =         		 0;
lambdaa =    		 1.724386196437486; % 1.724386196437486
tb =         		 0.020025734361830; % 0.020025734361830
ca =         		 0;
riskpremium =		 0;
r  =         		 -3.218875824868198; % -3.218875824868198


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);
for i = 1:NumberOfEndogenousVariables
  varname = deblank(M_.endo_names(i,:));
  eval(['ys(' int2str(i) ') = ' varname ';']);
end
else
ys=zeros(length(lgy_),1);
for i = 1:length(lgy_)
    ys(i) = eval(lgy_(i,:));
end
check = 0;
end