


function F = ssfun(x,OMEGAH,EPS,ZETA3,BETA,DH)

% Function to compute F in fsolve for seven central steady-state variables

% In the process, we make the following notation:

% x(1): HTHILDE     (TOTAL HOUSING, IMPATIENT)
% x(2): H1          (OWNER HOUSING, IMPATIENT)
% x(3): HZ          (RENTAL HOUSING, IMPATIENT) 


% HTHILDE
cell1 = ( OMEGAH^(1/EPS) * x(2)^((EPS-1)/EPS) + (1-OMEGAH)^(1/EPS)* x(3)^((EPS-1)/EPS) )^(EPS/(EPS-1)) - x(1);

% H1 
cell2 = (1-BETA*(1-DH))*(OMEGAH/(1-OMEGAH))^(1/EPS)*(x(3)/x(2))^(1/EPS)-x(1);

% HZ - All of impatient housing (HTHILDE) must be the sum of rental and
% impatient housing.

cell3 = x(3)+x(2)-x(1);

F = [cell1 ;
     cell2 ;
     cell3] ;     