function y=nss(xx,zeta,exss,delta_X,theta_ng,alpha,nu,k_ng, L_ng,d0,d1,d2)

y_ng=xx(1);

mu=0;
e=zeta*(1-mu)*y_ng;
ex=exss*e;
x=(e+ex)/(1-delta_X);

f1=y_ng-(k_ng^alpha*L_ng^(1-alpha))/(1+(d0+d1*x+d2*x^2));

y=[f1];