function yf=solve_yf_ss(yf0,tauf_la,alphaf_m,muf_la,pf_m,m_f,nuf,mcf,wf,alphaf,muf_l,phif,deltaf_k,rf_k,tauf_k,pf_d)
         options=optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
         yf=fsolve(@(yf)(1-tauf_la)*((1-alphaf_m)*muf_la*pf_m*m_f)*((nuf*mcf/wf+(1-alphaf_m)*(1-alphaf-nuf)*muf_l*mcf/wf)*yf+(1-alphaf_m)*muf_la*pf_m/wf*m_f)^phif-((alphaf*deltaf_k*mcf/rf_k)+tauf_la*(nuf*mcf+(1-alphaf_m)*(1-alphaf-nuf)*muf_l*mcf)+(tauf_k*alphaf*mcf)-((1-alphaf_m)*muf_la*(1-alphaf-nuf)*mcf)-pf_d)*yf*((nuf*mcf/wf+(1-alphaf_m)*(1-alphaf-nuf)*muf_l*mcf/wf)*yf+(1-alphaf_m)*muf_la*pf_m/wf*m_f)^phif-(1-tauf_la)*wf,yf0,options);
end