function [ys,check] = ncgm1_steadystate(ys,exe);

global M_ options_ oo_


%%-------------------------------------------------------------------------
%

NumberOfParameters = M_.param_nbr;

for ii = 1:NumberOfParameters
    paramname = deblank(M_.param_names(ii,:));
    eval([ paramname ' = M_.params(' int2str(ii) ');']);
end

% Code runs when using lines below:
% alphaa    = M_.params(1);
% betta     = M_.params(2);
% delta     = M_.params(3);
% rho       = M_.params(4);
% sigma_c   = M_.params(5);
% sigma_eps = M_.params(6);
%%-------------------------------------------------------------------------


k = (alphaa/(1/betta - (1 - delta)))^(1/(1 - alphaa));
y = k^(alphaa);
Inv = delta*k;
c = y - Inv;
a = 1;
w = (1-alphaa)*k^(alphaa);
R = alphaa*k^(alphaa-1);
r = (1/betta) - 1;

check = 0;

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

NumberOfEndogenousVariables = M_.endo_nbr;
ys = zeros(NumberOfEndogenousVariables,1);

for ii = 1:NumberOfEndogenousVariables
    varname = deblank(M_.endo_names(ii,:));
    eval(['ys(' int2str(ii) ') = ' varname ';']);
end
ysold = ys;

% Code runs when using lines below:
% ys = [log(y);
% log(Inv);
% log(k);
% log(a);
% log(c);
% log(w);
% log(R);
% r];
%%-------------------------------------------------------------------------