load param

Lams = beta;
infls=Pis;
Rds=Pis/beta;
RFs=Rds+sp;
RBs=Delta*(RFs-Rds)+Rds;
Qs=1/(RFs-kappa);
QBs=1/(RBs-kappa);
prs=((1-phip*Pis^((gammap-1)*(1-epsilon)))        /(1-phip)       )^(1/1-epsilon);
pms=prs*(epsilon-1)/epsilon;
us=1-Ls;
Ks=Ls*(  ( 1-Lams*(1-delta) )        /   alpha*pms*As*Qs*Pis*(1-kappa*Lams/Pis)            )^(1/(alpha-1));
Yms=As*Ks^alpha*Ls^(1-alpha);
Is=delta*Ks;
pks=1;
vps=(1-phip)*prs^(-epsilon)     /     (1-phip*Pis^((1-gammap)*epsilon));
Ys=Yms/vps;
x1s=pms*Ys/(1-phip*Lams);
x2s=Ys/(1-phip*Lams);
fms=pks*delta*Ks/(Qs*(1-kappa/Pis));
v1s=Lams/(Qs*Pis-Lams*kappa*Qs);
v2s=v1s-1;
Ihs=Is;
Gs=gy*Ys;
Cs=Ys-Is-Gs;
mus=(1-bh*beta)/(Cs*(1-bh));
Gammas=deltaN*Ls/Qgammas;
Us=1-(1-deltaN)*Ls;
mum=log(Us/Gammas)/log(Qgammas/mubar);
Ms=mubar*Us^(mum)*Gammas^(1-mum);
QUs=mubar*Us^(mum-1)*Gammas^(1-mum);
Thetas=xi/Qgammas;
ws=(1-alpha)*pms*Ys/Ls-Thetas+Lams*Thetas*(1-deltaN);
B=Bw*ws;
chi=((ws-eta*(1-alpha)*pms*Ys/Ls-eta*(1-deltaN)*Lams*QUs*Thetas)/(1-eta)-B)/(Ls^eta/mus);
RTRs=Rds;
Rres=Rds;
fcbs=fracf*Ys/Qs;
bcbs=fracb*Ys/QBs;
res=Qs*fcbs+QBs*bcbs;
fs=fms-fcbs;
bG=bgy*Ys/QBs;
bs=bG-bcbs;
ns=((RFs-Rds)*Qs*fs+(RBs-Rds)*QBs*bs+(Rres-Rds)*res)/(Pis-Rds);
phis=(Qs*fs+Delta*QBs*bs)/ns;
thetas=((1-sigma)*(Lams/Pis)*(Rds+res/ns*(Rres-Rds)+phis*(RFs-Rds)))      /        (phis*(1-Lams/Pis*(res/ns*(Rres-Rds)+RFs)));
Omegas=1-sigma+sigma*thetas*phis;
lams=(Lams*Omegas/Pis*(RFs-Rds)) /   (thetas-Lams*Omegas/Pis*(RFs-Rds));
omegas=-(1+lams)*Lams*Omegas/Pis*(Rres-Rds);
Ts=Gs+bG/Pis-QBs*bG*(1-kappa/Pis);
ds=Qs*fs+QBs*bs+res-ns;
levs=res/ds;
As=1;
divs=ns-ns*Rds/Pis+Ys-ws*Ls-fms/Pis+fms*Qs*(1-kappa/Pis)-xi*Gammas-Is;
X=(1-Ls)*B+ns*(1-Rds/Pis)+Is+Gs+fms*(1/Pis+Qs*(1-kappa/Pis))+(Rds-1)*ds+Ts-xi*Gammas;

