function res=steady_stateNS(X, parameters)

global C	N	B_P	Lambda	R	W_P	Pi	D_P	Y	Pstar_P	s	MC_P	Gamma1	Gamma2	eps_a	eps_r;

C       =	X(1);
N       =	X(2);
B_P     =	X(3);
Lambda	=	X(4);
R       =	X(5);
W_P     =	X(6);
Pi      =	X(7);
D_P     =	X(8);
Y       =	X(9);
Pstar_P	=	X(10);
s       =	X(11);
MC_P	=	X(12);
Gamma1	=	X(13);
Gamma2	=	X(14);
eps_a	=	X(15);
eps_r	=	X(16);

nu      =parameters(1);
upsilon	=parameters(2);
beta	=parameters(3);
theta	=parameters(4);
epsilon	=parameters(5);
phi     =parameters(6);
sigma_a	=parameters(7);
sigma_r	=parameters(8);
rho_a	=parameters(9);
rho_r	=parameters(10);


%% Steady State Equations
res(1,1)	= W_P - C*N^upsilon;
res(2,1)	= R - 1/beta;
res(3,1)	= beta - Lambda;
res(4,1)	= D_P - Y/epsilon;
res(5,1)	= MC_P - W_P;
res(6,1)	= MC_P - (epsilon-1)/epsilon;
res(7,1)	= Gamma1 - ((epsilon-1)/epsilon)*Y/(1-theta*beta);
res(8,1)	= Gamma2 - Y/(1-theta*beta);
res(9,1)	= Pstar_P - 1;
res(10,1)	= s - 1;
res(11,1)	= Pi - 1;
res(12,1)	= Y - C;
res(13,1)	= Y - N;
res(14,1)	= B_P;
res(15,1)	= C - W_P*N - D_P;
res(16,1)	= eps_a - 1;
res(17,1)	= eps_r - 1;

res = res';

end