* NLP written by GAMS Convert at 10/06/06 11:34:14 * * 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 * 32 32 0 0 0 0 0 0 * FX 0 0 0 0 0 0 0 0 * * Nonzero counts * Total const NL DLL * 32 1 31 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,objvar; Equations e1; e1.. - (sqr(0.01*(0.000826446280991736*x25*x2*exp(18.5751 - 3632.649/(239.2 + x1)) - x5*(693.37 + x25) + 693.37*x2)) + sqr(0.000869565217391304*exp( 18.5751 - 3632.649/(239.2 + x19))*x20 - x23) + sqr(0.01*(- 0.000826446280991736*x25*x2*exp(18.5751 - 3632.649/(239.2 + x1)) - x8*( 693.37 + x26) + x5*(693.37 + x25) + 0.000833333333333333*x26*x5*exp( 18.5751 - 3632.649/(239.2 + x4)))) + sqr((-4.5125) + 0.01*(- 0.000833333333333333*x26*x5*exp(18.5751 - 3632.649/(239.2 + x4)) - x11*(- 442.13 + x27) + x8*(693.37 + x26) + 0.000840336134453782*x27*x8*exp( 18.5751 - 3632.649/(239.2 + x7)))) + sqr(0.01*(-0.000840336134453782*x27* x8*exp(18.5751 - 3632.649/(239.2 + x7)) - x14*(-442.13 + x28) + x11*(- 442.13 + x27) + 0.000847457627118644*x28*x11*exp(18.5751 - 3632.649/(239.2 + x10)))) + sqr(0.01*(-0.000847457627118644*x28*x11*exp(18.5751 - 3632.649/(239.2 + x10)) - x17*(-442.13 + x29) + x14*(-442.13 + x28) + 0.000854700854700855*x29*x14*exp(18.5751 - 3632.649/(239.2 + x13)))) + sqr(0.01*(-0.000854700854700855*x29*x14*exp(18.5751 - 3632.649/(239.2 + x13)) - x20*(-442.13 + x30) + x17*(-442.13 + x29) + 0.000862068965517241* x30*x17*exp(18.5751 - 3632.649/(239.2 + x16)))) + sqr(0.01*(- 0.000862068965517241*x30*x17*exp(18.5751 - 3632.649/(239.2 + x16)) - x23*( -442.13 + x31) + x20*(-442.13 + x30) + 0.000869565217391304*x31*x20*exp( 18.5751 - 3632.649/(239.2 + x19)))) + sqr(0.01*(0.000826446280991736*x25* x3*exp(18.3443 - 3841.2203/(228 + x1)) - x6*(693.37 + x25) + 693.37*x3)) + sqr(0.000869565217391304*exp(18.3443 - 3841.2203/(228 + x19))*x21 - x24 ) + sqr(0.01*(-0.000826446280991736*x25*x3*exp(18.3443 - 3841.2203/(228 + x1)) - x9*(693.37 + x26) + x6*(693.37 + x25) + 0.000833333333333333*x26*x6 *exp(18.3443 - 3841.2203/(228 + x4)))) + sqr((-6.8425) + 0.01*(- 0.000833333333333333*x26*x6*exp(18.3443 - 3841.2203/(228 + x4)) - x12*(- 442.13 + x27) + x9*(693.37 + x26) + 0.000840336134453782*x27*x9*exp( 18.3443 - 3841.2203/(228 + x7)))) + sqr(0.01*(-0.000840336134453782*x27*x9 *exp(18.3443 - 3841.2203/(228 + x7)) - x15*(-442.13 + x28) + x12*(-442.13 + x27) + 0.000847457627118644*x28*x12*exp(18.3443 - 3841.2203/(228 + x10) ))) + sqr(0.01*(-0.000847457627118644*x28*x12*exp(18.3443 - 3841.2203/(228 + x10)) - x18*(-442.13 + x29) + x15*(-442.13 + x28) + 0.000854700854700855*x29*x15*exp(18.3443 - 3841.2203/(228 + x13)))) + sqr( 0.01*(-0.000854700854700855*x29*x15*exp(18.3443 - 3841.2203/(228 + x13)) - x21*(-442.13 + x30) + x18*(-442.13 + x29) + 0.000862068965517241*x30* x18*exp(18.3443 - 3841.2203/(228 + x16)))) + sqr(0.01*(- 0.000862068965517241*x30*x18*exp(18.3443 - 3841.2203/(228 + x16)) - x24*(- 442.13 + x31) + x21*(-442.13 + x30) + 0.000869565217391304*x31*x21*exp( 18.3443 - 3841.2203/(228 + x19)))) + sqr((-1) + 0.000826446280991736*exp( 18.5751 - 3632.649/(239.2 + x1))*x2 + 0.000826446280991736*exp(18.3443 - 3841.2203/(228 + x1))*x3) + sqr((-1) + 0.000833333333333333*exp(18.5751 - 3632.649/(239.2 + x4))*x5 + 0.000833333333333333*exp(18.3443 - 3841.2203/( 228 + x4))*x6) + sqr((-1) + 0.000840336134453782*exp(18.5751 - 3632.649/( 239.2 + x7))*x8 + 0.000840336134453782*exp(18.3443 - 3841.2203/(228 + x7)) *x9) + sqr((-1) + 0.000847457627118644*exp(18.5751 - 3632.649/(239.2 + x10 ))*x11 + 0.000847457627118644*exp(18.3443 - 3841.2203/(228 + x10))*x12) + sqr((-1) + 0.000854700854700855*exp(18.5751 - 3632.649/(239.2 + x13))*x14 + 0.000854700854700855*exp(18.3443 - 3841.2203/(228 + x13))*x15) + sqr((- 1) + 0.000862068965517241*exp(18.5751 - 3632.649/(239.2 + x16))*x17 + 0.000862068965517241*exp(18.3443 - 3841.2203/(228 + x16))*x18) + sqr((-1) + 0.000869565217391304*exp(18.5751 - 3632.649/(239.2 + x19))*x20 + 0.000869565217391304*exp(18.3443 - 3841.2203/(228 + x19))*x21) + sqr((-1) + 0.00087719298245614*exp(18.5751 - 3632.649/(239.2 + x22))*x23 + 0.00087719298245614*exp(18.3443 - 3841.2203/(228 + x22))*x24) + sqr((- 83.862) + 1e-5*(0.000826446280991736*x25*x2*exp(18.5751 - 3632.649/(239.2 + x1))*(9566.67 + 0.0422*x1*x1 - 1.59*x1) + 693.37*(0.0422*x1*x1 + 15.97* x1)*x2 - x5*(693.37 + x25)*(0.0422*x4*x4 + 15.97*x4) + 0.000826446280991736*x25*x3*exp(18.3443 - 3841.2203/(228 + x1))*(10834.67 + 8.74*x1) + 12549.997*x3*x1 - 18.1*x6*(693.37 + x25)*x4)) + sqr(1e-5*( 0.000833333333333333*x26*x5*exp(18.5751 - 3632.649/(239.2 + x4))*(9566.67 + 0.0422*x4*x4 - 1.59*x4) + x5*(693.37 + x25)*(0.0422*x4*x4 + 15.97*x4) - 0.000826446280991736*x25*x2*exp(18.5751 - 3632.649/(239.2 + x1))*( 9566.67 + 0.0422*x1*x1 - 1.59*x1) - x8*(693.37 + x26)*(0.0422*x7*x7 + 15.97*x7) + 0.000833333333333333*x26*x6*exp(18.3443 - 3841.2203/(228 + x4) )*(10834.67 + 8.74*x4) + 18.1*x6*(693.37 + x25)*x4 - 0.000826446280991736* x25*x3*exp(18.3443 - 3841.2203/(228 + x1))*(10834.67 + 8.74*x1) - 18.1*x9* (693.37 + x26)*x7)) + sqr((-18.9447111025) + 1e-5*(0.000840336134453782* x27*x8*exp(18.5751 - 3632.649/(239.2 + x7))*(9566.67 + 0.0422*x7*x7 - 1.59 *x7) + x8*(693.37 + x26)*(0.0422*x7*x7 + 15.97*x7) - 0.000833333333333333* x26*x5*exp(18.5751 - 3632.649/(239.2 + x4))*(9566.67 + 0.0422*x4*x4 - 1.59 *x4) - x11*(-442.13 + x27)*(0.0422*x10*x10 + 15.97*x10) + 0.000840336134453782*x27*x9*exp(18.3443 - 3841.2203/(228 + x7))*(10834.67 + 8.74*x7) + 18.1*x9*(693.37 + x26)*x7 - 0.000833333333333333*x26*x6*exp( 18.3443 - 3841.2203/(228 + x4))*(10834.67 + 8.74*x4) - 18.1*x12*(-442.13 + x27)*x10)) + sqr(1e-5*(0.000847457627118644*x28*x11*exp(18.5751 - 3632.649/(239.2 + x10))*(9566.67 + 0.0422*x10*x10 - 1.59*x10) + x11*(- 442.13 + x27)*(0.0422*x10*x10 + 15.97*x10) - 0.000840336134453782*x27*x8* exp(18.5751 - 3632.649/(239.2 + x7))*(9566.67 + 0.0422*x7*x7 - 1.59*x7) - x14*(-442.13 + x28)*(0.0422*x13*x13 + 15.97*x13) + 0.000847457627118644* x28*x12*exp(18.3443 - 3841.2203/(228 + x10))*(10834.67 + 8.74*x10) + 18.1* x12*(-442.13 + x27)*x10 - 0.000840336134453782*x27*x9*exp(18.3443 - 3841.2203/(228 + x7))*(10834.67 + 8.74*x7) - 18.1*x15*(-442.13 + x28)*x13) ) + sqr(1e-5*(0.000854700854700855*x29*x14*exp(18.5751 - 3632.649/(239.2 + x13))*(9566.67 + 0.0422*x13*x13 - 1.59*x13) + x14*(-442.13 + x28)*( 0.0422*x13*x13 + 15.97*x13) - 0.000847457627118644*x28*x11*exp(18.5751 - 3632.649/(239.2 + x10))*(9566.67 + 0.0422*x10*x10 - 1.59*x10) - x17*(- 442.13 + x29)*(0.0422*x16*x16 + 15.97*x16) + 0.000854700854700855*x29*x15* exp(18.3443 - 3841.2203/(228 + x13))*(10834.67 + 8.74*x13) + 18.1*x15*(- 442.13 + x28)*x13 - 0.000847457627118644*x28*x12*exp(18.3443 - 3841.2203/( 228 + x10))*(10834.67 + 8.74*x10) - 18.1*x18*(-442.13 + x29)*x16)) + sqr( 1e-5*(0.000862068965517241*x30*x17*exp(18.5751 - 3632.649/(239.2 + x16))*( 9566.67 + 0.0422*x16*x16 - 1.59*x16) + x17*(-442.13 + x29)*(0.0422*x16*x16 + 15.97*x16) - 0.000854700854700855*x29*x14*exp(18.5751 - 3632.649/(239.2 + x13))*(9566.67 + 0.0422*x13*x13 - 1.59*x13) - x20*(-442.13 + x30)*( 0.0422*x19*x19 + 15.97*x19) + 0.000862068965517241*x30*x18*exp(18.3443 - 3841.2203/(228 + x16))*(10834.67 + 8.74*x16) + 18.1*x18*(-442.13 + x29)* x16 - 0.000854700854700855*x29*x15*exp(18.3443 - 3841.2203/(228 + x13))*( 10834.67 + 8.74*x13) - 18.1*x21*(-442.13 + x30)*x19)) + sqr(1e-5*( 0.000869565217391304*x31*x20*exp(18.5751 - 3632.649/(239.2 + x19))*( 9566.67 + 0.0422*x19*x19 - 1.59*x19) + x20*(-442.13 + x30)*(0.0422*x19*x19 + 15.97*x19) - 0.000862068965517241*x30*x17*exp(18.5751 - 3632.649/(239.2 + x16))*(9566.67 + 0.0422*x16*x16 - 1.59*x16) - x23*(-442.13 + x31)*( 0.0422*x22*x22 + 15.97*x22) + 0.000869565217391304*x31*x21*exp(18.3443 - 3841.2203/(228 + x19))*(10834.67 + 8.74*x19) + 18.1*x21*(-442.13 + x30)* x19 - 0.000862068965517241*x30*x18*exp(18.3443 - 3841.2203/(228 + x16))*( 10834.67 + 8.74*x16) - 18.1*x24*(-442.13 + x31)*x22))) + objvar =E= 0; * set non default bounds * set non default levels x1.l = 120; x2.l = 0.09203; x3.l = 0.908; x4.l = 110; x5.l = 0.1819; x6.l = 0.8181; x7.l = 100; x8.l = 0.284; x9.l = 0.716; x10.l = 88; x11.l = 0.3051; x12.l = 0.6949; x13.l = 86; x14.l = 0.3566; x15.l = 0.6434; x16.l = 84; x17.l = 0.468; x18.l = 0.532; x19.l = 80; x20.l = 0.6579; x21.l = 0.3421; x22.l = 76; x23.l = 0.8763; x24.l = 0.1237; x25.l = 886.37; x26.l = 910.01; x27.l = 922.52; x28.l = 926.46; x29.l = 935.56; x30.l = 952.83; x31.l = 975.73; * 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;