% used for second revision of paper

$$$$$$$ model equations

pd1(t) = GM*( pd1(t-1) + BK*pd1id(t-1) ) + (1-GM)*Xd1(t)                                         
pd2(t) = GM*( pd2(t-1) + BK*pd2id(t-1) ) + (1-GM)*Xd2(t) 

px1(t) = GM*( px1(t-1) + BK*px1id(t-1) ) + (1-GM)*Xx1(t)                                         
px2(t) = GM*( px2(t-1) + BK*px2id(t-1) ) + (1-GM)*Xx2(t)  

pn1(t) = GM*( pn1(t-1) + BK*pn1id(t-1) ) + (1-GM)*Xn1(t)                                         
pn2(t) = GM*( pn2(t-1) + BK*pn2id(t-1) ) + (1-GM)*Xn2(t)

wg1(t) = GW*( wg1(t-1) + BW*wg1id(t-1) ) + (1-GW)*Xw1(t)                                         
wg2(t) = GW*( wg2(t-1) + BW*wg2id(t-1) ) + (1-GW)*Xw2(t)

Kt1(t) = (1-DP)*Kt1(t-1) + DP*It1(t)
Kt2(t) = (1-DP)*Kt2(t-1) + DP*It2(t)

Kn1(t) = (1-DP)*Kn1(t-1) + DP*In1(t)
Kn2(t) = (1-DP)*Kn2(t-1) + DP*In2(t)

BX(t) = (1/BT)*BX(t-1) + YNB*(yn1(t) + pn1(t) - p1(t)) + YTB*(yt1(t) + pa1(t) - p1(t)) - CB*c1(t) - ITB*It1(t) - INB*In1(t) - GY*(G1(t) + pg1(t) - p1(t))

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Z1d(t) = Z1(t)
Z2d(t) = Z2(t)
Z3d(t) = Z3(t)
Z4d(t) = Z4(t)

int1L(t) = int1(t)
int2L(t) = int2(t)

p1d(t) = p1(t)
p2d(t) = p2(t)

py1d(t) = py1(t)
py2d(t) = py2(t)

c1d(t) = c1(t)
c2d(t) = c2(t)

pd1id(t) = pd1i(t)
pd2id(t) = pd2i(t)
px1id(t) = px1i(t)
px2id(t) = px2i(t)
pn1id(t) = pn1i(t)
pn2id(t) = pn2i(t)
wg1id(t) = wg1i(t)
wg2id(t) = wg2i(t)

ys1(t) = CY*ys1(t-1) + (1-CY)*y1(t)
ys2(t) = CY*ys2(t-1) + (1-CY)*y2(t)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

At1(t+1) =  DA*At1(t) + DX5*At2(t) + DX1*An1(t) + DX2*An2(t) + ep1(t+1)  
At2(t+1) = DX5*At1(t) +  DA*At2(t) + DX2*An1(t) + DX1*An2(t) + ep2(t+1)
An1(t+1) = DX3*At1(t) + DX4*At2(t) +  DN*An1(t) + DX6*An2(t) + ep3(t+1)  
An2(t+1) = DX4*At1(t) + DX3*At2(t) + DX6*An1(t) +  DN*An2(t) + ep4(t+1)
TT(t+1) = DZ*TT(t) + ep5(t+1)
M1(t+1) = DM*M1(t) + ep6(t+1)  
M2(t+1) = DM*M2(t) + ep7(t+1)
U1(t+1) = DU*U1(t) + ep8(t+1)  
U2(t+1) = DU*U2(t) + ep9(t+1)
G1(t+1) = DG*G1(t) + ep10(t+1)  
G2(t+1) = DG*G2(t) + ep11(t+1)
V1(t+1) = DV*V1(t) + ep12(t+1)  
V2(t+1) = DV*V2(t) + ep13(t+1)
S1(t+1) = DS*S1(t) + ep14(t+1)  
S2(t+1) = DS*S2(t) + ep15(t+1)
D1(t+1) = DD*D1(t) + ep16(t+1)  
D2(t+1) = DD*D2(t) + ep17(t+1)
B1(t+1) = DB*B1(t) + ep18(t+1)  
B2(t+1) = DB*B2(t) + ep19(t+1)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

