function [ss,res,exitflag]=steadystate_nontradable_lunch(beta,theta,chi,psi,delta,m,nu,nun,omega,mu,mus,sigma_hf,sigma_nt,gy,sigmat_ss, sigman_ss,eta)

ss0(1,1)=0.3;%l
ss0(1,2)=0.3*(1-eta);%lt
ss0(1,3)=0.3;%ls
ss0(1,4)=0.3*(1-eta);%lts
ss0(1,5)=1.13;%w
ss0(1,6)=1.13;%ws

para(1,1)=beta;
para(1,2)=theta;
para(1,3)=chi;
para(1,4)=psi; 
para(1,5)=delta; 
para(1,6)=m; 
para(1,7)=nu; 
para(1,8)=omega; 
para(1,9)=mu;
para(1,10)=mus; 
para(1,11)=sigma_hf; 
para(1,12)=gy;
para(1,13)=sigmat_ss;
para(1,14)=nun;%
para(1,15)=sigma_nt;%
para(1,16)=sigman_ss;
para(1,17)=eta;
ztss=1;
ztsss=1;
znss=1;
znsss=1;

options=optimoptions('fsolve','algorithm','Levenberg-Marquardt');
[ss1,res,exitflag]=fsolve(@steadystate_nontradable_fun,ss0,options,para);


l=ss1(1,1);
lt=ss1(1,2);
ls=ss1(1,3);
lts=ss1(1,4);
w=ss1(1,5);
ws=ss1(1,6);


int=1/beta-1;
lnon=l-lt;
lns=ls-lts;
c=(w/chi*(1-l)^theta)^(1/theta);
cs=(ws/chi*(1-ls)^theta)^(1/theta);
rer=(cs/c)^(-theta);
kr=1/beta-1+delta;
kt=w/kr*nu/(1-nu)*lt;
kts=ws/kr*nu/(1-nu)*lts;
kn=w/kr*nun/(1-nun)*lnon;
kns=ws/kr*nun/(1-nun)*lns;
k=kt+kn;
ks=kts+kns;
inv=delta*(k);
invs=delta*(ks);
%g=gy*y;
y=m*(c+inv)/(1-m*gy);
%gs=gy*y;
ys=(1-m)*(cs+invs)/(1-(1-m)*gy);
ph=sigmat_ss/(sigmat_ss-1)*kr^nu*w^(1-nu)/(ztss^(1-nu)*nu^nu*(1-nu)^(1-nu));%w
pfs=sigmat_ss/(sigmat_ss-1)*kr^nu*ws^(1-nu)/(ztsss^(1-nu)*nu^nu*(1-nu)^(1-nu));%ws
phs=ph/rer;
pf=pfs*rer;
pt=(mu*(ph)^(1-sigma_hf)+(1-mu)*(pf)^(1-sigma_hf))^(1/(1-sigma_hf));
pts=((1-mus)*(phs)^(1-sigma_hf)+mus*(pfs)^(1-sigma_hf))^(1/(1-sigma_hf));
pn=sigman_ss/(sigman_ss-1)*kr^nun*w^(1-nun)/(znss^(1-nun)*nun^nun*(1-nun)^(1-nun));%w
pns=sigman_ss/(sigman_ss-1)*kr^nun*ws^(1-nun)/(znsss^(1-nun)*nun^nun*(1-nun)^(1-nun));%ws
yn=eta*pn^(-sigma_nt)*y;
yt=(1-eta)*pt^(-sigma_nt)*y;
yns=eta*pns^(-sigma_nt)*ys;
yts=(1-eta)*pts^(-sigma_nt)*ys;
yh=mu*(ph/pt)^(-sigma_hf)*yt;
yf=(1-mu)*(pf/pt)^(-sigma_hf)*yt;
yhs=(1-mu)*(phs/pts)^(-sigma_hf)*yts;
yfs=mu*(pfs/pts)^(-sigma_hf)*yts;

ss(1,1)=int;
ss(1,2)=c;
ss(1,3)=cs;
ss(1,4)=rer;
ss(1,5)=kr;
ss(1,6)=k;
ss(1,24)=kt;
ss(1,25)=kn;
ss(1,7)=ks;
ss(1,26)=kts;
ss(1,27)=kns;
ss(1,8)=inv;
ss(1,9)=invs;
ss(1,10)=y;
ss(1,11)=ys;
ss(1,12)=ph;
ss(1,13)=pfs;
ss(1,14)=phs;
ss(1,15)=pf;
ss(1,28)=pt;
ss(1,29)=pts;
ss(1,30)=pn;
ss(1,31)=pns;
ss(1,32)=yn;
ss(1,33)=yt;
ss(1,34)=yns;
ss(1,35)=yts;
ss(1,16)=yh;
ss(1,17)=yf;
ss(1,18)=yhs;
ss(1,19)=yfs;
ss(1,20)=l;
ss(1,36)=lt;
ss(1,37)=lnon;
ss(1,21)=ls;
ss(1,38)=lts;
ss(1,39)=lns;
ss(1,22)=w;
ss(1,23)=ws;