function z = SS_equations_R1(x,delta,omega,kappa,beta,theta,phi,alphaN,alphaT,rho,psii,bbar)
% 3 endogenous variables to minimize over
ctss = x(1);
htss = x(2);
cnss = x(3);

anss = 1; % SS equation 13
atss = 1; % SS equation 14
lamss = 0; % SS equation 16

iss = (1/beta)-1; % equation 3
css = ((omega^(1/kappa)*ctss^((kappa-1)/kappa))+((1-omega)^(1/kappa)*cnss^((kappa-1)/kappa)))^(kappa/(kappa-1)); % equation 9
pnss = (((1-omega)*ctss)/(omega*cnss))^(1/kappa); % equation 2
wss = ((1-alphaT)*atss*(htss^(-alphaT)))/(1+iss*theta); % euqation 8
hnss = ((wss*(1+iss*theta))/((1-alphaN)*pnss*anss))^(-1/alphaN); % euqation 7
ytss = atss*htss^(1-alphaT); % euqation 6
ynss = anss*hnss^(1-alphaN); % equation 5
hss = htss+hnss; % equation 10
uss = ((css-((hss^delta)/delta))^(-rho))*((omega*(css/ctss))^(1/kappa)); % equation 1
bss = (ctss-ytss)/iss; % equation 11

z(1) = (omega*(css/ctss))^(1/kappa)*wss*(1+((1-phi)/phi)*(lamss/uss))-hss^(delta-1); % equation 4
z(2) = ynss-cnss; % equation 12
z(3) = iss - (1/beta) + 1 - psii*(exp(bss-bbar)-1); % equation 15