pd1i(t) = pd1(t) - pd1(t-1)
pd2i(t) = pd2(t) - pd2(t-1)
px1i(t) = px1(t) - px1(t-1)
px2i(t) = px2(t) - px2(t-1)
pn1i(t) = pn1(t) - pn1(t-1)
pn2i(t) = pn2(t) - pn2(t-1)
wg1i(t) = wg1(t) - wg1(t-1)
wg2i(t) = wg2(t) - wg2(t-1)

(1-HB)*cx1(t) = c1(t) - HB*c1d(t-1) - ((1-HB)/RH)*B1(t)
(1-HB)*cx2(t) = c2(t) - HB*c2d(t-1) - ((1-HB)/RH)*B2(t)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Xd1(t) = (1-BT*GM)*xd1(t) + BT*GM*Et[Xd1(t+1)] - BT*GM*BK*pd1i(t)                              
Xd2(t) = (1-BT*GM)*xd2(t) + BT*GM*Et[Xd2(t+1)] - BT*GM*BK*pd2i(t)  

Xx1(t) = (1-BT*GM)*xx1(t) + BT*GM*Et[Xx1(t+1)] - BT*GM*BK*px1i(t)                               
Xx2(t) = (1-BT*GM)*xx2(t) + BT*GM*Et[Xx2(t+1)] - BT*GM*BK*px2i(t) 

Xn1(t) = (1-BT*GM)*xn1(t) + BT*GM*Et[Xn1(t+1)] - BT*GM*BK*pn1i(t)                               
Xn2(t) = (1-BT*GM)*xn2(t) + BT*GM*Et[Xn2(t+1)] - BT*GM*BK*pn2i(t)

Xw1(t) = (1-BT*GW)*w1(t) + BT*GW*Et[Xw1(t+1)] - BT*GW*BW*wg1i(t)                             
Xw2(t) = (1-BT*GW)*w2(t) + BT*GW*Et[Xw2(t+1)] - BT*GW*BW*wg2i(t)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

- (1-DC)*RH*cx1(t) + (1-DC)*RH*cx2(t) = - RH*Et[cx1(t+1)] + RH*Et[cx2(t+1)] - q(t) + Et[q(t+1)]

- RH*Et[cx1(t+1)] + (1-DC)*RH*cx1(t) + Et[re1(t+1)] = 0
- RH*Et[cx1(t+1)] + (1-DC)*RH*cx1(t) + Et[re2(t+1)] = 0 
- RH*Et[cx1(t+1)] + (1-DC)*RH*cx1(t) + Et[rb1(t+1)] = 0
- RH*Et[cx1(t+1)] + (1-DC)*RH*cx1(t) + Et[rb2(t+1)] = 0

re1(t) = (1-BT)*pr1(t) + BT*Z1(t) - Z1d(t-1)  
re2(t) = (1-BT)*pr2(t) + BT*Z2(t) - Z2d(t-1) 

rb1(t) = (1-BT*(1-1/BM))*( - p1(t)         ) + (BT*(1-1/BM))*Z3(t) - Z3d(t-1)
rb2(t) = (1-BT*(1-1/BM))*( - p2(t) + q(t)  ) + (BT*(1-1/BM))*Z4(t) - Z4d(t-1) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

ie1(t) = int1(t)
ie2(t) = int2(t)

int1(t) = DI*int1L(t-1) + DR*( DY*( y1(t) - yfx1(t) ) + (1.1+DT)*( CP*(1-FW)*pi1(t) + CP*FW*Et[pi1(t+1)] +  (1-CP)*(1-FW)*(py1(t) - py1d(t-1) ) +  (1-CP)*FW*(Et[py1(t+1)] - py1(t)) ) ) + M1(t)                         
int2(t) = DI*int2L(t-1) + DR*( DY*( y2(t) - yfx2(t) ) + (1.1+DT)*( CP*(1-FW)*pi2(t) + CP*FW*Et[pi2(t+1)] +  (1-CP)*(1-FW)*(py2(t) - py2d(t-1) ) +  (1-CP)*FW*(Et[py2(t+1)] - py2(t)) ) ) + M2(t)

ie1(t) = Et[rb1(t+1)] + Et[p1(t+1)] - p1(t)
ie2(t) = Et[rb2(t+1)] - Et[q(t+1)] + q(t) + Et[p2(t+1)] - p2(t)

