%y=zeros(43,1);
function [z,check] = FINAL_steadystate(junk,ys)

z=zeros(159,1);

global ex_ ex_det_ it_ recur_
global cadjc1 cadjc2 calfa1 calfa2 cbeta1 cbeta2 ccs1 ccs2 cexpbeta1 cexpbeta2 
global cfc1 cfc2 chabb1 chabb2 cim1 cim2 cimpbeta1 cimpbeta2 cindp1 cindp2 
global cindpc1 cindpc2 cindpm1 cindpm2 cindw1 cindw2 cinvs1 cinvs2 clandaw1 clandaw2 
global cmap1 cmap2 cmaq1 cmaq2 cmaw1 cmaw2 cmc1 cmc2 cmy1
global cmy2 coilc1 coilc2 coilm1 coilm2 coily1 coily2 constelab1 constelab2 
global constepinf1 constepinf2 constepminf1 constepminf2 constepoil conster1 conster2 cprobe cprobp1 cprobp2 
global cprobpc1 cprobpc2 cprobpm1 cprobpm2 cprobw1 cprobw2 crdpi1 crdpi2 crdy1 crdy2 
global crhoa1 crhoa2 crhob1 crhob2 crhog1 crhog2 crhoms1 crhoms2 crhont1 crhont2 
global crhopinf1 crhopinf2 crhopm1 crhopm2 crhopoil crhoqs1 crhoqs2 crhos crhow1 crhow2 
global crpi1 crpi2 crr1 crr2 cry1 cry2 csa1 csadjcost1 csadjcost2 csigl1 
global csigl2 csigma1 csigma2 csnt1 cspm1 cspm2 cspoil csrp csubst1 csubst2
global csw1 csw2 ctou1 ctou2 ctrend cydc1 cydc2 czcap1 czcap2 duma ctrend1 ctrend2


z(105)=constepoil;
z(23)=ctrend1;
z(9)=ctrend1;
z(16)=ctrend1;
z(21)=ctrend1;
z(106)=constepinf1;
z(110)=constepinf1+conster1;
z(15)=constelab1;
z(94)=constepminf1;
z(78)=constepinf1;
z(24)=ctrend2;
z(10)=ctrend2;
z(17)=ctrend2;
z(22)=ctrend2;
z(107)=constepinf2;
z(111)=constepinf2+conster2;
z(18)=constelab2;
z(95)=constepminf2;
z(79)=constepinf2;

check=0;
