function FF =solve_ss1(epsilon,pi_ss,beta,taut,sigma,m1,delta,alpha,xi,theta,omega1,omega2)

%% 

q_ss = 1; % eq.17
pw_ss = (epsilon-1)/epsilon; % eq.9
R_ss = pi_ss/beta; % eq.4

F_omega1 = normcdf((log(omega1)+sigma^2/2)/ sigma); 
g_omega1 = 1-normcdf((log(omega1)-sigma^2/2)/sigma)-(1-F_omega1)*omega1;
f_omega1 = 1/(omega1*sigma*sqrt(2*3.1416))*exp(-(log(omega1)+sigma^2/2)^2/(2*sigma^2));
int_omega1 = normcdf((log(omega1)-sigma^2/2)/sigma);
F_omega2 = normcdf((log(omega2)+sigma^2/2)/ sigma); 
int_omega2 = normcdf((log(omega2)-sigma^2/2)/sigma);

Lambda_ss=(1-F_omega1)/(1-F_omega1-m1*omega1*f_omega1); % eq.6
Rk_ss = -Lambda_ss*(R_ss*taut*F_omega2-R_ss)/(g_omega1+Lambda_ss*omega1*(1-F_omega1)+Lambda_ss*(1-m1)*(int_omega1-int_omega2));%eq.5
rk_ss=Rk_ss +delta-1; % eq.14
K_Y=(pw_ss*(1-alpha))/rk_ss; % eq.13
nb_Y = K_Y; %eq.15
n_Y = xi * Rk_ss*nb_Y*g_omega1+pw_ss*(1-theta)*alpha; % eq.18
b_Y = nb_Y - n_Y;


FF=[omega1*Rk_ss*nb_Y*(1-F_omega1)+(1-m1)*Rk_ss*nb_Y*(int_omega1-int_omega2)+R_ss*taut*b_Y*F_omega2-R_ss*b_Y;
    R_ss*taut*b_Y/((1-m1)*Rk_ss*nb_Y)-omega2];

end







