function F=P4_ss(x)
eta_c=0.6;          
eta_n=4;           
eta_x=0.01;
eta_h=0.02;         
nu_cash=0.5;              
beta=0.985;        
lambda_HD=0.5;       
lambda_FD=0.5;     
eta=2;             
alpha=1/3; 
alpha_F=1/3;
theta_K=18;
theta_prod1=10;
theta_prod2=10;
delta=0.02;%0.02        
mu_H_ss=0.2;          
mu_F_ss=0.2;           
psi_gov=x(10);                                  
theta_0_HF=0.2;
theta_0_FH=0.2;
theta_HB_0=0.1;%0.0001
theta_FB_0=0.1;%
i_R_ss=1/beta-1;
eps_1=1.5;
eps_2=0.2;
eps_q1=0.02;%0.02
eps_q2=0.2;
kappa=0.2;
Y_H_ss=1;
Y_F_ss=1;
b=x(11);
sub=0.99997563;%



F=[x(1)-1/beta+1;%i_R
x(2)-(1+x(18))*(x(1)+delta); %r_H
x(3)-1+(eta_n*(0.3)^(1/eta_c))/x(4); %N_H
x(5)-eta_x*nu_cash*(0.3)^(1/eta_c)*(1+x(1))/(x(1)); %m_HP
x(6)-eta_x*(1-nu_cash)*(0.3)^(1/eta_c)*(1+x(1))/(x(1)-x(15)); %d_H
(x(7)/x(3))-(alpha/(1-alpha))*(x(4)/x(2)); 
2-(x(7))^alpha*(x(3))^(1-alpha); 
x(8)-delta*x(7);
x(9)-x(5)-x(8);
(1-x(10))-0.3-x(8);
x(8)-x(6)+0.5*x(11)-x(9);
x(12)-(1/(1-beta))*eta_h*(0.3)^(1/eta_c);
x(13)-(kappa*(x(12)/x(8)))^eps_q1;
1+x(14)-(1+x(18))/x(13);
x(15)-(1-mu_H_ss)*x(18);
x(16)-x(8)+(1-mu_H_ss)*x(6);
x(17)-theta_HB_0*(x(16)/x(6));
%1+x(18)-(1+x(1))*(1+x(17));];  
x(18)-x(1)*(1+x(17));];  


  
% Columns 1 through 7

%   0.0152    0.0358    0.6905   1.7378    0.0448    0.2282    16.7762     

%    i_R       r_H       N_H      w_H       m_HP      d_H       K_H

% Columns 8 through 14

%    0.3355    0.3803    0.3645    0.5459     0.1793    0.9779    0.0383   

%     l_H       m_H      psi_gov     b         z_HA       q_H      i_HL

% Columns 15 through 21

%    0.0122    0.1530    0.0001     0.0153     

%     i_HD      l_HB     theta_HB    i_HC