function n = BenchmarkModel_steadystate_helper(n0,VARTHETA,BETA,TAO_SS,w,PSI,DELTA,r_k, MU, g, ETA);
        options = optimset('Display','off','TolX',1e-10,'TolFun',1e-20);
       
        n = fsolve(@(n) n-log((-(1-VARTHETA)/(1-BETA*VARTHETA)/(1-VARTHETA)*(1-BETA*VARTHETA)*(1-TAO_SS)*exp(w)/PSI/((DELTA*exp(w)*MU/(r_k*(1-MU)))-exp(0)*(exp(w)*MU/(r_k*(1-MU)))^MU+exp(g)/exp(n)))^(1/(1+ETA))), n0,options);
      
end