pi1(t) = p1(t) - p1d(t-1)
pi2(t) = p2(t) - p2d(t-1)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

rkt1(t) = BT*BX*( MKt1(t) + qt1(t) - p1(t) ) + BT*(1-DP)*PS*(It1(t) + D1(t))
rkt2(t) = BT*BX*( MKt2(t) + qt2(t) - p2(t) ) + BT*(1-DP)*PS*(It2(t) + D2(t))

- (1-DC)*RH*cx1(t) + PS*(It1(t) + D1(t)) = - RH*Et[cx1(t+1)] + Et[rkt1(t+1)] + U1(t)
- (1-DC)*RH*cx2(t) + PS*(It2(t) + D2(t)) = - RH*Et[cx2(t+1)] + Et[rkt2(t+1)] + U2(t)

rkn1(t) = BT*BX*( MKn1(t) + qn1(t) - p1(t) ) + BT*(1-DP)*PS*(In1(t) + D1(t)) 
rkn2(t) = BT*BX*( MKn2(t) + qn2(t) - p2(t) ) + BT*(1-DP)*PS*(In2(t) + D2(t)) 

- (1-DC)*RH*cx1(t) + PS*(In1(t) + D1(t)) = - RH*Et[cx1(t+1)] + Et[rkn1(t+1)] + U1(t)
- (1-DC)*RH*cx2(t) + PS*(In2(t) + D2(t)) = - RH*Et[cx2(t+1)] + Et[rkn2(t+1)] + U2(t)

Iv1(t) = KP*In1(t) + (1-KP)*It1(t)
Iv2(t) = KP*In2(t) + (1-KP)*It2(t)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

PTB*prt1(t) = YTB*( yt1(t) + pa1(t) ) - WB*AF*YTB*( Lt1(t) + wg1(t) ) - ITB*( It1(t) + p1(t) ) - (1-KP)*TX*GY*( G1(t) + pg1(t) ) - PTB*p1(t) 
PTB*prt2(t) = YTB*( yt2(t) + pa2(t) ) - WB*AF*YTB*( Lt2(t) + wg2(t) ) - ITB*( It2(t) + p2(t) ) - (1-KP)*TX*GY*( G2(t) + pg2(t) ) - PTB*p2(t) + PTB*q(t) 

PNB*prn1(t) = YNB*( yn1(t) + pn1(t) ) - WB*AX*YNB*( Ln1(t) + wg1(t) ) - INB*( In1(t) + p1(t) ) - KP*TX*GY*( G1(t) + pg1(t) ) - PNB*p1(t)
PNB*prn2(t) = YNB*( yn2(t) + pn2(t) ) - WB*AX*YNB*( Ln2(t) + wg2(t) ) - INB*( In2(t) + p2(t) ) - KP*TX*GY*( G2(t) + pg2(t) ) - PNB*p2(t) + PNB*q(t)

(PTB+PNB)*pr1(t) = PTB*prt1(t) + PNB*prn1(t)   
(PTB+PNB)*pr2(t) = PTB*prt2(t) + PNB*prn2(t)  

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

xd1(t) = qt1(t) + V1(t)                                         
xd2(t) = qt2(t) + V2(t) 

xx1(t) = qt1(t) + V1(t) - (1-E1)*S(t)                                          
xx2(t) = qt2(t) + V2(t) + (1-E2)*S(t) 

MLt1(t) = (AF-1)*Lt1(t) + (1-AF)*Kt1(t-1) + (1-AF)*cpt1(t) + At1(t)
MLt2(t) = (AF-1)*Lt2(t) + (1-AF)*Kt2(t-1) + (1-AF)*cpt2(t) + At2(t)

MKt1(t) = AF*Lt1(t) - AF*Kt1(t-1) + (1-AF)*cpt1(t) + At1(t)  
MKt2(t) = AF*Lt2(t) - AF*Kt2(t-1) + (1-AF)*cpt2(t) + At2(t)

MLt1(t) = wg1(t) - qt1(t)
MLt2(t) = wg2(t) - qt2(t) 

