function [ys,check]=first1_steadystate(junk,ys);

alpha = get_param_by_name('alpha');
mu = get_param_by_name('mu');
x_ss = get_param_by_name('x_ss');
beta_i = get_param_by_name('beta_i');
beta_p = get_param_by_name('beta_p');
m_ibar = get_param_by_name('m_ibar');
ee_hbar = get_param_by_name('eeh_bar');
lev = get_param_by_name('lev');
delta = get_param_by_name('delta');
nk = get_param_by_name('nk');
phi = get_param_by_name('phi');
gamma = get_param_by_name('gamma');
gamma_p = get_param_by_name('gamma_p');
chi = get_param_by_name('chi');
pi_ss = get_param_by_name('pi_ss');


check=0;

r_n=1/beta_p;
r=r_n;
sp = ((1-alpha)*(1-mu) + x_ss - 1)/x_ss; %income share of patient household
si = (1-alpha)*mu/x_ss; %income share of impatient household
ci_Y = (1-beta_i-m_ibar*(beta_p-beta_i))*si/(1-beta_i-m_ibar*(beta_p-beta_i)+ee_hbar*m_ibar*(1-beta_p));
qhi_Y = ee_hbar*si/(1-beta_i-m_ibar*(beta_p-beta_i-ee_hbar*(1-beta_p)));
r_k=lev/beta_p;
k_Y=alpha/(x_ss*(r_k-(1-delta)));
k_n=1/nk;
cp_Y  = sp-(1-1/beta_p)*((1-gamma_p)*beta_p*m_ibar*qhi_Y+k_Y-(k_Y/k_n))/gamma_p;
qhp_Y = ee_hbar*cp_Y/(1-beta_p);
l_p=((1-mu)*(1-alpha)/(x_ss*cp_Y))^(1/(1+phi));
l_i=(mu*(1-alpha)/(x_ss*ci_Y))^(1/(1+phi));
y=(k_Y)^(alpha/(1-alpha))*(l_i)^mu*(l_p)^(1-mu);
bi_Y = beta_p*m_ibar*ee_hbar*si/(1-beta_i-m_ibar*(beta_p-beta_i) +ee_hbar*m_ibar*(1-beta_p));
b_Y= k_Y-k_Y/k_n;
k=k_Y*y;
n=nk*k;
c_e=(1-gamma)*n/gamma;
b_i= bi_Y*y;
c_p=cp_Y*y;
c_i=ci_Y*y;
b=b_Y*y;
w_p=1/x_ss*(1-mu)*(1-alpha)*y/l_p;
w_i=1/x_ss*mu*(1-alpha)*y/l_i;
d_p=((1-gamma_p)*b_i+b)/gamma_p;
lambda=(beta_p-beta_i)/c_i;
h_i=c_i*(1-beta_p)/((1-gamma_p)*c_i*(1-beta_p)+gamma_p*c_p*(1-beta_i-m_ibar*(beta_p-beta_i)));
h_p=(1-h_i*(1-gamma_p))/gamma_p;
q_h=qhi_Y*y/h_i;
I=delta*k_Y*y;
v=lev*r*n;
dbar=(n-gamma*v)/(1-gamma);
c=gamma_p*c_p+(1-gamma_p)*c_i+c_e;
G=y-c-I;
CC_SS= log(c);
I_SS= log(I);
QQ_SS= log(q_h);
Y_SS= log(y);

% data_Y = 0;
% data_I = 0;
% data_P = 0;
% data_hp = 0;
% data_r = 0;

xxx =[1
    b
    b_i
    c
    c_e
    c_i
    c_p
    d_p
    ee_hbar
    G
    h_i
    h_p
    I
    k
    lambda
    l_i
    l_p
    m_ibar
    n
    1
    q_h
    1
    r
    r_k
    r_n
    1
    v
    w_i
    w_p
    x_ss
    y
    ];
ys=log(xxx);
;

end