function [ys, check] = usmodel2_steadystate(ys,exe)
 global M_

	   
	  
gammaa=M_.params(1);
Xiur = M_.params(2);
Cur = M_.params(3);
chiwu = M_.params(4);
pia=M_.params(5);
DiscF=M_.params(6);
zetaa=M_.params(7);
h=M_.params(8);
sigmau=M_.params(9);
sigmar=M_.params(10);
dzetap=M_.params(11);
omegau=M_.params(12);
nu=M_.params(13);
zetaw=M_.params(14);
zetap=M_.params(15);
phiT=M_.params(16);
rhor=M_.params(17);
phipi=M_.params(18);
phiy=M_.params(19);
rhoz=M_.params(20);
rhomu=M_.params(21);
rhob=M_.params(22);
rhophi=M_.params(23);
rhoB=M_.params(24);
rhozeta=M_.params(25);
rhog=M_.params(26);
sigmaz=M_.params(27);
sigmalambdaf=M_.params(28);
sigmamu=M_.params(29);
sigmab=M_.params(30);
sigmaphi=M_.params(31);
sigmaB=M_.params(32);
sigmaT=M_.params(33);
sigmam=M_.params(34);
sigmazeta=M_.params(35);
sigmag=M_.params(36);
gamma = M_.params(37);
piss = M_.params(38);
zetass = M_.params(39);
betau = M_.params(40);
BondDuration = M_.params(41); 
dzeta = M_.params(42);
 alpha=M_.params(43);
lambdaw=M_.params(44);
lambdafpara = M_.params(45);
rholambdaf=M_.params(46);
delta=M_.params(47);
BLB =M_.params(48);
Gzss=M_.params(49);
Bzss=M_.params(50);
Piss = M_.params(51);
Rss = M_.params(52);
RLss = M_.params(53);
kappa = M_.params(54);
kappaEH = M_.params(55);
betar = M_.params(56);
wztilde = M_.params(57);
Kztilde = M_.params(58);
Ltilde = M_.params(59);
betabar = M_.params(60);
rkss = M_.params(61);
wzss = M_.params(62);
Kzss = M_.params(63);
Kzbarss = M_.params(64);
Lss = M_.params(65);
Izss = M_.params(66);
Czss = M_.params(67);
Czrss = M_.params(68);
Czuss = M_.params(69);
chipu = M_.params(70);
qu = M_.params(71);
Tzss = M_.params(72);
d2s=M_.params(73);
d2a=M_.params(74);
bpara=M_.params(75);
sigmaa=M_.params(76);
BLMVzss = M_.params(77);
phirss = M_.params(78);
buss = M_.params(79);
phiuss = M_.params(80);
BLMVB= M_.params(81);
Yzss = M_.params(82);
uss = M_.params(83);
qss = M_.params(84);
muss = M_.params(85);
mcss = M_.params(86);
Xiuss = M_.params(87);
Xirss = M_.params(88);
Xpnuss = M_.params(89);
Xpnrss = M_.params(90);
Xpduss = M_.params(91);
Xpdrss = M_.params(92);
Xwnuss = M_.params(93);
Xwnrss = M_.params(94);
Xwduss = M_.params(95);
Xwdrss = M_.params(96);
lambdafss = M_.params(97);
brss = M_.params(98);
BLzss = M_.params(99);
BTotMVzss = M_.params(100);

check = 0;

br = 0;
 bu = 0;
BLMVz = 0;
BLz = 0;
BTotMVz = 0;
 Bz = 0;
 Gz  = 0;
 lambdaf  =0;        
rk = 0;
mc =0;
Kz = 0;
L  = 0;
Kzbar =0;
Iz = 0;
Czr = 0;
Czu = 0;
Xiu =0;
Xir = 0;
Xpnu = 0;
Xpnr =  0;
Xpdu = 0;
Xpdr = 0;
phiu = 0;
phir = 0;
Xwnu = 0;
Xwnr = 0;
Xwdu = 0;
Xwdr = 0;
 pi  = 0;
 q  = 0;
 R  = 0;
 RL =0;
 Tz  = 0;
 u  =0;
 Xim  = 0;
 XiT  = 0;
 Xizeta  = 0;
Yz  =0;
z=0;
zeta = 0;
wz = 0;
mu = 0;
XiB  = 0;


 
	  
ys = [br
bu
BLMVz 
BLz 
BTotMVz 
Bz 
Gz 
lambdaf 
rk 
mc 
Kz 
L 
Kzbar 
Iz 
Czr 
Czu 
Xiu 
Xir 
Xpnu 
Xpnr 
Xpdu 
Xpdr 
phiu 
phir 
Xwnu 
Xwnr 
Xwdu 
Xwdr 
pi 
q
R 
RL 
Tz 
u 
Xim 
XiT 
Xizeta 
Yz 
z 
zeta 
wz 
mu 
XiB];