yt1(t) = AF*Lt1(t) + (1-AF)*Kt1(t-1) + (1-AF)*cpt1(t) + At1(t)  
yt2(t) = AF*Lt2(t) + (1-AF)*Kt2(t-1) + (1-AF)*cpt2(t) + At2(t)

AF*Lt1(t) + (1-AF)*Kt1(t-1) - AF*cpt1(t) + At1(t) + qt1(t) - p1(t) = FZ*cpt1(t)
AF*Lt2(t) + (1-AF)*Kt2(t-1) - AF*cpt2(t) + At2(t) + qt2(t) - p2(t) = FZ*cpt2(t)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

xn1(t) = qn1(t) + V1(t)
xn2(t) = qn2(t) + V2(t)

MLn1(t) = (AX-1)*Ln1(t) + (1-AX)*Kn1(t-1) + (1-AX)*cpn1(t) + An1(t)
MLn2(t) = (AX-1)*Ln2(t) + (1-AX)*Kn2(t-1) + (1-AX)*cpn2(t) + An2(t)

MKn1(t) = AX*Ln1(t) - AX*Kn1(t-1) + (1-AX)*cpn1(t) + An1(t)  
MKn2(t) = AX*Ln2(t) - AX*Kn2(t-1) + (1-AX)*cpn2(t) + An2(t)

MLn1(t) = wg1(t) - qn1(t)
MLn2(t) = wg2(t) - qn2(t)

yn1(t) = AX*Ln1(t) + (1-AX)*Kn1(t-1) + (1-AX)*cpn1(t) + An1(t)  
yn2(t) = AX*Ln2(t) + (1-AX)*Kn2(t-1) + (1-AX)*cpn2(t) + An2(t)

AX*Ln1(t) + (1-AX)*Kn1(t-1) - AX*cpn1(t) + An1(t) + qn1(t) - p1(t) = FZ*cpn1(t)
AX*Ln2(t) + (1-AX)*Kn2(t-1) - AX*cpn2(t) + An2(t) + qn2(t) - p2(t) = FZ*cpn2(t)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

L1(t) = (1-KP)*Lt1(t) + KP*Ln1(t)
L2(t) = (1-KP)*Lt2(t) + KP*Ln2(t)

w1(t) = p1(t) + RH*cx1(t) + (MU-1)*L1(t) + S1(t)
w2(t) = p2(t) + RH*cx2(t) + (MU-1)*L2(t) + S2(t)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

yn1(t) = CB*c1(t) + ITB*It1(t) + INB*In1(t) + GY*G1(t) - (1-GY)*FI*( pn1(t) - p1(t) ) 
yn2(t) = CB*c2(t) + ITB*It2(t) + INB*In2(t) + GY*G2(t) - (1-GY)*FI*( pn2(t) - p2(t) ) 

dt1(t) = ( CB*c1(t) + ITB*It1(t) + INB*In1(t) )/(1-GY) - FI*( pt1(t) - p1(t) )
dt2(t) = ( CB*c2(t) + ITB*It2(t) + INB*In2(t) )/(1-GY) - FI*( pt2(t) - p2(t) )

dh1(t) = dt1(t) - TH*( pd1(t) - pt1(t) )
dh2(t) = dt2(t) - TH*( px1(t) - E1*S(t) - pt2(t) )                                       
df1(t) = dt1(t) - TH*( px2(t) + E2*S(t) - pt1(t) )  
df2(t) = dt2(t) - TH*( pd2(t) - pt2(t) )   

yt1(t) = (1-GY)*( ha1*dh1(t) + ha2*dh2(t) ) + TT(t) + GY*G1(t)                                                
yt2(t) = (1-GY)*( fa1*df1(t) + fa2*df2(t) ) - TT(t) + GY*G2(t)  

y1(t) = KP*(yn1(t) + pn1(t) - py1(t)) + (1-KP)*(yt1(t) + pa1(t) - py1(t)) 
y2(t) = KP*(yn2(t) + pn2(t) - py2(t)) + (1-KP)*(yt2(t) + pa2(t) - py2(t))

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

p1(t) = KP*pn1(t) + (1-KP)*pt1(t)                                           
p2(t) = KP*pn2(t) + (1-KP)*pt2(t)

pt1(t) = hs1*pd1(t) + fs1*( px2(t) + E2*S(t) )                                           
pt2(t) = hs2*( px1(t) - E1*S(t) ) + fs2*pd2(t) 

