function n = solve_n_ss(alpha,omega,beta,gamma,phi,rho,delta,a_ss,B,chi);
n_ss=0;
func_hand = @employ;   
options = optimset('Display','off','TolFun',1e-50);
n = fsolve(func_hand, 0.5, options, alpha, omega, beta, gamma, phi, rho, delta, a_ss, B, chi);