gridneta.gms
<![CDATA[
* NLP written by GAMS Convert at 06/26/06 16:29:35
*
* Equation counts
* Total E G L N X C
* 1 1 0 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 73 73 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 73 1 72 0
*
* Solve m using NLP minimizing objvar;
Variables x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19
,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36
,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47,x48,x49,x50,x51,x52,x53
,x54,x55,x56,x57,x58,x59,x60,x61,x62,x63,x64,x65,x66,x67,x68,x69,x70
,x71,x72,objvar;
Positive Variables x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17
,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34
,x35,x36,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47,x48,x49,x50,x51
,x52,x53,x54,x55,x56,x57,x58,x59,x60,x61,x62,x63,x64,x65,x66,x67,x68
,x69,x70,x71,x72;
Equations e1;
e1.. - (0.01*(x1*x1 + x2*x2 + x3*x3 + 2*x4*x4 + 2*x5*x5 + x6*x6 + x7*x7 + 4*x8
*x8 + 4*x9*x9 + x10*x10 + x12*x12 + x13*x13 + 2*x14*x14 + 2*x15*x15 + 2*
x16*x16 + x17*x17 + 2*x18*x18 + 4*x19*x19 + 4*x20*x20 + x21*x21 + 4*x22*
x22 + x24*x24 + x25*x25 + x26*x26 + 2*x27*x27 + 2*x28*x28 + 2*x29*x29 + 2*
x30*x30 + x31*x31 + 4*x32*x32 + 2*x33*x33 + 4*x34*x34 + x36*x36 + 2*x37*
x37 + 4*x38*x38 + 2*x39*x39 + 2*x40*x40 + 2*x41*x41 + 2*x42*x42 + 2*x43*
x43 + 4*x44*x44 + 2*x45*x45 + 4*x46*x46 + x48*x48 + x49*x49 + x50*x50 + 4*
x51*x51 + 2*x52*x52 + 4*x53*x53 + x54*x54 + 4*x55*x55 + 4*x56*x56 + 4*x57*
x57 + 4*x58*x58 + x60*x60 + x61*x61 + x63*x63 + x65*x65 + x67*x67 + x69*
x69) + 0.01*(sqrt(1 + sqr(x1) + sqr(x1 - x14)) + sqrt(1 + sqr(x2) + sqr(x2
- 6*x14)) + sqrt(1 + sqr(x3) + sqr(x3 - x16)) + sqrt(1 + sqr(x4) + sqr(x4
- 6*x16)) + sqrt(1 + sqr(x5) + sqr(x5 - x18)) + sqrt(1 + sqr(x6) + sqr(x6
- 6*x18)) + sqrt(1 + sqr(x7) + sqr(x7 - x20)) + sqrt(1 + sqr(x8) + sqr(x8
- 6*x20)) + sqrt(1 + sqr(x9) + sqr(x9 - x23)) + sqrt(1 + sqr(x10) + sqr(
x10 - 6*x22)) + sqrt(1 + sqr(x12) + sqr(x12 - 6*x24)) + sqrt(1 + sqr(x13)
+ sqr(x13 - x26)) + sqrt(1 + sqr(x14) + sqr(x14 - 6*x26)) + sqrt(1 + sqr(
x15) + sqr(x15 - x28)) + sqrt(1 + sqr(x16) + sqr(x16 - 6*x28)) + sqrt(1 +
sqr(x17) + sqr(x17 - x30)) + sqrt(1 + sqr(x18) + sqr(x18 - 6*x30)) + sqrt(
1 + sqr(x19) + sqr(x19 - x32)) + sqrt(1 + sqr(x20) + sqr(x20 - 6*x32)) +
sqrt(1 + sqr(x21) + sqr(x21 - x35)) + sqrt(1 + sqr(x22) + sqr(x22 - 6*x34)
) + sqrt(1 + sqr(x24) + sqr(x24 - 6*x36)) + sqrt(1 + sqr(x25) + sqr(x25 -
x38)) + sqrt(1 + sqr(x26) + sqr(x26 - 6*x38)) + sqrt(1 + sqr(x27) + sqr(
x27 - x40)) + sqrt(1 + sqr(x28) + sqr(x28 - 6*x40)) + sqrt(1 + sqr(x29) +
sqr(x29 - x42)) + sqrt(1 + sqr(x30) + sqr(x30 - 6*x42)) + sqrt(1 + sqr(x31
) + sqr(x31 - x44)) + sqrt(1 + sqr(x32) + sqr(x32 - 6*x44)) + sqrt(1 +
sqr(x33) + sqr(x33 - x47)) + sqrt(1 + sqr(x34) + sqr(x34 - 6*x46)) + sqrt(
1 + sqr(x36) + sqr(x36 - 6*x48)) + sqrt(1 + sqr(x37) + sqr(x37 - x50)) +
sqrt(1 + sqr(x38) + sqr(x38 - x62)) + sqrt(1 + sqr(x39) + sqr(x39 - x52))
+ sqrt(1 + sqr(x40) + sqr(x40 - x64)) + sqrt(1 + sqr(x41) + sqr(x41 - x54
)) + sqrt(1 + sqr(x42) + sqr(x42 - x66)) + sqrt(1 + sqr(x43) + sqr(x43 -
x56)) + sqrt(1 + sqr(x44) + sqr(x44 - x68)) + sqrt(1 + sqr(x45) + sqr(x45
- x59)) + sqrt(1 + sqr(x46) + sqr(x46 - x70)) + sqrt(1 + sqr(x48) + sqr(
x48 - x72)) + sqrt(1 + sqr(x49) + sqr(x49 - x62)) + sqrt(1 + sqr(x50) +
sqr(x50 - x3)) + sqrt(1 + sqr(x51) + sqr(x51 - x64)) + sqrt(1 + sqr(x52)
+ sqr(x52 - x5)) + sqrt(1 + sqr(x53) + sqr(x53 - x66)) + sqrt(1 + sqr(x54
) + sqr(x54 - x7)) + sqrt(1 + sqr(x55) + sqr(x55 - x68)) + sqrt(1 + sqr(
x56) + sqr(x56 - x9)) + sqrt(1 + sqr(x57) + sqr(x57 - x71)) + sqrt(1 +
sqr(x58) + sqr(x58 - x11)) + sqrt(1 + sqr(x60) + sqr(x60)) + sqrt(1 + sqr(
x61) + sqr(x61)) + sqrt(1 + sqr(x63) + sqr(x63)) + sqrt(1 + sqr(x65) +
sqr(x65)) + sqrt(1 + sqr(x67) + sqr(x67)) + sqrt(1 + sqr(x69) + sqr(x69)))
+ 0.000833333333333333*POWER(10 - x1 + x2 + x3 - x4 - x5 + x6 + x7 - x8
- x9 + x10 - x13 + x14 + x15 - x16 - x17 + x18 + x19 - x20 - x21 + x22 -
x25 + x26 + x27 - x28 - x29 + x30 + x31 - x32 - x33 + x34 - x37 + x38 +
x39 - x40 - x41 + x42 + x43 - x44 - x45 + x46 - x49 + x50 + x51 - x52 -
x53 + x54 + x55 - x56 - x57 + x58,4)) + objvar =E= 0;
* set non default bounds
* set non default levels
* set non default marginals
Model m / all /;
m.limrow=0; m.limcol=0;
$if NOT '%gams.u1%' == '' $include '%gams.u1%'
Solve m using NLP minimizing objvar;
]]></plaintext></body>