pg1(t) = KP*pn1(t) + (1-KP)*pd1(t)                                           
pg2(t) = KP*pn2(t) + (1-KP)*pd2(t)

pa1(t) = (1-GY)*( ha1*pd1(t) + ha2*( px1(t) + (1-E1)*S(t) ) ) + GY*pd1(t) 
pa2(t) = (1-GY)*( fa1*( px2(t) - (1-E2)*S(t) ) + fa2*pd2(t) ) + GY*pd2(t)

py1(t) = KP*pn1(t) + (1-KP)*pa1(t)
py2(t) = KP*pn2(t) + (1-KP)*pa2(t)


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

cd(t) = cx1(t) - cx2(t) - (1/RH)*q(t)

q(t) = S(t) + p2(t) - p1(t) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

yd(t) = CB*yy1(t) - CB*yy2(t) - CB*(1/RH)*( B1(t) - B2(t) ) + CB*((RH-1)/RH)*q(t)
yx(t) = yd(t) + BT*Et[yx(t+1)] 

CB*yy1(t) = YNB*(yn1(t) + pn1(t) - p1(t)) + YTB*(yt1(t) + pa1(t) - p1(t)) - GY*( G1(t) + pg1(t) - p1(t) ) - ITB*It1(t) - INB*In1(t)  
CB*yy2(t) = YNB*(yn2(t) + pn2(t) - p2(t)) + YTB*(yt2(t) + pa2(t) - p2(t)) - GY*( G2(t) + pg2(t) - p2(t) ) - ITB*It2(t) - INB*In2(t) + CB*q(t) 

rxb(t) = rb2(t) - rb1(t)
rxe(t) = re2(t) - re1(t)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

yfx1(t) = 0 
yfx2(t) = 0

$$$$$$$ pre-determined variables (other than shocks)

pd1(t-1)
pd2(t-1)
px1(t-1)
px2(t-1)
pn1(t-1)
pn2(t-1)
wg1(t-1)
wg2(t-1)

Kt1(t-1)
Kt2(t-1)
Kn1(t-1)
Kn2(t-1)

BX(t-1)

Z1d(t-1)
Z2d(t-1)
Z3d(t-1)
Z4d(t-1)

int1L(t-1)
int2L(t-1)

p1d(t-1)
p2d(t-1)

py1d(t-1)
py2d(t-1)

c1d(t-1)
c2d(t-1)

pd1id(t-1)
pd2id(t-1)
px1id(t-1)
px2id(t-1)
pn1id(t-1)
pn2id(t-1)
wg1id(t-1)
wg2id(t-1)

ys1(t-1)
ys2(t-1)

$$$$$$$ shocks

At1(t)
At2(t)
An1(t)
An2(t)
TT(t)
M1(t)
M2(t)
U1(t)
U2(t)
G1(t)
G2(t)
V1(t)
V2(t)
S1(t)
S2(t)
D1(t)
D2(t)
B1(t)
B2(t)

$$$$$$$ all other variables

pd1i(t)
pd2i(t)
px1i(t)
px2i(t)
pn1i(t)
pn2i(t)
wg1i(t)
wg2i(t)

Z1(t)
Z2(t)
Z3(t)
Z4(t)

Xd1(t)
Xd2(t) 
Xx1(t)
Xx2(t)
Xn1(t)
Xn2(t)
Xw1(t)
Xw2(t)

S(t)

c1(t)
c2(t)

cx1(t)
cx2(t)

re1(t) 
re2(t)

rb1(t) 
rb2(t)

ie1(t)
ie2(t)

int1(t)
int2(t)

It1(t)
It2(t)

In1(t)
In2(t)

Iv1(t)
Iv2(t)

rkt1(t)
rkt2(t)

rkn1(t)
rkn2(t)

cpn1(t)
cpn2(t)

cpt1(t)
cpt2(t)

prt1(t)
prt2(t)

prn1(t)
prn2(t)

pr1(t)
pr2(t)

xd1(t)  
xd2(t)

xx1(t)  
xx2(t)

qt1(t)
qt2(t)

MLt1(t)
MLt2(t)

MKt1(t)
MKt2(t)

Lt1(t)
Lt2(t)

xn1(t)
xn2(t)

