MODEL:
sets:
numb/1..20/:y,z,Y0,Z0,I0,S0,P0,Ir,Sr,Pr;
endsets
data:
P0=1285,682,56,87,46,78,218,52,81,83,42,41,72,128,220,78,90,32,58,121;
Pr=0.00204,0.00303,0.00915,0.0059,0.00587,0.00612,0.00541,0.0045,0.00369,0.00661,0.0080,0.0061,0.00601,0.00692,0.00512,0.00656,0.00661,0.00644,0.0046,0.0059;
I0=737,739,193,268,480,256,464,189,721,110,36,97,104,67,310,72,114,106,83,211;
Ir=11.1,11.7,10.0,12.5,9.8,7.6,10.2,10.9,10.0,6.9,9.6,8.7,10.3,8.0,12.9,11.1,8.8,10.0,9.0,10.4;
S0=1.16,0.83,0.30,0.23,0.22,0.2,0.44,0.15,0.22,0.16,0.18,0.14,0.22,0.18,0.53,0.17,0.18,0.12,0.15,0.13;
Sr=13.2,12.2,10.0,12.1,8.6,7.6,11.8,12.0,11.4,11.7,12.3,13.5,8.9,9.4,10.2,9.2,8.4,11.3,10.3,8.8;
Z0=160,140,180,360,315,318,235,315,320,310,320,352,280,310,220,320,310,340,280,320;
Y0=143,72,102,96,110,120,86,131,126,186,210,170,205,180,88,210,189,210,200,180;
enddata
[OBJ]max=@sum(numb(j):(((y(j)(1+Ir(j))^10)/Y0(j))+(((Z(j)(1+Sr(j))^10)(1+Pr(j))^10/Z0(j))));
@sum(numb(j):y(j))=0.38110e8;
@sum(numb(j):z(j))=0.22*110e8;
@for(numb(j):y(j)/(Y0(j)*I0(j)10^4)>=0.50.63169;);
@for(numb(j):y(j)/(Y0(j)*I0(j)10^4)<=1.5*0.63169;);
@for(numb(j):z(j)/(Z0(j)*S0(j)*P0(j)*10^4)>=0.50.56856;);
@for(numb(j):z(j)/(Z0(j)*S0(j)*P0(j)*10^4)<=1.5*0.56856;);
@for(numb(j):y(j)>=0;);
@for(numb(j):z(j)>=0;);
END