model:
sets:

S1/1..7/:LS1,XLS1;
S2/1..7/:CS2,FCS2,BCS2,ECS2,XCS2,YCS2,pCS2,aCS2,qCS2,FVS2,BVS2;
S3/1..7/:KS3,FKS3,BKS3,EKS3,XKS3,YKS3,pKS3,aKS3,qKS3,FGS3,BGS3;
linka(S1,S2):DLS1CS2,SLS1CS2,XLS1CS2;
linkb(S2,S3):DCS2KS3,SCS2KS3,XCS2KS3;

ENDSETS

data:
XLS1=320 152.5 659.5 180 1041.4 320 192;
DLS1CS2=0        54.1        46.4        45.4        40.3        65.2        26.8
54.1        0        13.1        50.8        51.8        52.0        59.3
46.4        13.1        0        38.5        40.0        33.3        47.2
45.4        50.8        38.5        0        17.4        16.0        36.1
40.3        51.8        40.0        17.4        0        30.4        21.2
65.2        52.0        33.3        16.0        30.4        0        42.9
  26.8      59.3   47.2        36.1        21.2        42.9   0;
SLS1CS2=35;

DCS2KS3=0        54.1        46.4        45.4        40.3        65.2        26.8
54.1        0        13.1        50.8        51.8        52.0        59.3
46.4        13.1        0        38.5        40.0        33.3        47.2
45.4        50.8        38.5        0        17.4        16.0        36.1
40.3        51.8        40.0        17.4        0        30.4        21.2
65.2        52.0        33.3        16.0        30.4        0        42.9
  26.8      59.3   47.2        36.1        21.2        42.9   0;
SCS2KS3=45;

FCS2=5000000;
BCS2=450;
ECS2=40;
pCS2=1172        1105        1494        1531        1218        1297        837;
aCS2=0.2;
qCS2=0.1;
FVS2=100000;
BVS2=30;

FKS3=48000000;
BKS3=1200;
EKS3=1089;
pKS3=1172        1105        1494        1531        1218        1297        837;
aKS3=1;
qKS3=0.2;
FGS3=300000;
BGS3=80;

ENDDATA


Min=Z1+Z2+Z3+R;

Z1=@sum(S2:YCS2*(FCS2+FVS2))+@sum(S3:YKS3*(FKS3+FGS3));
Z2=@sum(linka:DLS1CS2*ELS1CS2*XLS1CS2)+@sum(linkb:DCS2KS3*ECS2KS3*XCS2KS3);
Z3=@sum(S2:YCS2*XCS2*ECS2)+@sum(S3:YKS3*XKS3*EKS3);
R=@sum(S2:@PI()*pCS2*pCS2*aCS2*qCS2)+@sum(S3:@PI()*pKS3*pKS3*aKS3*qKS3);

@FOR(S1(I1):
    XLS1(I1)=@sum(S2(J1):XLS1CS2(I1,J1)));
@FOR(S2(J1):
    XCS2(J1)=@sum(S1(I1):XLS1CS2(I1,J1)));
@FOR(S2(J1):
    XCS2(J1)=@sum(S3(K1):XCS2KS3(J1,K1)));
@FOR(S3(K1):
    XKS3(K1)=@sum(S2(J1):XCS2KS3(J1,K1)));

@FOR(S2(J1):
    XCS2(J1)<=BCS2(J1)*YCS2(J1));
@FOR(S3(H1):
    XKS3(H1)<=BKS3(H1)*YKS3(H1));
@FOR(LINKa(I1,J1):
     XLS1CS2<=BVS2);
@FOR(LINKb(J1,H1):
     XCS2KS3<=BGS3);


@sum(S2(J1):
    YCS2(J1))>=1;
@sum(S2(J1):
    YCS2(J1))<=4;
@sum(S3(H1):
    YKS3(H1))=1;


@FOR(S2(J1):
    @BIN(YCS2(J1)));
@FOR(S3(H1):
    @BIN(YKS3(H1)));

END