
function [ys] = STEADYBtesi()
	

global M_

alpha= M_.params(9);
rhop = M_.params(1);
rhog = M_.params(2);
pi=    M_.params(7);
etap=  M_.params(11);
etag=  M_.params(12);
delta= M_.params(5);
fp0=   M_.params(16);
fg0=   M_.params(17);

l0=0.781;
g0=0.153;
u0=1-l0-g0;
hp0=rhop*l0;
hg0=rhog*g0;
frac=hp0/hg0;
D=delta.*(1-rhop);
cost=0.32;
qp0=0.35;
qg0=0.224;
vp0=hp0/qp0;
vg0=hg0/qg0;
costr=(1-D).*cost./qp0;
omega0=(alpha.*l0.^(alpha-1))./(1+costr);
w0=omega0;
wg0=pi*omega0;
r0=cost*omega0;
xp0=w0/(1-delta*(1-rhop-fp0));
xg0=wg0/(1-delta*(1-rhog-fg0));
sigma0=(1+frac*(xp0/xg0))^(-1);
y0=l0^alpha;
c0=y0-rhop*l0*(r0/qp0)-rhog*g0*(r0/qg0);
betat0=delta;
mup0=hp0*(((1-sigma0)*u0)^(-etap))*((vp0)^(etap-1));
mug0=hg0*((sigma0*u0)^(-etag))*((vg0)^(etag-1));

	
ys =[1;u0;vp0;vg0;hp0;hg0;xp0;xg0;sigma0;l0;c0;g0;betat0;qp0;omega0;mup0;mug0;r0];

