# PAVIANI: Himmelblau 20
# written in AMPL by Yu-Ju Kuo and Hans D. Mittelmann, 6/2001
param N;
param NH;
param NG;
param F := 0.7302D0*530.D0*14.7D0/40.D0;
param B{1..12};
param C{1..12};
param D{1..12};
param A{1..24};
param E{1..6};
var X{1..N}:=0;
var S1= sum{j in 13..24}X[j]/B[j-12];
var S2= sum{j in 13..24}X[j-12]/B[j-12];
var S3 = sum{j in 1..12} X[j]/D[j];
minimize obj:100*sum{i in 1..N}A[i]*X[i];
s.t. hxi{i in 1..12}: X[i+12]/(B[i]*S1)-C[i]*X[i]/(40*B[i]*S2)=0;
s.t. hx13 : sum{i in 1..N}X[i]-1=0;
s.t. hx14 : S3+F*S1-1.671=0;
s.t. gxi{i in 1..N} : X[i]*100>=0;
s.t. gx25{i in 1..3}: -(X[i]+X[i+12])/sum{j in 1..N}X[j]+E[i]>=0;
s.t. gx29{i in 4..6}: -(X[i+3]+X[i+15])/sum{j in 1..N}X[j]+E[i]>=0;
s.t. gx32 : 100*sum{i in 1..12}X[i]>=0;
s.t. gx33 : 100*sum{i in 13..N}X[i]>=0;



data;
param N:=24;
param NH:=14;
param NG:=32;
param A:=
1 0.0693
2 0.0577
3 0.05
4 0.2
5 0.26
6 0.55
7 0.06
8 0.1
9 0.12
10 0.18
11 0.1
12 0.09
13 0.0693
14 0.0577
15 0.05
16 0.2
17 0.26
18 0.55
19 0.06
20 0.1
21 0.12
22 0.18
23 0.1
24 0.09;
param B:=
1 44.094
2 58.12
3 58.12
4 137.4
5 120.9
6 170.9
7 62.501
8 84.94
9 133.425
10 82.507
11 46.07
12 60.097;
param C:=
1 123.7
2 31.7
3 45.7
4 14.7
5 84.7
6 27.7
7 49.7
8 7.1
9 2.1
10 17.7
11 0.85
12 0.64;
param D:=
1 31.244
2 36.12
3 34.784
4 92.7
5 82.7
6 91.6
7 56.708
8  82.7
9 80.8
10 64.517
11 49.4
12 49.1;
param E:=
1 0.1
2 0.3
3 0.4
4 0.3
5 0.6
6 0.3;