qn1(t)
qn2(t)

Ln1(t)
Ln2(t)

MLn1(t)
MLn2(t)

MKn1(t)
MKn2(t)

L1(t)
L2(t)

w1(t)
w2(t) 

pi1(t)
pi2(t)

dt1(t)
dt2(t)

dh1(t)
dh2(t)  
df1(t)
df2(t) 

yn1(t)  
yn2(t)

yt1(t)  
yt2(t)

y1(t)
y2(t)

p1(t)  
p2(t)  

pa1(t)
pa2(t)

pg1(t)
pg2(t)

pt1(t)
pt2(t)

py1(t)
py2(t)

cd(t)

q(t)

yd(t)
yy1(t)
yy2(t)
 
yx(t)

rxb(t)
rxe(t)

yfx1(t)
yfx2(t)

$$$$$$$ iid shocks

ep1(t+1)
ep2(t+1)
ep3(t+1)
ep4(t+1)
ep5(t+1)
ep6(t+1)
ep7(t+1)
ep8(t+1)
ep9(t+1)
ep10(t+1)
ep11(t+1)
ep12(t+1)
ep13(t+1)
ep14(t+1)
ep15(t+1)
ep16(t+1)
ep17(t+1)
ep18(t+1)
ep19(t+1)

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

DC = 0.00001

GM = 0.75 
GW = 0.75             
BT = 0.99
TH = 1.5
FF = 10.0
MU = 2.5
NA = 0.5
RH = 1.5        
VV = 2-2*0.58
ET = 1.0
AF = 0.67
GY = 0.2        
KP = 0.4        
FI = 0.45       
TX = 0.15 

FZ = 0.2

HB = 0.0
BK = 0.5
BW = 0.5

BM = 10000

DP = 0.025     
PS = 0.25

DT = 3.9
DI = 0.85       
DY = 0.1       

CP = 1
FW = 0
CY = 1.0

DA = 0.9          
DN = 0.9          
DX1 = 0.6         
DX2 = 0.6       
DX3 = 0
DX4 = 0
DX5 = 0
DX6 = 0

DZ = 0.90
DM = 0.00
DG = 0.90
DV = 0.00
DS = 0.90
DB = 0.90
DU = 0.00
DD = 0.90

SA = 0.009        
SN = 0.005

ST = 0.00          
SM = 0.0012           
SG = 0.003          
SV = 0.0015          
SS = 0.025            
SB = 0.000
SU = 0.006 
SD = 0.001
 
CX1 = 0     
CX2 = 0        
CX3 = 0        
CX4 = 0        

AX = AF

E1 = ET
E2 = ET

DR = 1-DI

hs1 = 1-(1-NA)*VV
hs2 = NA*VV
fs1 = (1-NA)*VV
fs2 = 1-NA*VV

ha1 = 1-(1-NA)*VV
ha2 = (1-NA)*VV
fa1 = NA*VV
fa2 = 1-NA*VV

YB = 1

YTB = (1-KP)*YB
YNB = KP*YB

WB = (FF-1)/FF
BX = (1/BT)-1+DP

KTB = (1-AF)*WB*YTB/BX
ITB = DP*KTB

KNB = (1-AX)*WB*YNB/BX
INB = DP*KNB

CB = YB-ITB-INB-GY

PTB = (1-WB*AF)*YTB - ITB - (1-KP)*TX*GY
PNB = (1-WB*AX)*YNB - INB - KP*TX*GY

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

SA*SA     CX1*SA*SA CX3*SA*SN CX4*SA*SN 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
CX1*SA*SA SA*SA     CX4*SA*SN CX3*SA*SN 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
CX3*SA*SN CX4*SA*SN SN*SN     CX2*SN*SN 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
CX4*SA*SN CX3*SA*SN CX2*SN*SN SN*SN     0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       ST^2 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  SM^2 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  SM^2 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  SU^2 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  SU^2 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  0.0  SG^2 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  0.0  0.0  SG^2 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  0.0  0.0  0.0  SV^2 0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  SV^2 0.0  0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  SS^2 0.0  0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  SS^2 0.0  0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  SD^2 0.0  0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  SD^2 0.0  0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  SB^2 0.0
0.0       0.0       0.0       0.0       0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  SB^2

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$