load param

Lams = beta;
Rds = Pis/beta;
RFs = (1+sp)^(1/4)*Rds;
RBs = Deltas*(RFs - Rds) + Rds;
QFs = (RFs - kappaf)^(-1);
QBs = (RBs - kappab)^(-1);
RLFs = QFs^(-1) + kappaf;
RLBs = QBs^(-1) + kappab;
QEHs = (Rds - kappab)^(-1);
REHs = QEHs^(-1) + kappab;
M2s = Lams*Pis^(-1)/(QFs*(1-kappaf*Lams*Pis^(-1)));
Ms = (M2s-1)*psi+1;
pms = ((epsip-1)/epsip); % steady state marginal cost
Ks = (alpha*pms/(Ms/Lams - (1-delta0)*Ms))^(1/(1-alpha));
Yms = Ks^(alpha);
ws = pms*(1-alpha)*Ks^(alpha);
Ihs = delta0*Ks;
Is = Ihs;
Ys = Yms;
Gs = gy*Ys;
Cs = Ys - Is - Gs -rotemberg/2*(Pis-1)^2*Ys;
wss = (ws^(1-epsiw)*(1-phiw*Pis^((1-epsiw)*(gammaw-1)))/(1-phiw))^(1/(1-epsiw)); % reset wage
mrss = ((epsiw-1)/epsiw)*wss*(1-phiw*Lams*Pis^(epsiw*(1-gammaw)))/(1-phiw*Lams*Pis^((1-epsiw)*(gammaw-1)));
vws = (1-phiw)*(wss/ws)^(-epsiw)/(1-phiw*Pis^(epsiw*(1-gammaw)));
Ls = vws;
mus = (1/Cs)*(1/(1-bh) - beta*bh/(1-bh));
chi = mus*mrss*Ls^(-eta);
delta1 = pms*alpha*Ks^(alpha-1)/Ms;
f1s = mrss*ws^(epsiw)/(1-phiw*Lams*Pis^(epsiw*(1-gammaw)));
f2s = ws^(epsiw)/(1-phiw*Lams*Pis^((1-epsiw)*(gammaw-1)));
fprs = psi*Ihs/(QFs*(1-kappaf*Pis^(-1)));
bGs = bgy*Ys/QBs;
bcbs = Ys*fracb/QBs; 
fcbs = Ys*fracf/QFs;
bs = bGs - bcbs;
res = QFs*fcbs + QBs*bcbs;
fs = fprs - fcbs;
ns = (QFs*fs + QBs*bs + res)/levs;
ds = QFs*fs + QBs*bs + res - ns;
phis = (QFs*fs + Deltas*QBs*bs)/ns;
thetas = ( (1-sigma)*Rds*Pis^(-1)*Lams + (1-sigma)*(RFs-Rds)*phis*Pis^(-1)*Lams)/(phis - sigma*(RFs - Rds)*phis^2*Pis^(-1)*Lams - sigma*phis*Rds*Pis^(-1)*Lams);
X = ns*(1-sigma*((RFs-Rds)*phis + Rds)*Pis^(-1));
as = thetas*phis;
Omegas = 1-sigma + sigma*as;
lams = ((thetas/((RFs - Rds)*Lams*Omegas*Pis^(-1))) - 1)^(-1);
Res = Rds;
omegas = (Rds - Res)*(1+lams)/Pis*Lams*Omegas;
Tcbs = (1+kappaf*QFs)*Pis^(-1)*fcbs + (1+kappab*QBs)*Pis^(-1)*bcbs - Res*Pis^(-1)*res;
Ts = Gs + bGs*Pis^(-1) - Tcbs - QBs*(bGs - kappab*Pis^(-1)*bGs);