% computes the steady state of rbc_baseline.mod (Hansen model). 


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

ctheta = M_.params(1);
cdelta = M_.params(2);
crho = M_.params(3);
cbeta = M_.params(4);
cgamma = M_.params(5);
a = M_.params(6);
nstar = M_.params(7);
%%
%% END OF THE FIRST MODEL INDEPENDENT BLOCK.


%% THIS BLOCK IS MODEL SPECIFIC.
%%
%% Here the user has to define the steady state.
%%
%z   = 0;
nstar   = 1/3;
r   = (1-cbeta)/ cbeta;
y_k = (r+cdelta)/ctheta;
y   = (cgamma*y_k^(-ctheta)*nstar^(1-ctheta))^(1/(1-ctheta));
k   = (y_k^(-1))*y;
c   = y-cdelta*k;
a   = 1/c*(1-ctheta)*(y/nstar);
n   = nstar;


check = 0;
%%
%% 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.
%%
% NumberOfEndogenousVariables = M_.endo_nbr;                    % Number of endogenous variables.
% ys = zeros(NumberOfEndogenousVariables,1);                    % Initialization of ys (steady state).
% for i = 1:NumberOfEndogenousVariables                         % Loop...
%   varname = deblank(M_.endo_names(i,:));                      %    Get the name of endogenous variable i.                     
%   eval(['ys(' int2str(i) ') = ' varname ';']);                %    Get the steady state value of this variable.
% end                                                           % End of the loop.
% %%
% %% END OF THE SECOND MODEL INDEPENDENT BLOCK.