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

global M_

paramfile_deficit

nparams = size(M_.param_names,1);
for icount = 1:nparams
    eval(['M_.params(icount) = ',M_.param_names(icount,:),';'])
end


check=0;


l=1;
B=0;
Delta=1;
y=1;
h=1;
pi=piss;
w=alph;
c_s=0.447214;
c_e=0.552786;
r_d=r_dss;
r_b=r_d;
i=iss;
q_h=0.6;
b=0.326700003;
j_b=0.0297;
k=0.029403;
d=0.297297002;
%G=0.2;
T=G;
t_p=0.200297;
lb=(1-beta_e*r_b)/c_e^gamma_e;


ys = [ l w c_s c_e d b T r_d r_b pi y t_p q_h a eps_h gfd k i B j_b h G lb Delta] 



