function [ys,check] = in_class_dynare_parameters_saved_steadystate(ys,exe)

global M_

alpha = M_.params(1);
beta = M_.params(2);
delta = M_.params(3);
phi = M_.params(4);
theta = M_.params(5);
chi = M_.params(6);
rhoa = M_.params(7);
rhopsi = M_.params(8);
sa = M_.params(9);
spsi = M_.params(10);
psiss = M_.params(11);
kn_ss = M_.params(12);
mu_ss = M_.params(13);


q = 1;
%mu_ss = 0;
mu_ss = ((1-alpha)/psiss)*((1/beta)- (1-delta))-alpha;
kn_ss = (alpha/(1/beta -(1-delta)-mu_ss*psiss))^(1/(1-alpha));
n = ((1/theta)*(1/(1+mu_ss))*((1-alpha)*(kn_ss)^(alpha))/((kn_ss)^(alpha) - delta*(kn_ss)))^(1/(1+chi));
k = n*(alpha/(1/beta - (1-delta) - mu_ss*psiss))^(1/(1-alpha));
psi = psiss;
A = 0;
r = (1/beta) - 1;
y = k^alpha*n^(1-alpha);
i = delta*k;
c = y - i;
lambda = 1/c;
w = psiss*kn_ss;
mu = mu_ss;
%w = c*theta*n^chi;


check = 0;

ys = [log(y) 
    log(i) 
    log(k) 
    log(A) 
    log(c)
    log(n)];

% psi mu lambda w q r

end



%%% SS conditions when kn_ss and mu_ss are plugged in. 
% q = 1;
% n = ((1/theta)*(1/(1+(((1-alpha)/psiss)*((1/beta)- (1-delta))-alpha)))*((1-alpha)*(((alpha/(1/beta -(1-delta)-(((1-alpha)/psiss)*((1/beta)- (1-delta))-alpha)*psiss))^(1/(1-alpha))))^(alpha))/((((alpha/(1/beta -(1-delta)-(((1-alpha)/psiss)*((1/beta)- (1-delta))-alpha)*psiss))^(1/(1-alpha))))^(alpha) - delta*(((alpha/(1/beta -(1-delta)-(((1-alpha)/psiss)*((1/beta)- (1-delta))-alpha)*psiss))^(1/(1-alpha))))))^(1/(1+chi));
% k = n*(alpha/(1/beta - (1-delta) - (((1-alpha)/psiss)*((1/beta)- (1-delta))-alpha)*psiss))^(1/(1-alpha));
% psi = psiss;
% A = 0;
% r = (1/beta) - 1;
% y = k^alpha*n^(1-alpha);
% i = delta*k;
% c = y - i;
% lambda = 1/c;
% w = psiss*((alpha/(1/beta -(1-delta)-(((1-alpha)/psiss)*((1/beta)- (1-delta))-alpha)*psiss))^(1/(1-alpha)));
% mu = (((1-alpha)/psiss)*((1/beta)- (1-delta))-alpha);
