load param

infls=Pis;
Rs=infls/betal;
Lams=betal;
Rrs=Rs;
ss=1;
pmFs=ss;
Rfs=sp+Rs;
Qfs=1/(Rfs-kappa);
Rbs=Delta*sp+Rs;
Qbs=1/(Rbs-kappa);
vprs=(1-Phip)*((1-Phip*inflF^(1-epsilon))/(1-Phip))^(epsilon/(epsilon-1))        /          (1-Phip*inflF^(-epsilon));
pFs=pmFs*(epsilon/(epsilon-1))*(1-Phip*inflF^(1-epsilon)*Lams*((1-Phip*inflF^(1-epsilon))/(1-Phip)))...
    /     (       ( (1-Phip*inflF^(1-epsilon))   / (1-Phip)   )^(1/(epsilon-1))     -     Phip*inflF^(1-epsilon)*Lams*(  (1-Phip*inflF^(1-epsilon))  /  (1-Phip)  )      );
pFrs=pFs*(   (1-Phip*inflF^(1-epsilon))  /  (1-Phip)  )^(1/(1-epsilon));
pHs=gamma*(pFs/(1-gamma))^((gamma-1)/gamma);
Rhms=Rs*exp(-zeta*ss);
Xs=(pHs/ss)^v;
EXs=Xs*(pHs/ss)^(-v);
hms=EXs* (fracYH*gamma/(1-gamma)*vprs-(1-fracYH)*pmFs) / (fracYH*gamma/(1-gamma)*vprs*Rhms);
YmFs=(EXs-Rhms*ss*hms)/pmFs;
YFs=vprs*YmFs;
x1rs= (YFs*(1-Phip*inflF^(1-epsilon))/(1-Phip))        /       (1-Phip*inflF^(1-epsilon)*Lams* (1-Phip*inflF^(1-epsilon)) / (1-Phip) );
x2rs=(epsilon-1)/epsilon*x1rs;
YhHs=gamma*YFs/(1-gamma);
YHs=YhHs+EXs;
Ys=(YhHs)^gamma*YFs^(1-gamma);
prs=( (1-Phip*inflH^(1-epsilon)) / (1-Phip) )^(1/(1-epsilon))*pHs;
vps=(1-Phip)*( (1-Phip*inflH^(1-epsilon)) / (1-Phip) )^(epsilon/(epsilon-1))  /  (1-Phip*inflH^(-epsilon));
x1s=YHs*(prs/pHs)^(1-epsilon)/(1-Phip*inflH^(1-epsilon)*Lams*(prs/pHs)^(1-epsilon));
x2s=x1s*(epsilon-1)/epsilon;
pms=x2s*pHs/YHs*( (prs/pHs)^epsilon -Phip*inflH^(-epsilon)*Lams );
Yms=vps*YHs;
theta1s=Lams/(infls*Qfs*(1-kappa*Lams/infls))-1;
theta2s=theta1s+1;
Ks=alpha*Lams*pms*Yms/(theta2s*(1-Lams*(1-delta)));
Ls=(Yms/(As*Ks^alpha))^(1/(1-alpha));
ws=(1-alpha)*pms*Yms/Ls;
Is=delta*Ks;
fs=Is/(Qfs*(1-kappa/infls));
Cs=Ys-Is;
Cls=fracCl*Cs;
lamls=1/Cls;
Cbs=Cs-Cls;
lam1s=1/Cbs;
lam2s=lam1s*(Qbs-betab*(1+kappa*Qbs)/infls)  /  (Qbs*(1-kappa*betab/infls));
%lam2s=lam1s*(Qbs-betab*(1+kappa*Qbs)/infls);
qhs=gammah/(lamls*(1-betal))+gammah/(lam1s*(1-betab)-xi*lam2s);
Hls=gammah/(qhs*lamls*(1-betal));
Hbs=gammah/(qhs*(lam1s*(1-betab)-xi*lam2s));
bs=xi*qhs*Hbs/(Qbs*(1-kappa/infls));
Phi=(Ls/  (  (ws*lamls)^(1/eta)  +  (ws*lam1s)^(1/eta)  )  )^(-eta);
Lbs=(lam1s*ws/Phi)^(1/eta);
Lls=(lamls*ws/Phi)^(1/eta);
ns = (Qfs*fs + Qbs*bs )/levs;
phis = (Qfs*fs + Delta*Qbs*bs)/ns;
M=ns*(1-sigma*Rs/infls)-sigma/infls*((Rfs-Rs)*Qfs*fs+(Rbs-Rs)*Qbs*bs-Rhms*ss*hms);
ds=Qfs*fs+Qbs*bs-ns-ss*hms;
theta = ( (1-sigma)*Rs*Pis^(-1)*Lams + (1-sigma)*(Rfs-Rs)*phis*Pis^(-1)*Lams)/(phis - sigma*(Rfs - Rs)*phis^2*Pis^(-1)*Lams - sigma*phis*Rs*Pis^(-1)*Lams);
Omegas = 1-sigma + sigma*theta*phis;
lams=Lams*Omegas/infls*(Rfs-Rs)/(theta-Lams*Omegas/infls*(Rfs-Rs));